r - 

F)D-ni41  960 

UNCLflSSIFIEl 

PHASE  CONJUGATE  OPTICS(U) 
CA  G  J  DUNNING  ET  AL  PEG 
F19628-80-C-0185 

HUGHES  RESEARCH  LABS  HALIBU 
84  RADC-TR-84-18 

F/G  20/6 

1/ 

NL 

■1 

w 

a 

□ 

-- 

■ 

m 

— 

Hi 

■■ 

AD-A141  960 


RADC-TR-84-18 
Final  Tochnlcal  Report 
February  1984 


PHASE  CONJUGATE  OPTICS 


Hughes  Research  Laboratories 


G.  J.  Dunning 
M.  B.  Klein 
R.  C.  Lind 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBOTIOIII  URLIMITEO  _  .  _ 

DTIC 


Selects 

JUN  1  1  1984 


ROME  AIR  DEVELOPMENT  CENTER 
Air  Force  Systems  Command 
Griff iss  Air  Force  Base,  NY  13441 


£ 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Office  (PA)  and 
Is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At  NTIS 
it  will  be  releasable  to  the  general  public,  including  foreign  nations. 


RADC  TR-84-18  has  been  reviewed  and  is  approved  for  publication. 


APPROVED: 


CHARLES  L.  WOODS 
Project  Engineer 


HAROLD  ROTH,  Director 
Solid  State  Sciences  Division 


FOR  THE  COMMANDER: 


JOHN  A.  RITZ 

Acting  Chief,  Plans  Office 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  (  ESO  )  Hanscom  AFB  MA  01731.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  recjuircs  that  it  be  returned. 


SKCuntTV  CLASUriCATION  or  This  PAOK  flWiMl  Omt*  enftmd) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  REPORT  number  U.  COVT  ACCESSION  NO.I  1  RECIPlEs'^'S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


1.  report  number 

RADC-TR-84-18 


4.  title  (Mti  Submit) 

PHASE  CONJUGATE  OPTICS 


S.  PERFORMING  OROANIZATION  NAME  AnO  ADDRESS 

Hughes  Research  Laboratories 
3011  Malibu  Canyon  Road 
Malibu  CA  90265 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Rome  Air  Development  Center  (ESO) 
Hanscom  AFB  MA  01731 


5  type  of  REPORT  4  PERIOD  COWERED 

Final  Technical  Report 
1  Aug  80  -  31  Jan  83 


S.  PERFORMING  010.  REPORT  NUMBER 

N/A  _ 


S.  contract  or  grant  num0CR'’ii 

F19628-80-C-0185 


10.  PROGRAM  element,  project.  TASK 
AREA  S  WORK  UNIT  NUMBERS 

61102F 

2306J236 


!Z.  REPORT  DATE 

February  1984 

'3.  NUMBER  OF  pages 

234 


14.  MOHITORING  AGSNCv  name  4  AOCPESS'i/  diitormnt  /rom  Canfro/fintf  0?/»ce'  j  'S.  SECURITY  CLASS  thi%  rvporr> 

o  _  :  UNCLASSIFIED 


CECLASSjP'CaTICN  downgrading 
.  schedule 
N/A 


*6.  Distribution  statemen*^  ot  Report; 


Approved  for  public  release;  distribution  unlimited 


1?.  Distribution  statement  i  in  Block  JO,  H  dirtorcnt  from  Report) 


IB.  SUPPLCMENTARY  notes 

RADC  Project  Engineer:  Charles  L.  Woods  (ESO) 


19.  KEY  WORDS  (Contlmto  on  rovocoo  otdo  1/  nocoooorr  end  Idonttty  bf  ttoek  ntm»bor) 


Phase  Conjugation 
Optical  Data  Processing 


Degenerate  Four-Wave  Mixing 
Photorefractive  Materials 


20.  abstract  (Continuo  on  roooroo  oido  It  nmcmooorr  end  IdontUr  kr  block  ntmtbor) 

The  purpose  of  this  program  has  been  the  study,  both  theoretically  and 
experimentally,  of  the  use  of  phase  conjugate  optics  for  information 
processing  and  image  transmission  applications.  A  significant  portion 
of  effort  was  directed  toward  an  improved  understanding  of  the  photo- 
refractive  effort  and  its  application  to  two-wave  and  four-wave  mixing. 
We  have  developed  a  model  for  grating  formation  in  photorefractive 
materials  and  have  derived  expressions  for  figures  of  merit,  which  aid 


mm 


UNCLASSIFIED 


—  .  /'  - -  .  - 

-in  the  selection  of  materials  and  the  optimization  of  efficiency  and 
speed.  -We  l»ave  also  developed  a  model  for  two-wave  mixing  in  the 
transient  regime,  and^ljave  measured  response  times  for  both  two-wave 
and  four-wave  mixing.  Other  devices  and  materials  considered  for  ODP 
applications  include  the  Hughes  liquid  crystal  light  valve  and  long-chain 
organic  molecules. 

Measurements  on  several  ODP  configurations  were  performed.  First,  we 
luuia  used  optical  phase  conjugation  to  transmit  an  image  through  a 
1.75  m  multimode  fiber,  with  a  resolution  of  better  than  15  llnes/mm. 
Secondly,  we  have  demonstrated  several  static  convolution  and  correlation 
operations,  using  character  masks.  Finally,  we  Wve  investigated  the 
use  of  phase  conjugation  for  real-time  optical  correlation.  Results 
of  each  of  these  experiments  will  be  described. 

■j  > 


Accession  F 

NTTS  G”\&I 
DTIG  TAB 
Unannounced 


TABLE  OF  CONTENTS 


SECTION  i>aCE 

1  INTRODUCTION .  7 

2  MECHANISMS  AND  MATERIALS  FOR  PHASE 

CONJUGATION .  11 

A.  Photoref ractive  Materials .  13 

B.  Liquid  Crystal  Light  Valve .  35 

C.  Organic  Materials .  39 

3  IMAGE  TRANSMISSION  THROUGH  FIBERS  USING 

PHASE  CONJUGATION .  49 

A.  Transmission  of  Images  Through  Fibers 

Using  Phase  Conjugation .  49 

B.  Phase  Conjugation  Using  Optical  Fibers 

as  the  Non-Linear  Medium .  57 

4  TWO-WAVE  MIXING .  61 

5  OPTICAL  CONVOLUTION  AND  CORRELATION .  7  1 

A.  Static  Two  Dimensional  Optical 

Convolution  and  Correlation .  71 

B.  Real  Time  Correlation .  76 

6  CONCLUSIONS  AND  RECOMMENDATIONS .  91 

REFERENCES .  95 

APPENDICES 

A  PRELIMINARY  REPORT  ON  THE  MATHEMATICS  OF 

PHOTOREFRACTIVITY .  A- 1 

B  SELECTION  OF  PHOTOREFRACTIVE  MATERIALS  FOR 

DEGENERATE  FOUR  WAVE  MIXING .  B-1 

C  OPTIMAL  PROPERTIES  OF  PHOTOREFRACTIVE 

MATERIALS  FOR  OPTICAL  DATA  PROCESSING .  C-1 

D  TWO-WAVE  MIXING  IN  THE  TRANSIENT  REGIME .  J- 1 

E  DEMONSTRATION  OF  IMAGE  TRANSMISSION 

THROUGH  FIBERS  BY  OPTICAI.  PHASE 

CONJUGATION .  E- 1 


3 


LIST  OF  ILLUSTRATIONS 


FIGURE  PAGE 

2-1  Typical  geometry  for  hologram  writing  or  DFWM 

in  a  nonlinear  medium .  15 

2-2  Geometry  for  hologram  writing  or  UFWM  in  a 

uniaxial  crystal  such  as  BaTi03 .  16 

2-3  Experimental  geometry.  Ef,  and  Ep 

are  the  forward  pump,  backward  pump,  and 

probe  waves,  respectivity .  24 

2-4  Charge  density  for  an  incident  probe  pulse 

width  of  30  25 

2-5  Charge  density  for  an  incident  pulse  width  of 

9  ^di .  26 

2-6  Charge  density  for  an  incident  probe  pulse 

width  of  1.5  27 

2-7  Geometry  for  beam  coupling  experiments 

in  BaTi03  at  4416A . 29 

2-8  Weak  beam  (probe)  gain  versus  grating  spacing 

for  3  BaTi0  3  samples .  30 

2-9  Schematic  of  LCLV .  36 

2-10  Configurations  for  microwatt  pliase  con j uga ion 

using  liquid  crystal  valve .  37 

2-11  Measured  MTF  of  liquid  crystal  light  valve  .  40 

2-12  Long  conduction  paths  have  large  dipole  momenta.  42 

2-13  Advantage  of  bonds .  43 

2- 14  Free  electron  model:  One-dimensions  box  with 

reflecting  walls .  44 

3- 1  Demonstration  of  image  transmission  through  a 

fiber .  50 

3-2  Schematic  of  experiment  used  for  image 

transmission  through  fibers .  52 

3-3  Photograph  of  experimental  apparatus  used  i or 

image  transmission  through  fibers .  53 


PREVIOUS  PAGE 
IS  BLANK 


FIGURE 

3-4  Reconstructed  image  after  double  passing 

the  fiber . 

3- 5  Schematic  of  degenerate  four  wave  mixing 

experiments  using  silica  fiber  as  nonlinear 
medium . 

4- 1  Two-beam  coupling  experiment  schematic. . 

4-2  Single  pulse  data  for  two-beam  coupling 

experiment  incident  probe  pulse  lower  trace, 
amplified  probe  pulse  upper  trace . 

4-3  Repetitive  waveform  data  for  two-beam  coupling 

experiment . 

4- 4  Degenerate  four  wave  mixing  signal  when  probe  is 

(i)  A  single  pulse,  (ii)  a  repetitive  waveform, 
(iii)  constant  value . . 

5- 1  Geometry  used  in  two-dimension  correlation 

experiment . 

5-2  Character  mask  used  in  2-D  correlation 

experiment . 

5-3  Normalized  cross-correlation  matrix . 

5-4  Two-dimensional  correlation  experiment 

schematic . 

5-5  Relative  indensity  and  frequency  of  diffracted 

orders  produced  by  acousto-optic  devices . 

5-6  Optical  multiplication  experiment  schematic . 

5-7  Electronics  used  to  drive  acousto-optic  devices 

block  diagram . 

5-8  Fringe  pattern  created  by  pump  beam  and  probe 

beam  at  position  of  conjugator . 

5-9  Data  for  optical  multiplication  of  the  function, 
'^i(t)  and  F2(t).  A.  Pi{t)  =  0.01  Hz  sawtooth, 
Fp(t)  =  0.9  Hz  square  wave.  B.  F^lt)  =  0.5  Hz 
square  wve ,  F2(t)  =  0.05  Hz  sqare  wave . 

5-10  .'.tability  of  various  mounting  techniques . 

5-11  Real  time  correlation  experiment  schematic . 


5-12 


Real  time  auto  correlation  data 


SKCTION  1 


INTRODUCTION 


The  purpose  of  this  program  has  been  the  study,  both 
theoretically  and  experimentally,  of  the  use  of  phase  conjugate 
optics  for  information  processing  and  image  transmission 
applications.  The  objectives  are  to  analyze  and  develop  phase 
conjugation  techniques  for  optical  signal  processing,  to 
determine  the  usefulness  of  phase  conjugation  for  restoring 
spatial  information  transmitted  over  multimode  optical  fibers, 
and  to  study  real-time  optical  convolution  and  correlation  of  two 
spatially  or  temporally  modulated  input  waves.  The  contract  was 
modified  after  its  inception  to  include  a  task  covering  the 
theoretical  and  experimental  study  of  materials  for  use  in 
performing  phase  conjugation  of  spatially  modulated  optical 
waves. 

The  key  theoretical  and  experimental  results  obtained  under 
this  program  have  been  to; 

•  Develop  a  model  of  grating  formation  in  photo- 
refractive  materials  and  derive  expressions  for  figures 
of  merit,  which  aid  in  the  selection  of  materials  and 
tue  optimization  of  efficiency  and  speed. 

•  Measure  efficiency  and  response  time  for  four-wave 
mixing  and  two-wave  mixing  in  photoref ract ive 
materials. 

•  Demonstrate  real-time  holographic  capabilities  of  the 
Hughes  liquid  crystal  light  valve. 

•  Analyze  nonlinear  optical  behavior  of  long-chain 
organic  molecules  and  estimat'*  limiting  values  of  t'ue 
third  order  nonlinear  susceptibility  in  t!ie  visible  and 
microwave  spectral  regions. 

•  Demonstrate  image  transmission  through  a  multimode 
fiber  using  optical  phase  conjugation. 

•  Perform  two-dimensional  optical  correlation  using  phase 
conjugation . 

®  Investigate  the  application  of  phase  conjugation  to 
real-time  optical  information  processing. 
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During  the  course  of  this  program,  we  have  studied  three 
separate  materials  or  devices  for  phase  conjugation  via 
degenerate  four-wave  mixing  (DFWM).  Photoref ractive  materials 
have  the  advantage  of  very  large  efficiences  at  cw  power  levels, 
and  have  moderate  response  times  sec).  We  have 

developed  a  model  for  grating  formation  in  photoref ractive 
materials  and  have  derived  expressions  for  figures  of  merit  for 
comparison  and  optimization  of  materials.  At  present,  BaTi03  or 
3BN  are  preferred  when  high  efficiency  is  required,  while  BSO  is 
preferred  when  speed  of  response  is  most  important.  Our  analysis 
suggest  that  significant  improvements  in  the  response  times  of 
BaTiO^  and  SBN  can  be  obtained  by  improved  control  of  the 
quantity  and  valence  state  of  the  dopants. 

The  second  device  of  interest  is  the  Hughes  liquid  crystal 
light  valve,  a  high-sensitivity,  high-gain  device  for  large 
screen  projection.  We  have  demonstrated  that  the  light  valve  can 
be  considered  as  a  real-time,  thin  holographic  medium,  and  is 
capable  of  correcting  for  aberration  in  optical  beams.  We  have 
sliown  how  tills  device  can  be  adapted  for  conventional  four-wave 
mixing . 

Pinal Ly,  we  have  considered  the  possible  use  of  long-chain 
organic  molecules  for  DFWM.  We  have  developed  a  simple  model  for 
calculating  third-order  nonlinear  susceptibilities  in  these 
materials.  The  important  result  of  this  model  is  that  the 
susceptibility  scales  as  the  third  power  of  the  length  of  this 
molecule,  suggesting  that  giant  values  of  susceptibility  may  be 
obtainable  in  sufficiently  long  molecules. 

The  transmission  of  spatial  information  through  a  multimode 
fiber  is  degraded  for  several  reasons,  two  of  which  are  modal 
dispersion  and  mode  mixing.  The  process  of  phase  conjugation  — 
that  is,  passing  a  spatially  modulated  wave  through  a  fiber, 
conjugating,  and  then  retransmitting  the  beam  back  through  the 
same  fiber  —  can  compensate  for  these  effects  and  produce  a 
restored  image.  We  have  used  optical  phase  conjugation  to 
transmit  an  image  tiirough  a  1.75  meter  multimode  optical  fiber 
with  a  resolution  of  better  than  15  iines/mm. 


The  ability  of  one  nonlinear  technique,  (DFWM),  to  perform 
the  multiplication  of  two  (or  three)  input  signals  in  a  nonlinear 
medium  provides  the  capability  of  generating  the  optical  con¬ 
volution  or  correlation  of  spatially  modulated  waves.  The  use  of 
a  nonlinear  optical  process  to  generate  the  optical  correlation/ 
convolution  is  esentially  a  generalized  Van  der  Lugt  Holographic 
filter  technique  in  which  a  nonlinear  medium  operating  as  a  four- 
wave  mixer  replaces  a  conventional  holographic  film  recording. 

The  key  advantage  of  such  a  system  is  that  the  "hologram"  may  be 
written  and  played  back  almost  instantaneously  (i.e.,  real-time 
signal  processing  can  be  performed).  Because  of  this  potential, 
we  have  investigated  experimentally  both  spatial  and  temporal 
correlation  geometries. 

We  have  performed  two  dimensional  correlation  between  a 
single  character  and  matrix  of  characters  using  phase 
conjugation.  The  phase  conjugator  is  placed  at  the  common  back 
focal  plane  of  three  lenses.  The  fields  at  this  position  are 
proportional  to  the  Fourier  transform  of  the  front  focal  plane 
fields.  A  mask  of  characters  is  placed  at  the  front  focal  plane 
of  two  of  the  lenses  and  the  focal  point  of  a  fourth  lens  is 
positioned  at  the  remaining  front  focal  plane.  The  phase 
conjugate  signal  generated  is  proportional  to  the  correlation  of 
the  illuminated  characters. 

In  contrast  to  the  static  spatial  masks  used  in  the  above 
experiments  we  also  investigated  the  use  of  acousto-optic  Bragg 
cells  for  imposing  temporal  information  to  be  correlated.  in 
this  case,  the  signals,  f^(t)  and  f^(t),  to  be  correlated  are 
encoded  on  the  optical  beam  by  two  acousto-optic  (A.O. ) 
modulators.  The  first  A.O.  modulator  is  used  to  encode  f|{t)  as 
a  temporal  intensity  variation.  The  second  A.O.  modulator  is 
used  as  a  delay  line  which  produces  a  spatial  variation 
proportional  to  f ^ ( t ) • f ^ ( t-i ) .  The  modulated  probe  then 
propagates  to  the  phase  conjugator.  The  phase  conjugate  signal 
produced  is  proportional  to  the  correlation  of  tj(t)  and  f  (t). 

In  this  manner  we  investigated  the  application  of  p'uase 
conjugation  to  real-time  optical  correlation. 
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MeCHANISi^S  AND  MATKKIAL.S  (•'OK  P:1ASI';  CON J UGA'I’ i  Oi J 

A  large  number  of  nonlinear  mechanisms  have  been  eropo -.e  I  or 
demonstrated  which  have  promise  for  phase  conjujation  oi 
spatially-modulated  optical  signals  via  degenerate  Onjr-wave 
mixing  (DFWM).  We  have  performed  a  preli;ninar/  sarvr/  of 
mechanisms,  along  with  the  perforinance  of  specific  a.it.  e  r  Lai  s ;  rate 
results  are  given  in  Table  1.  Several  of  the  f  ibnlaiei  Knt  ‘r  i  a  1  r, 
have  the  desirable  properties  of  reipiLrinq  molest  (  -w}  pump  po^er 
levels  and  having  response  times  in  the  '  1 1 is  M.’ond  or  sub¬ 
millisecond  region.  It  is  expected  th  ■  most  ol:  thi-  mate'-ials 
should  operate  at  spatial  frequenices  r  .i.oaching  10J0  liner,  .nm 
(1  pm  resolution  in  the  material),  and  "s  be  saLtai'lr-  :  o’- 
handling  high  spatial  bandwidth  signals 

Among  the  mechanisms  and  materials  listed  in  J'a  O  e  1,  v.’e 
have  chosen  three  for  further  study:  photore  f  r act  L  re  ,nai.i> r  l a  1  s  , 
the  liquid  crystal  light  valve  an<3  long-chain  organic  oolj’.ier":. 

Photoref ract ive  materials  operate  over  a  ri.*latively  wide 
wavelength  range  at  modest  cw  power  levels,  and  have  r'-sponra} 
times  in  the  sec  range.  Photoref  ract  i  re  na'.erijls  with 

high  promise  include  BaTi03  and  SliW  (wit:)  demorist rated  liiga  DFv'.’M 
reflectivity),  KNb03  (moderate  reflectivity  with  fast  res;)ons:), 
and  3S0  (lower  reflectivity  with  very  fast  response).  .L)ni 
characteristics  of  BaTi03,  KNb03  and  BSO  are  included  ii  Table  1. 

Although  BSO  is  an  electro-optic  crystal^^  with  optical 
16 

activity  ,  it  has  been  used  to  demonstrate  high  sensitivity 
holographic  storage  ,  joint  transform  processing  ,  and 
incoherent  optical  conversion  , 

The  liquid  crystal  light  valve  is  a  state-of-Lhe-ar t  i^'vLce 
developed  at  ilughes  Research  Laboratories  for  lar  icr-'Mi 
projection  at  TV  scan  rates.  It  is  b.isically  an  cp.ical-to- 
optical  image  transducer  that  ii.as  high  input  L  t  i  v  i  '  y 

(~50  MW/cm^)  and  high  gain,  allowing  prf)jecf.  Lon  >.’1!:  i  i  it  vci  ' 
light  sources.  Although  the  LCLV  was  dev'-lopi’d  f  rr.,)!-'  i  ,i;  ' 
imaging  applications,  it  is  re,j<li  ly  a!Lapta!)lr‘  for  r'eil-tb'K' 
holography  and  the  optical  data  orrici.'s:.  Ln  ;  Minot  L;iri,:  d  oive  i  ‘  ro: 
its  holographic  properties. 
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Organic  molecules  in  the  form  of  long-chain  polymers  witli 
conjugated  bonds  are  very  promising  nonlinear  materials  because 
of  the  large  induced  dipole  which  results  from  the  long  length  ol: 
the  molecule.  Experiments  performed  to  date  show  the  expected 
scaling  of  the  nonlinear  coefficient  with  chain  length.  The 
analysis  described  below  allows  the  estimation  of  the  required 
length  to  yield  values  of  on  the  order  of  10”^  esu. 

Such  a  nonlinearity  would  allow  high  efficiency  DH’WM  with  low- 
power  cw  lasers. 

In  this  section  we  describe  the  results  of  analysis  and 
experiments  on  photoref ractive  materials,  allowing  the  design  and 
selection  of  materials  for  a  given  application.  We  next  describe 
measurements  intended  to  characterize  the  real-time  holograp'nic 
properties  of  the  liquid  crystal  light  valve.  Finally,  we 
describe  our  analysis  of  the  nonlinear  optical  behavior  of  long- 
chain  polymers,  leading  to  the  estimation  of  the  requirerneiits  for 
third-order  susceptibilities  on  the  order  of  esu. 

A.  PHOTOREFRACTIVE  MATERIALS 

1 .  Analysis 

Photoref ractive  materials  are  the  tnost  promising  candidate 
materials  for  optical  data  processing  at  the  present  time.  In 
this  program  we  have  carefully  analyzed  and  expanded  upon 
existing  models  for  grating  formation  in  photorefractive 
materials,  and  we  have  established  figures  of  merit  for  selecting 
an  optimum  material  for  a  given  experiment. 

The  major  part  of  our  analysis  of  photorefractive  materials 
is  summarized  in  Appendices  A,  B  and  C.  Appendix  A  is  a  report 
by  Marshall  Sparks  which  details  the  physical  assumptions  used  in 
this  and  subsequent  analytical  work  and  derives  expressions  Cor 
the  steady  state  space  charge  electric  field  in  a  number  of 
operating  regimes.  The  time  response  of  the  space  charge  field 
was  also  calculated  for  certain  limiting  assumptions. 


13 


Appendix  B  is  a  manuscript  which  was  prepared  for  the  SPIli 
San  Diego  Meeting  in  August,  1982.  It  reviews  the  analysis  of 
grating  formatior,  and  presents  exact  expressions  for  the  time 
constants  in  the  limit  of  small  modulation  depth.  Several 
performance  parameters  were  also  discussed,  and  the  application 
to  materials  selection  was  considered. 

Appendix  C  is  a  preprint  of  a  draft  manuscript  to  be 
submitted  to  a  special  issue  of  Optical  Engineering  on  Optical 
Data  Processing.  In  this  manuscript  the  analysis  of  grating 
formation  is  again  reviewed.  Various  material  figures  of  merit 
are  then  considered  in  detail,  and  these  parameters  are  evaluated 
for  BaTi03  and  BSO  in  order  to  compare  their  performance  in  a 
number  of  applications. 

The  simple  holographic  picture  of  grating  formation  and  DFWM 
is  outlined  in  Figure  2-1.  In  this  scheme  a  holographic  grating 
is  written  by  the  beams,  and  E2.  A  third  beam,  E3,  scatters 
from  this  grating  and  creates  the  signal  beam,  Ei^,  which 
propagates  in  the  backward  E2  direction.  The  nonlinear 
reflectivity  is  defined  as  R  =  \Ei^/E2  |^. 

It  is  also  important  to  keep  in  mind  that  the  nonlinear 
crystals  used  for  photoref ractive  applications  are  highly 
anisotropic.  In  particular,  the  effective  electro-optical 
coefficient  depends  strongly  on  the  polarization  of  the  beams  and 
on  the  angle,  0.  between  the  c-axis  and  the  grating  wave  vector 
{see  Figure  2-2).  For  extraordinary  polarized  beams  in  BaTi03, 
tne  effective  electro-optic  coefficient  for  9=0  is  r33.  However, 
for  o>o  there  is  an  added  contribution  from  r^2*  Since  r^2 
much  larger  than  r33  (see  Table  3),  the  grating  efficiency  or 
DFVJM  reflectivity  increases  rapidly  with  0  up  to  value  of  ~40°. 
Because  of  the  large  refractive  indices  of  BaTi03,  this  optimum 
angle  cannot  be  obtained  unless  the  crystal  is  immersed  in  index 
matching  fluid.  In  all  our  experiments  with  BaTiO^  the  crystal 
is  in  air  and  the  maximum  value  of  9  is  ~20°.  Even  at  this  value 
of  0,  two-wave  mixing  gains  on  the  order  of  20,  and  four-wave 
mixing  reflectivities  on  the  order  of  200-400%,  can  easily  be 
obse  rved . 


11368  6H1 


Figure  2-1 


Typical  geometry  for  hologram  writing  or  DFWM 
nonlinear  medium.  In  holographic  language,  '*1 
the  reference  writing  beam,  Ep  is  the  image  writing 
beam,  5  3  is  the  radout  beam,  and  is  the  signal 
beam.  In  the  language  of  DPWM,  h:  1  and  l]  <  are  the 
pump  beams,  R2  is  the  probe  beam,  and  R ,,  is  the 

the  grating 


signal  beam, 
spacing . 


The  parameter 
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leometr/  for  hologram  writing  or  DFWM  in  a  uniaxial 
crystal  such  as  Bali03.  The  beams  have  the  same 
identification  as  in  Figure  1.  The  parameter  kg 
is  the  grating  wave  vector,  given  by 


=  2’'/A, 


g- 


1  6 


Two  separate  models  for  describing  charge  transport  in 
photoref ractive  materials  have  been  formulated.  The  hopping 
model ^  assumes  that  carrier  transport  occurs  via  hopping  from  a 
filled  donor  site  to  a  neighboring  empty  trap.  This  model  was 
originally  developed  to  describe  electrical  conduction  in  semi- 
insulating  materials  and  amorphous  semiconductors.  The  band  trans¬ 
port  model^'^'^^  assumes  that  electrons  (or  holes)  are  optically 
excited  from  filled  donor  (or  acceptor)  sites  to  the  conduction 
(or  valence)  band  where  they  migrate  to  dark  regions  in  the 
crystal  by  drift  or  diffusion  before  recombining  into  an  empty 
trap.  In  both  models  the  transported  charges  result  in  an  ionic 
space  charge  grating  which  is,  in  general,  out  of  phase  with  the 
incident  irradiance.  The  space  charge  grating  is  balanced  by  a 
periodic  space  charge  electric  field  in  accordance  with  Poisson's 
equation.  This  space  charge  field  modulates  the  refractive  index 
through  the  electro-optic  effect.  The  band  transport  model  has 
been  developed  to  a  greater  degree  and  is  in  wider  use  in  the 
literature;  this  is  the  model  which  we  have  used  in  our 
analysis. 

In  order  to  analyze  the  DFWM  properties  of  ferroelectr ics , 
we  first  require  expressions  for  the  amplitude  and  spatial  phase 
of  the  space  charge  electric  field.  In  the  Appendices  the 
required  expressions  have  been  derived  (using  the  band  transport 
model)  for  specific  limiting  cases.  One  case  of  relevance  for 
many  experiments  is  m<<1,  where  m  is  the  fractional  modulation  of 
the  input  irradiance.  For  this  case^'^'^^  the  solution  for  the 
steady  state  space  charge  field  amplitude  is 


=  m 

sc  q 


1/2 


(  1  ) 
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where 


=  Applied  drift  field  (normal  to  grating  planes), 

=  —  k  (diffusion  field), 
eg' 

4TTeN- 

A 

=  ——r —  (limiting  spacecharge  field). 


In  the  above  expressions,  k  is  the  Boltzman  constant,  kg  is 
the  grating  wave  number,  is  the  density  of  empty  traps,  and 
e  is  the  dielectric  constant.  The  field  Eg  is  that  which 
results  from  the  separation  of  all  available  charges  by  one 
grating  period.  The  spatial  phase,  ip  of  Egj,  (relative  to 
the  irradiance),  is  given  by 


tan  '!>  = 


The  general,  Eo<<Eq,  so  that  with  no  applied  drift  field 
^sc  ”  *^0'  'I'  “  ^/2.  This  is  the  diffusion  limit, 

for  which  a  characteristic  spatial  phase  shift  between  the 
i.rradiance  and  the  space  charge  field  is  observed.  As  the 
applied  drift  field  is  increased  from  zero,  we  have  Eg^.  “  m 
(Eq^  +  Eq^)^^^,  and  the  spatial  phase  is  intermediate  between  0 
and  Tf/2.  Finally,  for  large  values  of  drift  field  (Eo>Eq),  the 
space  charge  field  saturates  at  Eg^,  “  m  Eq,  and  the  phase 
shift,  P,  is  once  again  if/2. 

Note  that  for  small  values  of  the  drift  field  (EQ<Eq),  the 
space  charge  field  is  independent  of  all  material  parameters.  In 
this  case  the  induced  refractive  index  variation  is  determined 
entirely  by  the  relevent  electro-optic  coefficient.  Thus,  for  large 
values  of  index  modulation,  materials  with  large  values  of  electro¬ 
optic  coefficient  are  desirable.  However,  materials  in  this 


category  (e.g.,  BaTi0  3,  S3N)  tend  to  be  slow.  Specifically,  t'ne 

band  transport  solutions  for  small  modulation  fraction  give  the 
write  or  erase  time^  as 


where  is  the  dielectric  relaxation  time, 

_  e  _  _ _ 

di  4  no  4iTneP  ' 


(3) 


(4) 


where  a  is  the  conductivity,  n  is  the  free  carrier  density,  and  tJ 
the  mobility.  The  function,  f(kg,EQ),  accounts  for  t!ie  spatial 
variation  of  the  free  carrier  density  (see  Apoeniiices  B  and  C).  T 
many  cases  of  physical  interest,  f(kg,Eo)  “  1  and  = 

■^di- 

If  the  conductivity,  a,  is  dominated  by  the  contribution  from 
photocarriers,  then  a  simple  rate  equation  solution  gives 


n 


(5) 


where  s  is  the  photoionization  cross  section,  is  tlie 
recombination  rate  coefficient  ,  is  the  density  of  filled 

traps,  (as  defined  earlier)  is  the  density  of  empty  traps,  and 
Iq  is  the  average  irradiance.  By  combining  Equations  (4)  and 
(5),  we  obtain 
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di 


e  Y. 


4Tf  eiis 


J_  1 

I  R 
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(6) 


where  the  reduction  ratio  R  is  defined  as 


The  quantity  in  parentheses  in  P^quation  (6)  can  be  considered  as  a 
material  figure  of  merit,  having  units  of  energy/unit  area.  Once  a 
material  is  chosen,  only  limited  experimental  control  over  this 
quantity  is  possible  (see  Appendix  C).  However,  the  quantity  K  can 
be  varied  over  many  orders  of  magnitude  by  reduction  or  oxidation  of 
a  given  sample.  It  is  clear  from  Kquation  (6)  that  increasing  R 
through  chemical  reduction  can  significantly  reduce 

Values  of  R  as  high  as  10  are  not  hard  to  achieve  in  many  cases,  and 
in  some  materials  (e.g.,  8SO)  values  of  R  on  the  order  of  10^  are 
typical.**  Note,  however,  that  large  increases  in  R,  if  accompanied 
by  a  reduction  in  W^,  can  reduce  the  limiting  grating  efficiency 
by  reducing  Eq  (see  Appendix  C).  Note  also  the  dependence 

in  Equation  (6),  indicating  the  direct  relation  between  speed  and 
irradiance . 

In  this  review  of  grating  formation  in  photoref ractive 
materials,  it  was  assumed  that  the  incident  fields  in  each  volume 
element  are  known.  In  an  actual  beam  coupling  or  DFWM  experiment, 
only  the  incident  <)ptical  fields  on  the  sample  boundaries  are  known. 
The  complete  grating  solution  requires  the  suostitution  of  the 
refractive  index  change  (as  a  function  of  the  local  optical  fields) 
into  the  wave  equations  which  describe  the  variation  of  the  optical 
fields  with  position  in  the  material.  This  results  in  four  coupled- 
wave  equations  for  Dfvm,  and  two  equations  for  the  simpler  case  of 
beam  coupling  (or  two-wave  mixing).  In  these  equations,  the 
contribution  from  both  the  small-period  and  large-period  gratings 
must  be  considered,  as  well  as  self-interaction  effects.  Several 
approaches  to  the  solution  of  the  coupled-wave  equations  have  been 
presented .  ^  In  the  liinit  of  low  pump  depletion  and  small  beam 

coupling,  the  coupled  wave  approach  is  not  necessary,  and  the  volume 
grating  solutions  are  adequate  for  analysis.  Even  when  large  beam 
oupling  is  present,  the  volume  grating  solutions  for  the  write  and 
erase  times  are  still  correct. 

On  the  basis  of  the  above  analysis,  we  may  define  four 
P'jr  i  (srmanco  param<jters  (or  figures  of  merit)  wliich  characterize  a 
given  material  for  a  particular  application.  These  parameters 
(described  in  detail  in  Appendix  C)  are:  steady  state  index  change 


(Angs)/  response  time  C^w^'^e)'  write  energy  (w),  and 
photorefractive  sensitivity  {S=dn/dw).  The  relevance  oE  each 
parameter  will  vary  with  the  application.  For  example,  when  high 
diffraction  efficiency  or  four-wave  mixing  reElectivity  is  required, 
the  steady-state  index  change  will  be  the  most  important  parameter. 
If  high  speed  operation  is  desirable,  then  the  response  time  will  be 
most  important.  If  limited  optical  energy  is  available,  then  the 
write  energy  will  be  most  relevant.  ^nd  finally,  if  sensitivity  to 
small  changes  is  critical,  then  the  sensitivity  will  be  most 
relevant . 

In  Table  2  we  have  listed  performance  parameters  for  BaTiOi  and 
BSO  for  two  values  of  grating  spacing  and  applied  drift  field.  The 
relative  advantages  of  these  two  materials  is  discussed  in  Appen¬ 
dix  C.  Essentially,  the  large  electro-optic  coefficient  of  BaTiOj 
makes  this  material  preferred  when  high  efficiency  is  required.  On 
the  other  hand,  BSO  is  preferred  where  fast  response  time  or  low 
writing  energy  is  required.  Note  that  neither  material  has  been 
optimized  for  photorefractive  applications.  In  particular,  it  is 
expected  that  control  of  the  dopant  species,  valence  state  and 
concentration  in  BaTi03  should  reduce  the  response  time  by  a 
significant  amount. 

The  volume  grating  solutions  summarized  above  provide  the 
steady  state  and  step  function  response  for  the  charge  density  and 
space  charge  field,  but  are  not  directly  applicable  for  pulsed 
excitation.  In  Appendix  D  we  present  a  model  for  grating  formation 
with  pulsed  excitation.  This  model  (summarized  below)  also  applies 
to  two-wave  mixing  and  four-wave  mixing  for  the  case  of  nondeplete<3 , 
strong  pumps  and  a  weak  probe  wave.  This  work  was  supported  by 
Hughes  IR&D  funds  and  the  purpose  of  including  this  model  is  to 
allow  comparision  with  the  pulsed-probe  two-wave  mixing  studies 
described  in  Section  5. 

Our  model  relies  on  the  following  asumptions: 

•  The  pumps  are  cw. 

•  The  incident  probe  is  a  pulse  of  duration,  fp,  and 
magnitude,  •^’q. 


•  The  interaction  geometry  is  as  shown  in  Figure  2-3. 

For  DFWM,  all  three  input  waves  are  present.  For  two- 
wave  missing,  we  shall  set 

•  The  pump  waves  are  nondepleted  during  the  nonlinear 
interaction. 

Using  the  charge  transport  model  of  Kukhtarev^'^  and 
Laplace  transform  techniques,  we  have  obtained  solutions  for  tbis 
charge  density  and  the  space  charge  electric  field.  We  then 
calculated  the  two-wave  response  of  BaTi0  3  to  a  sijuare-wave  probe 
pulse,  with  three  possible  time  durations.  The  assumed  conditions 
are : 

Dielectric  Relaxation  Rate:  3.0  sec"' 

Diffusion  Rate:  I.OxlO'^  sec"' 

Photoionization  Rate:  1.0x10*^  sec"' 

Electron  Recombination  Rate:  1.0x10“^  sec"' 

Ion  Recombination  Rate:  1.0x10^-  sec"' 

In  Figure  2-4  we  have  plotted  the  charge  density  versus  time  for 
a  pulse  width  of  30  Because  of  the  long  pulse  width,  the 

material  is  in  the  steady  state  during  most  of  the  pulse  duration. 
In  Figure  2-5  the  pulse  width  is  only  10  and  the  material 

is  only  briefly  in  the  steady  state.  For  a  sufficiently  short  pul 
(Figure  2-6)  the  system  does  not  reach  steady  state  and  the  peak 
signal  is  reduced.  The  results  of  section  5  clearly  follow  the 
behavior  plotted  in  Figures  2-4-2-6. 
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Rxper  iinents 

In  the  experimental  phase  of  our  materials  study,  we  have 
measured  the  photorefractive  properties  of  a  number  of  materials  via 
two-wave  mixing  (beam  coupling)  and  four-wave  mixing.  We  have  also 
developed  a  technique  for  poling  the  BaTi03  crystals  required  for 
our  experiments. 

In  our  beam  coupling  experiments,  we  have  measured  the  steady- 
state  weak-beam  gain  at  4416A  as  a  function  of  grating  spacing  for 
several  BaTi03  samples.  The  geometry  for  our  experiment  is  shown  in 
Pigure  2-7.  Both  input  beams  were  derived  from  a  Liconix  He-Cd 
laser;  the  pump  power  was  25  mW  and  the  probe  power  was  3  mW.  The 
beams  were  unfocused  and  the  spot  diameter  at  the  crystal  was  ~2  mm. 
The  crystal  was  oriented  so  that  the  grating  wave  vector,  kg,  was 
aligned  parallel  to  the  crystal  c-axis.  For  p-polarized  input  beams 
(extraordinary  polarization  inside  the  crystal)  the  induced 
refractive  index  change  depends  in  a  complicated  jaanner  on  the 
electro-optic  coefficients,  ri3  and  r33.  Furthermore,  extra¬ 
ordinarily  polarized  beams  are  subject  to  beam  fanning,®  which 
creates  significant  beam  distortion.  We  have  thus  chosen  to  use 
beams  with  ordinary  polarization  in  the  crystal;  the  appropriate 
electro-optic  coefficient  is  thus  r^3.  For  this  orientation  the 
beam  coupling  is  relatively  small  (since  i^i3<i^33)»  so  that  the 
effects  (jf  pump  depletion  are  not  significant. 

In  Figure  2-8  we  plot  the  weak  beam  gain,  G  (output  intensity 
divided  b/  input  intensity),  as  a  function  of  grating  spacing,  Ag, 
for  three  BaTi03  samples  (all  purchased  from  Banders  Associates), 
the  gain  can  be  written  as 

G  =  e"^*"  (3) 

wiiiire  r  is  the  gain  coefficient  and  L  is  the  interaction  length. 
Kukiitarev'^  has  shown  that  the  gain  coefficient,  Y,  is  proportional 
to  the  space  charge  field,  l^s^.  Thus,  from  Equation  (1) 


we  L  i.  nfi 


since  Eq=0. 

By  substituting  the  definitions  for  an;l 
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Thus,  for  small  values  of  A  ,  Y  ~  E  ~ 

1  '3  q 

Y  ~  Ej~  A  The  peak  value  of  Y(or  G) 

Q  g 

for  or  =  tt  e  kT/e^N^.  Thus, 

^  -^o 

allows  the  determination  of  N^. 


A  .  For  large  values  o^.  A 

q  q 

as  a  function  of  A  occurs 

q 

measurement  of  A 

■^o 


The  data  of  Figure  2-8  are  consistent  with  the  above  model. 

For  large  Ag,  the  gain  for  all  crystals  approaches  the  same  value 
(Ej  is  independent  of  N;^).  For  small  Ag,  the  curves  rise  with 
different  slopes  and  peak  at  different  values  of  Ag ,  consistent 
with  different  values  of  N;^.  For  crystal  GB3  we  find 

=  2x10^^/cin^,  while  N;[^  =  4x10^^/cm^  in  crystal  Gb.  The  value 
for  GB3  is  close  to  that  measured  in  a  separate  sample^  by  a  similar 
technique.  This  consistency  lends  support  to  the  numerical  value 
=  2x10^^/cm^,  which  was  chosen  for  the  figure  of  merit 
calculations  in  Appendix  C. 

In  a  separate  set  of  experiments,  we  have  measured  the 
reflectivity  and  response  time  for  DFWM  in  four  different  crystals 
at  6471  A.  In  our  experiments  the  two  counterpropagat ing  pump  oeams 
both  had  incident  intensities  of  ~10  mW/cm^,  while  the  probe 
intensity  was  ~1  mW/cm^.  The  angle  between  the  probe  and  forward 
pump  was  40°,  corresponding  to  a  grating  spacing  of  1.1  pm.  I'iu' 
results  of  our  experiments  are  summarized  in  Table  3.  For  each 
material  the  laser  polarization  and  crystal  vor  ient  it  ion  were 
adjusted  to  exploit  a  particular  component  of  tfie  electr(i-opt  i  c 


coefficient.  Note  that  the  DFWM  reflectivity  varies  inoriotDn  i  .'a  I  1  y 
with  the  value  of  the  electro-optic  coefficient,  as  expected  i:  t'le 
charge  transport  is  dominated  by  diffusion.  The  material  with  the 
shortest  response  time  is  BGO,  which  is  consistent  with  tti>'  large 
values  of  reduction  ratio  expected  in  this  material.  The 
particularly  slow  response  times  in  BNN  and  LiNbOj  suggest  that  tin? 
reduction  ratios  are  much  lower  in  these  materials. 

To  conclude  this  section,  we  will  describe  our  experiments  on 
the  poling  of  BaTi03.  An  as-grown  crystal  of  BaTiOj  consists  <jf 
small  spontaneously  polarized  regions  called  domains.  These  domain;, 
can  be  polarized  in  one  of  6  possible  directions:  along  one  of  tiie 
three  cubic  axes  and  with  positive  or  negative  polarity.  I'he  do;:iain 
generally  arrange  themselves  so  that  the  net  macroscopic 
polarization  is  zero  in  any  direction.  The  possible  number  of 
polarization  directions  can  be  reduced  to  two  by  pressing  the 
crystal  on  2  pairs  of  cubic  faces.  This  is  known  as  mechanical 
poling,  and  it  leaves  the  crystal  with  only  one  axis  of 
polarization,  but  with  both  possible  senses  along  this  axis.  Tlie 
remaining  domains  (called  180°  domains)  must  be  removed 
electrically;  i.e.,  an  applied  electric  field  will  cause  the 
parallel-polarized  doinains  to  grow  and  the  antiparallel-polarized 
domains  to  shrink  and  disappear.  As-purchased  crystals  from  dander.; 
Associates  have  been  mechanically  poled  by  the  supplier;  thus,  the 
required  task  is  electrical  poling. 

There  are  two  well-known  problems  associated  with  the 
electrical  poling  of  BaTiOj:  (1)  the  crystals  are  very  .susceptible 
to  fracture  if  they  are  heated  or  cooled  tlirough  the  Curie 
transition;  (2)  the  crystal  face  in  contact  wit’u  the  positive  p<?ling 
electrode  suffers  severe  localized  fractures,  the  extent  of  which 
are  proportional  to  the  temperature  of  the  crystal  and  the  amount  of 
time  the  poling  voltage  is  applied.  To  circumvent  tliese  (problems, 
we  only  heat  the  crystal  to  "'S"  below  the  Curie  temperature,  and  we 
start  cooling  the  crystal  as  soon  as  the  poling  field  (5-10  kV/'cm) 
is  applied.  We  use  silver  paint  as  electrodes,  and  we  im.ierse  the 
crystal  in  a  heated  beaker  of  high  viscosity  dielectric  oil 


in  order  to  provide  a  stable  temperature  environment.  A  complete 
heating/cooling  cycle  requires  approximately  3  hours. 

After  the  completion  of  a  poling  run,  it  is  important  that  an 
accurate  system  for  measuring  the  degree  of  poling  be  available.  We 
use  two  techniques  for  measuring  the  degree  of  poling.  The 
preferred  technique  is  the  measurement  of  the  electro-optic 
coefficient,  r^=r 3 3- ( no/ng ) ^r i 3 .  This  coefficient  is  most 
easily  measured  by  applying  a  modulating  field  along  the  c-axis  and 
measuring  the  amplitude  modulation  induced  on  a  beam  propagating  at 
90°  to  the  c-axis,  and  polarized  at  45°  to  the  c-axis.  In  an 
unpoled  crystal  the  induced  retardation  from  domains  polarized  along 
+C  approximately  cancels  the  retardation  from  domains  polarized 
along  -C,  and  very  small  amplitude  modulation  is  observed.  When  the 
crystal  is  100%  poled,  the  amplitude  modulation  is  large  and  the 
measured  value  of  r^  is  at  a  maximum.  We  have  poled  two  crystals 
since  the  poling  technique  described  above  was  developed.  In  each 
case  the  measured  value  of  r^,  before  poling  was  approximately  100 
times  smaller  tiian  the  literature  value.  After  electrical  poling, 
the  measured  coefficient  was  a  factor  of  2.4  smaller  than  the 
literature  value  for  both  crystals.  It  is,  of  course,  possible  that 
both  crystals  are  only  ~40%  poled.  However,  the  similarity  in  the 
two  measured  values  suggests  that  the  literature  value  may  be  in 
error . 

A  second  technique  for  measuring  degree  of  poling  is  to  measure 
the  weak  beam  gain  in  a  beam  coupling  experiment.  This  should  be 
lone  for  values  of  >  2  pm,  in  order  to  eliminate  the 
contribution  of  to  the  results  (i.e.,  the  diffusion  limit  is 
preferred).  The  weak-beam  gain  coefficient  is  proportional  to  the 
degree  of  poling.  This  technique  has  the  advantage  of  not  requiring 
an  electroded  crystal;  the  major  disadvantage  is  that  the 
experimental  arrangement  is  more  complicated. 


B. 


LIQUID  CRYSTAL  LIGHT  VALVR 


As  mentioned  earlier,  the  liquid  crystal  light  valve  (LCLV)  is 
an  optical-to-optical  image  transducer  which  has  high  input 
sensitivity  {~50  uW/cm'^)  and  high  gain,  allowing  projection  with 
intense  light  sources.  As  shown  in  Figure  2-9,  the  device  consists 
of  a  sandwich  of  thin-films  that  electrically  control  the  optical 
birefringence  of  a  thin  liquid  crystal  layer.  With  no  writing 
illumination  ("off"  state)  the  applied  voltage  does  not  appear 
across  the  liquid  crystal  layer  because  the  resistivity  of  the 
silicon  photoconductive  layer  is  high.  Polarizers  are  arranged  in 
the  projection  beam  so  that  with  no  writing  illumination  there  is  no 
output  signal  (dark  output).  When  a  spatially  varying  writing  beam 
illuminates  the  silicon  photoconductive  layer,  the  resistivity  o^ 
this  layer  is  reduced  in  the  illuminated  regions,  and  part  of  the 
applied  voltage  appears  across  the  liquid  crystal  in  proportion  ta 
the  strength  of  the  illumination.  This  applied  voltage  alters  the 
birefringences  of  the  liquid  crystal  layer,  and  thus  a  transmitted 
signal  appears  in  projection.  One  important  feature  of  the  LCLV  is 
that  the  dielectric  mirror  and  light  blocking  layer  provide  high 
optical  isolation  between  input  and  output,  thereby  allowing  large 
optical  gain.  The  performance  characteristics  of  the  LCLV  are 
summarized  below: 

•  Resolution:  70  lines/mm  (35  cycles/m;n) 

•  Sensitivity:  20  uW/cm^ 

•  Rise/decay  times:  5  mS  (livquid  crystal  limited) 

•  Wavelength  range:  400-1100  nm 

®  Aperture:  2  inches  diameter. 

The  important  advantages  of  the  LCLV  for  phase  conjugation  are 
the  low  operating  power  levels  and  the  relatively  fast  response 
time.  Two  schemes  for  utilizing  the  light  valve  for  phase 
conjugation  are  shown  in  Figure  2-10.  In  its  standard 
configuration,  the  light  valve  is  analogous  to  a  conventional 
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holographic  medium.  The  two  write  beams  produce  a  grating  pattern 
in  the  photoconductor  which  is  transferred  to  a  thin  phase  grating 
in  the  liquid  crystal,  which  is  read  out  on  the  output  side.  The 
advantage  of  this  scheme  is  the  large  optical  gain;  however,  the 
scheme  is  awkward  for  aberration  compensation.  On  the  other  hand, 
by  re\noval  of  the  blocking  layer  and  dielectric  coating  all  beams 
can  be  brought  in  from  the  output  side  of  the  device,  resulting  in  a 
more  conventional  four-wave  mixing  geometry.  The  only  disadvantage 
is  that  no  optical  gain  is  present;  all  beams  must  have  intensities 
in  the  50  y'W/cm^  range.  The  phase  conjugation  mechanism  in  this 
configuration  is  as  follows:  The  reference  (pump)  and  image  (probe) 
beams  interfere  in  the  photoconductive  layer,  producing  a  spatial 
variation  in  conductivity.  This  spatial  pattern  is  then  transferred 
to  the  liquid  crystal  as  a  spatial  pattern  in  the  refractive  index, 
which  is  then  read  out  by  the  reflected  pump  beam,  producing  a 
signal  (phase  conjugate)  beam  in  the  retro  probe  direction.  The 
major  advantages  of  such  a  device  are  its  high  sensitivity  (20 
uW/cm^)  and  broad  wavelength  coverage.  Because  of  the  thin-film 
nature  of  the  liquid  crystal  layer,  the  DFWM  reflectivity  is  not 
expected  to  exceed  ~20%.  The  resolution  (70  lines/mm)  is  not  a 
severe  limitation  because  of  the  large  aperture  available. 

We  have  recently  performed  several  holographic  grating 
experiments  with  an  "off-the-shelf"  Hughes  LCLV.  In  the  first 
experiment,  tlie  two  write  beams  (derived  from  a  He-Ne  laser)  were 
collimated  beams  wilth  ~5  mm  diameter  incident  at  an  adjustable 
relative  angle.  The  readout  beam  (from  a  separate  He-Ne  laser)  were 
also  collimated  and  had  a  diameter  of  ~2  mm.  The  locations  of  beams 
on  opposite  faces  of  the  LCLV  were  set  for  optimum  spatial  overlap. 
With  no  voltage  applied  to  the  LCLV,  only  a  specular  reflected 
readout  beam  was  observed.  With  the  voltage  set  for  optimum 
operation,  a  holographic  grating  is  produced  in  the  liquid  crystal 
layer,  leading  to  diffraction  orders  in  the  reflected  readout  beam. 
We  have  measured  the  first  order  scattering  efficiency  as  a  function 
of  angular  separation  of  the  write  beams.  This  measures  the 
modulation  transfer  function  (MT^)  of  the  light  valve  since  the 


spatial  frequency  is  uniquely  determined  by  the  angular  separation 
and  the  operating  wavelength.  An  experimental  plot  of  HTF  is  given 
in  Figure  2-11.  The  limiting  spatial  frequency  is  ~25  cycles/mm, 
compared  with  ""SB  cycles/mm  obtained  from  improved  Hughes  light 
valves . 

In  a  second  experiment  with  the  LCLV,  the  principle  of 
aberration  correction  was  demonstrated.  The  writing  beams  consisted 
of  a  collimated  beam  and  a  beam  from  a  lens.  The  readout  beam  was 
also  collimated,  and  the  diffraction  signal  beam  was  shown  to  be 
converging,  with  a  convergence  angle  equal  to  the  divergence  angle 
of  the  image  writing  beam.  If  we  consider  a  lens  as  a  simple  phase 
aberrator,  then  a  second  lens  placed  in  the  signal  beam  (at  the  sam<i 
distance  from  the  LCLV  as  the  lens  in  the  image  writing  beam)  would 
collimate  the  signal  beam.  This  was,  in  fact,  demonstrated,  showing 
that  the  aberration  of  the  lens  was  corrected. 


C.  ORGANIC  MATERIALS 

As  mentioned  earlier,  long-chain  organic  polymers  have  been 
considered  as  possible  candidates  for  OFWM  (a  third  order  nonlinear 
mechanism)  by  virtue  of  their  large  third-order  nonlinear 
coefficients.  The  required  value  of  the  third  order  nonlinear 
susceptibility,  for  efficient  DFWM  with  typical  cw  lasers 

is  >  10“^  cm^/erg  =  10~^  esu.  This  value  will  be  used  as  a 

standard  in  evaluating  candidate  organic  materials. 

In  analyzing  the  nonlinear  optical  properties  of  organic 
materials,  several  questions  should  be  considered.  The  first  is, 

"Do  solid  and  liquid  organics  that  have  large  values  of  the  second- 
order  nonlinear  susceptibility  also  have  large  values  of 

the  third-order  nonlinear  susceptibility  A  literature 

search  revealed  that  the  answer  is  yes,  and  that  tiiis  result  is  well 
known . 

The  second  question  is,  "Can  organics  with  large  x^'^  be 
predicted?"  Again,  the  answer  was  found  in  the  literature.  The 
requirements  to  obtain  large  x^^^  include  bonds,  long  straight 
molecules,  and  large  charge  displacement.  A  very  simple,  tree- 
electron  model  is  developed  to  illustrate  these  results. 


Additional  questions  t.tiat  arose  during  the  investigation  are 
follows:  What  jnaximurn  values  of  and  time  constant  are 

attainable?  The  preliminary  answer  is  that  the  limit  is  r;et  oy  tn 
conditions  =  L^,  =  v/m  and  , 

very  roughly,  where  Lgff  is  the  effective  molecular  length,  m 
is  the  optical  infrared,  or  microwave  frequency,  and  v  and  t  are  th 
electron  velocity  and  relaxation  time-  These  condition  appear  to 
preclude  the  possibility  of  attaining  exceptionally  large  values  o 
X^^^  in  the  visible  region,  but  large  values  of  are 

predicted  at  microwave  frequencies  below  100  GHz. 

Donors  and  acceptors,  which  increase  the  value  of 
also  increase  the  value  of  x^^^.  The  final  question  remains 
unanswered:  Are  the  individual  x^'^/  etc. 

meaningful  when  the  x's  have  such  large  values  as  x  ^ ^  =  10“" 
c:n^/erg?  A  different  approach,  not  treating  the  x  ^ ^  as 
small,  may  be  necessary. 

1 .  Physical  Description  of  Large  Susceptibilities 

As  illustrated  in  Figure  2-12,  large  dipole  moments  and  large 
susceptibilities  result  from  displacements  of  charge  over  great 
distances.  In  organic  molecules,  the  strongly  overlapping  tt  hon.ls 
afford  good  conduction  paths,  as  illustrated  in  Figure  2-13.  Thus 
long,  IT -bonded  molecules  have  large  susceptibilities. 

It  is  known  that  the  periodic  structure  of  molecules  that  are 
not  straight  impedes  the  electron  motion  and  slightly  reduces  the 
susceptibilities.  However,  this  effect  is  less  important  t’nan  n 
bonds  and  great  lengths. 


2 .  Free-Electron  Model 

A  simple  f ree-electron  model  illustrates  the  central  feat  ires 
of  the  large  nonlinear  susceptibilities.  The  molecule  is 
represented  by  a  box  of  length  L  with  perfectly  reflecting  walls 
(see  Figure  2-14).  The  applied  electric  fiel<l,  accelerntos 
the  electron  when  it  moves  in  the  direction  opposit -o  to 
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E  (force  F--|e|E  =  eE),  and  slows  it  down  when  it  mover3  in  the 
direction  of  E.  The  equation  of  motion  is 


in 


dv 

dt 


G  E 


1  I 


The  solution  is 


V  =  V  +  — ^  =  V  (  1  +  >5 )  ,  (12) 

o  in  o  '  '  ' 

where  6  =  e  Et/mv^,  Vq  =  (2e/m)^^^,  and  c  is  the  initial 
electron  energy.  The  probability  that  an  electron  is  between  x  ani 
x+dx  is  proportional  to  v“ ^ : 


P  dx 

X 


V  ^dx 
fL  -1  , 

J  V  dx 
o 


(13) 


The  dipole  moment  is  easily  calculated  from  using  the 
expression 


e  P  (x-L/2)dx 

ox 

(  M) 

aE  +  3E^  +  YE^  , 

(15) 

where  ot,  3  and  Y  are  known  as  hyperpolarizabilities.  The  induced 
polarization  is 


P  =  Ru/v  +  +  .  (16) 

where  R  is  a  form  factor,  and  V  is  the  volume,  proportional  to  L, 
the  length  of  the  box. 

To  continue  with  our  solution  for  m  ,  we  consider  ttie  Unit  oi 


small  E  (  <5 < <  1  )  : 


We  cat!  relate  t  to  z  by  integrating  the  equation  of  motion 


for  V : 


/  vdt 
o 


=  vot  +(eE/2m)t‘ 


We  thus  find  that 

«  =  ^  =  -I  +  (W  aex|)1/2  _  ,,,, 

o  mv 

o 

Finally,  by  combining  Equations  (13),  (14),  (15),  (17),  and  (19),  we 
find  that 


1  4,  4 

1  e  L 


=  3.6  X  10  e  ^  cm^/erg. 


where  r.  is  in  A  and  e  is  in  eV.  Using  typical  values  for  the 
local  field  correlation  factor,  and  taking  into  account  the  linear 
dependence  of  molecular  volume  on  length,  we  find  that 

~  3.6  X  10  cm^/erg  .  (21) 

Experiments  by  Hermann  and  Queuing  have  confirmed  the  predicted  L** 
dependence  of  Y.  Molecules  with  as  many  as  30  double  bonds  were 
considered,  corresponding  to  Y  ”  10“^^  cm^/erg.  The  model  also 
gives  the  correct  order  of  magnitude  for  Y  if  e  »  )  eV  is  used. 

3 .  Giant  Values  of  Nonlinear  Suscepibilities 

The  formal  use  of  the  above  model,  which  gives  such  good 
agreoiTient  with  experimental  values  for  molecules  up  to  ~20  bonds  in 
length,  predicts  enormous  values  of  for  sufficiently  long 

molecules.  For  example,  we  find  from  Equation  (18)  that  a  length  of 
L  =  2.3  Mm  gives  rise  to  x^ =  10“^  cm^/erg. 
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Intuitively,  such  giant  values  of  seem  unlikely.  In 

order  to  determine  if  these  are  fundamentally  possible,  we  ask  that 
the  practical  limit  is  to  L.  Polymer  physicists  believe  that 
molecules  3  pm  long  are  possible;  so  the  physical  length  is  not 
believed  to  be  the  limit. 

The  first  limit  is  the  electron  mean  free  path,  ^rnfp*  If 
A[„fp<<L,  then  the  velocity  is  v  =  pE  (where  the  mobility  is 
p  =  et/m,  and  t  is  the  electron  relaxation  time  for  moinentu.n 
reversal),  rather  than  v  =  Vq  +  eEt/m,  and  given  in  Kcjuation  (12). 
In  this  case  of  the  velocity  is  nearly  constant  so 

that  the  probability  is  nearly  constant,  and  p  is  small. 

As  an  estimate  of  the  value  of  tor  organics,  use 


With 


V  =  (2  e/m)^/^  s  lo®  cm/S 


for  free  electron  and 


1 


mp/e 

2.8  :<  10 


(23) 


for  P  =  500  cm^/Vs  =  1.5x10^  cm^/statvolts ,  Eijuation  (22)  gives 


A  ~  0.28  pm 
mfp 

With  L  =  Aj^£p  =  0.28  Pm,  we  find  that 

-  (10  ^  cm^/erg ) ( 0 . 28  pm/2.8  pm) 
^  10  ^  cm^/erg  , 


(24) 


as  a  very  rough  estimate. 


Next,  consider  the  frequency  restriction.  The  distance  the 
electron  travels  before  the  electric  field  reverse  is,  roughly,  for 
A  =  0.53  urn. 


=  v/w 


=  2.Sx  10  ^cm=  2.8x  10 


Formally  setting  L  =  Lu  =  2.8  x  10  Pm  gives 


=  10  cm^/erg  , 


which  is  not  large  because  L  =  2.8  A  is  short.  This  example 
illustrates  that  giant  values  of  are  precluded  at  optical 

frequencies  by  the  frequency  limitation. 

For  a  microwave  frequency  of,  say,  10  GHz,  the  value  of  Lu) 


L, ,  =  1  . 6  X  10  cm  =  16  pm 


.'lit'i  L  =  Lf.)  =  16  pm,  we  find  that 


(3)  ,  o  3  , 

X  =  1 . 9  cm  /erg 


This  large  value  of  x^^^  indicates  that  the  frequency  limitation 
for  x^^^  in  the  microwave  region  is  not  the  final  limitation. 

The  mean-f ree-path  liinitacion  restricts  the  value  of  x^^^  in  the 
microwave  region  to  x^‘^^  ^  10“^  cm^/erg,  according  to  Equation 
(24),  which  is  a  very  rough  estimate.  These  very  rough  estimates 
suggest  ttiat  values  of  x^^^  -  10“^  cm^/erg  may  be  attainable  at 


microwave  frequencies, 


,*,**w'  •  •  • 
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SECTION  3 


IMAGE  TRANSMISSION  THROUGH  FIBERS  USING  PHASE  CONJUGATION 
AND  THE  USE  OF  FIBERS  AS  THE  PHASE  CONJUGATOR 

In  this  portion  of  the  program  we  deinonstated  the 
tifansmisSion  of  images  through  optical  fibers  using  phase 
conjugation.  We  also  investigated  the  suitability  of 
Commel?cially  available  fibers  to  serve  as  the  nonlinear  medium. 

A.  TRANSMISSION  OF  IMAGES  THROUGH  FIBERS  USING  PHASE 

CONJUGATION 

During  this  portion  of  the  program  we  demonstrated  the 
transmission  of  images  through  multimode  fiber  waveguides  using 
optical  phase  conjugation^^  or  wavefront  reversal  techniques. 

The  transmission  of  spatial  information  through  fibers  is 
degraded  due  to  mode  mixing  and  modal  dispersion.  By  making  use 
of  the  wavefront  reversal  property  of  nonlinear  optical 
techniques'^  we  can  compensate  for  these  effects.  In  our 
preliminary  experiments  an  image  was  transmitted  through  a 
fiber,  conjugated,  and  made  to  retraverse  the  same  fiber, 
resulting  in  the  reconstruction  of  the  image  at  the  input  end  of 
the  fiber.  For  this  to  be  a  useful  forward  transmission 
technique,  operation  with  two  nearly  identical  fibers  will  be 
required.  The  image  will  be  transmitted  through  the  first 
fiber,  conjugated,  and  then  transmitted  through  the  second  fiber 
which  reconstructs  the  image  at  its  output. 

The  preliminary  demonstration^  was  accomplished  by 
utilizing  a  multimode  (85  urn  diameter,  1.75  m  long)  step-index 
fiber.  The  results  of  this  experiment  were  published  in  Optics 
Letters  1_,  No.  11,  p.  558,  in  an  article  entitled  "Demonstration 
of  image  transmission  through  fibers  by  optical  phase 
conjugation".  The  manuscript  is  included  as  Appendix  E.  The 
output  of  an  argon  laser  (5145  A)  was  transmitted  through  a 
portion  of  an  Air  Force  resolution  chart  and  imaged  at  the  input 
end  of  the  fiber.  The  fiber  output  was  directed  into  a  crystal 
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of  barium  titanate  (BaTi03)  which  DFWM  was  used  to  generate 
the  phase  conjugate  signal.  The  phase  conjugate  signal  then 
retraversed  the  same  fiber  and  reconstructed  the  image  of  the 
resolution  chart.  The  results  of  this  experiment  are  shown  in 
Figure  3-1.  The  first  photograph  shows  the  input  image  of  an 
^ir  Force  resolution  chart.  The  second  photograph  shows  the 
output  at  the  end  of  the  fiber.  The  output  is  typical  of  a 
highly  multimode  fiber  where  a  large  number  of  modes  are 
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F igure  3-1 . 

Demonstration  of  image 
transmission  through  a 
fiber.  (a)  Image  of 
Air  Force  resolution 
chart  at  input  end  of 
fiber.  (b)  Output  after 
a  single  pass  through  the 
fiber.  (c)  Reconstruction 
of  resolution  chart  after 
output  has  been  phase 
conjugated  and  retraverses 
fiber. 


excited.  From  the  output,  it  io  impossible  to  discern  tnat  uny 
spatial  information  has  been  transmitted  throuyh  the  fiber.  L'lie 
reconstructed  image  after  double  passing  the  fiber  is  shown  in 
the  last  photograph. 

Geometry  employed  for  the  four-wave  mixing  experiments  in 
this  work  is  shown  in  Figure  2-3.  It  consists  of  two 
counterpropagating  pumps  (Ef,h’^j)  and  a  probe  wave  (hp) 

(which  is  the  beam  with  the  image  information  imposed  on  it) 
incident  on  the  nonlinear  medium.  These  beams  interact  to 
generate  the  phase  conjugate  signal  (Eg). 

A  schematic  of  the  experiment  is  shown  in  Figure  3-2.  and 
the  experimental  apparatus  is  shown  in  Figure  3-3.  'I'he  argon- 
ion  laser  (5145  A)  used  as  the  illumination  source  has  an 
intracavity  prism  and  etalon  for  single  longitudinal  mode 
operation  giving  coherence  lengths  on  the  order  of  10  m.  In 
order  to  have  useful  interference  between  the  pump  and  probe  the 
coherence  length  of  the  laser  must  be  comparable  to  the  pat'n 
length  differences  encountered  in  the  experiment.  For  our 
experiments  the  probe  path  length  is  longer  than  the  path  length 
of  the  pump  beams  by  at  least  the  length  of  the  fiber,  thus 
requiring  the  ~10  m  coherence  length. 

The  laser  output  is  continually  monitored  by  a  scanning 
Fabry-Perot  etalon.  Because  of  crystal  orientation,  p-polariza- 
tion  of  the  input  beams  is  required  in  our  experiments.  There¬ 
fore,  a  A/2  plate  and  Faraday  rotator  were  used  to  rotate  the 
plane  of  polarization  of  the  argon  laser.  The  Faraday  rotator 
was  also  used  in  conjunction  with  a  polarizer  to  prevent  any 
return  beams  from  entering  the  laser  and  destroying  single  mode 
operation.  The  beam  splitter,  BS2,  picks  oft  the  pump  beams.  A 
100  cm  f.l.  lens  was  used  to  reduce  the  pump  beam  size  in  the 
crystal.  The  beam  splitter,  BS3,  picked  off  a  s.econdary  proho 
signal  which  was  used  to  optimize  the  crystal  orientation.  Thn 
beam  splitter,  BS4,  was  used  to  generate  the  indivitJual 
counterpropagating  pumps.  The  pumpr.  typurally  nave  a  povv'o' 
ratio  of  2:1.  The  throughput  <if  the  boain  splltcer,  i1S2,  i:-.  n-.e.l 
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to  obtain  the  probe  used  in  the  fiber.  This  probe  beam  was 
expanded  by  a  beam  expanding  telescope  to  illuminate  the  Air 
Force  resolution  chart.  Following  the  resolution  chart  is  a 
beam  splitter,  BS6.  This  element  is  used  to  pick  off  the 
return  signal  which  forms  an  image  on  a  sheet  of  cardboard.  The 
image  plane  is  the  exact  same  distance  from  the  beam  splitter  as 
the  original  Air  Force  chart.  An  f.l.,  7.6  cm  focal  length 
relay  lens  was  used  to  capture  the  light  diffracted  by  the 
patterns  of  the  resolution  chart.  This  diffracted  light  would 
otherwise  miss  the  entrance  aperture  at  the  microscopic 
objective.  The  microscope  objectives  are  7X,  0.20 
N.A.  The  first  microscope  objective  couples  the  light  into  the 
fiber.  The  fibers  are  held  by  a  unique  Hughes-designed  fiber 
holder  and  fiber  assembly  which  has  five  degrees  of  freedom: 

X  -  y  -  z  translation  and  azimuthal  and  polar  rotations. 

In  our  studies  we  have  used  four  types  of  fibers  which  are 
representative  of  the  various  fiber  designs.  The  first  fiber  is 
a  "single"  or  lower  mode  number  fiber  from  Bell  r_,aborator ies . 

It  has  a  step  index  with  a  core  diameter  of  8.5  urn  and  a 
numerical  aperture,  N.A.,  of  0.052.  Tliere  are  several  lengths, 
the  longest  being  60  m.  At  the  wavelengths  of  interest,  ~5000 
we  can  propagate  three  discrete  modes,  of  wiiich  two  are 
degenerate.  The  number  of  modes,  M,  is  given  by  M  = 
k '^-a  ^  (  N .  A .  )  ^  ,  where  a  is  the  core  diameter.  The  second  fiber  is 
a  multimode  graded  index  fiber  supplied  by  IT']’.  The  core 
diameter  is  50  urn,  the  N.A.  is  0.25  and  the  length  is  6  m. 
Theoretically,  approximately  10"*  modes  can  propagate  in  this 
fiber  at  the  wavelengths  of  interest.  The  third  fiber  is  a 
Corning  85  urn  step  index  fiber  (N.A.  =  0.134)  which  is  1.75  m  in 
length.  It  also  can  propagate  ~10"  modes.  The  fourth  fiber  is 
a  inultimode  step  index  fiber  supplied  by  Quartz  Products.  The 
core  diameter  is  200  urn  and  the  N.A.  is  0.20.  This  Fiber  can 
su[5port  on  the  order  of  10  '  mo'les.  We  have  measured  ttie 
properties  of  fiber  transmission  at  three  di'derent  wavelengths, 
4131  C' ,  5145  ,  and  6471  A  so  that  we  can  ooti  ni/o  both  the 


wavelength  and  fiber  for  a  particular  experiment.  As  mentioned 
for  the  initial  image  transmission  studies,  we  used  the  Corning 
85  urn  fiber. 

A  second  microscope  objective  (see  Figure  3-2) 
quasicollimates  the  output  beam  from  the  fiber.  A  5  cm  focal 
length  lens  is  then  used  to  reduce  the  probe  beam  to  a  diameter 
smaller  than  the  diameter  of  the  pump  beams  inside  the  crystal. 
The  pumps  typically  make  an  angle  of  75°  with  respect  to  the 
c-axis  and  the  probe,  which  has  an  angle  of  68°.  This  implies 
that  the  angle  of  the  grating  normal  is  18°  with  respect  to  the 
c-axis.  The  highest  reflectivity  is  obtained  when  the  grating 
vector  is  inclined  to  the  positive  c-axis  direction,  as 
predicted  by  Feinberg  and  Hellwarth.^  Using  this  geometry  we 
have  obtained  reflectivities  in  excess  of  133%. 

These  large  returns  have  enabled  us  to  easily  reconstruct 
the  image  and  show  image  amplif iction.  At  present  we  have  been 
able  to  resolve  better  than  15  lines/mm  (67  urn  spacing),  see 
Figure  3-4.  (The  resolution  of  the  optical  system  without  the 
fiber  is  >20  lines/mm;  and  we  have  seen  in  a  separate  high 
resolution  set-up  greater  than  114  lines/mm  (9  urn  spacing)  in 
BaTi03  by  itself.)  The  imaging  geometry  used  for  these 
preliminary  demonstrations  is  not  optimized.  In  addition, 
effects  such  as  damping  of  the  higher  order  modes  relative  to 
the  lower  order  modes  (which  could  cause  a  scrambling  or  loss  of 
informtion)  and  polarization  scrambing  can  effect  the  image 
quality.  In  the  latter  case  typical  commercial  fibers  do  not 
preserve  polarization.  Hence,  the  linearly  polarized  input 
probe  beam  will  emerge  randomly  polarized.  For  our  DFWM 
experiments  the  pumps  are  linearly  polarized  in  the  plane  of  the 
experiment  and  only  this  component  of  polarization  is 
conjugated . 
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Reconstructed  image  after  double  passing 
the  fiber.  These  experiments  shov/  the 
ri.esolution  is  in  excess  of  IS  lines/mm. 


B.  PHASE  CONJUGATION  USING  OPTICAL  FIBERS  AS  THE  NON-LINEAR 

MEDIUM 

Because  of  the  possible  systems  applications  it  was  of 
interest  to  see  if  a  fiber  could  be  used  as  the  nonlinear  aiedium 
itself.  The  experimental  arrangement  for  the  DFWM  experiments 
using  fused  silica  fibers  as  the  phase  conjugate  medium  is  shown 
in  Figure  3-5.  The  forward  pump  and  probe  were  brought  in 
parallel  to  each  other  and  were  symmetrically  displaced  from  the 
optic  axis  of  the  first  lens.  This  lens  was  designed  to  couple 
the  forward  pump  and  probe  into  the  low  order  modes  of  the 
fiber.  The  backward  pump  was  coupled  into  the  other  end  of  the 
fiber  by  a  microscope  objective.  The  backward  pump  was  also 
aligned  to  launch  only  low  order  inodes.  The  fibers  used  were 
typically  5  m  in  length  and  have  been  discussed  in  Section  3A. 

In  performing  these  tests  we  encountered  difficulties  in 
determining  if,  in  fact,  a  DFWM  signal  was  observed  because  of 
s ignal-to-noise  problems  fundamental  to  a  DFWM  experiment  in 
fibers.  There  are  two  basic  noise  sources,  one  arising  from  the 
backward  traveling  pump  and  the  other  due  to  Fresnel  reflection 
of  the  probe  from  the  input  end  of  the  fiber.  Both  signals 
limit  the  DFWM  reflectivity  that  can  be  detected  and  are  larger 
than  the  expected  signal  based  on  a  simple  calculation  (see 
below).  (To  understand  these  noise  sources  refer  to  Figure  3-5.) 
The  experimental  arrangement  allowed  us,  in  principle,  to 
spatially  separate  the  forward  pump,  backward  purap  and  probe. 

The  signal  was  separated  from  the  probe  path  by  the  signal-beam 
splitter.  There  were  two  competing  design  considerations. 

For  best  coupling  it  was  necessary  for  the  forward  pump  and 
probe  to  be  physically  collinear  with  the  optic  axis  of  the 
first  lens.  But  for  good  phase  conjugate  signal  separation  the 
probe  should  be  as  far  from  the  optic  axis  as  is  physically 
possible;  i.e.,  as  close  to  the  edge  of  the  lens  aperture  as 


gure  3  5.  Schematic  of  degenerate  four  wave  mixing  experiments 
using  silica  fiber  as  nonlinear  medium. 


possible.  In  our  preliminary  work  we  attempted  a  compromise 
between  these  conditions.  Unfortunately,  because  oE  mode 
coupling  in  the  fiber  there  was  enough  energy  coupled  from  the 
backward  pump  into  higher  order  modes  that  the  energy  in  these 
modes  traveled  back  along  the  same  path  as  the  signal.  Phase 
lock  detection  was  used  (while  chopping  the  probe),  but  the 
scattered  backward  pump  prevented  a  conclusive  measurement.  The 
second  source  was  the  probe  retro-reflection  off  the  end  of  the 
fiber  and  coupling  optics.  We  were  able  to  redu.'e  this  retro- 
reflection  by  incorporating  a  polarizer/quarter  wave  plate 
combination.  The  linearly  polarized  light  from  the  polarizer 
passed  through  the  ^^4  plate  at  45°,  with  respect  to  the 
ordinary  and  extraordinary  axis.  The  emergent  light  was  then 
circularly  polarized  and  impinged  on  the  fiber  end.  The 
reflected  light  was  also  circularly  polarized.  The  light  then 
retraversed  the  ^/4  plate  and  emerged  linearly  polarized  but  now 
perpendicular  to  the  input  polarization  and  was  attenuated  by 
the  polarizer.  The  signal  polarization  is  not  preserved  in  the 
fiber  and  emerged  with  a  random  polarization  and  the  component 
parallel  to  the  polarizer  was  passed.  Even  with  this  solution 
to  the  retro  reflection  problem,  we  were  unable  to  detect  a  DFWil 
signal,  which  is  not  surprising  if  one  does  a  simple  calculation 
of  the  expected  conjugate  signal  strength,  as  given  below. 

For  a  nonresonant  material  the  phase  conjugate  reflectivity 
is  given  by 

R  =  tanM  K|  ,  (5) 


where 


For  silica  glass  the  linear  and  nonlinear  indices  of  refraction 
are; 

and 

n2  =  1.8  X  10"^^  esu, 

giving 

=  4.3  ,  10-14  ^3^  ^ 

For  fiber  diameters  of  100  um  and  for  an  input  powr  of  1  W,  I  = 

2  X  10^.  Thus  , f or  a  wavelength  of  0.5  pm  and  a  fiber  length  of 
100  cm,  R  =  10“^%.  In  order  to  increase  this  reflectivity  one 
should  dope  the  fiber,  thus  increasing  the  nonlinear 
susceptibility.  The  use  of  optical  fibers  as  the  nonlinear 
medium  for  phase  conjugation  will  require  the  development  of 
fibers  with  a  higher  nonlinear  susceptibility  than  are  presently 
available . 


SECTION  4 


TWO-WAVE  MIXING 

One  technique  for  studying  the  photoref ractive  effect  is  to 
examine  the  energy  coupling  between  a  probe  beam  and  a  pump  beam  in 
a  two-wave  mixing  experiment*.  The  probe  and  pump  interfere  in  the 
crystal,  forming  a  fringe  pattern  throughout  the  volume.  These 
fields  cause  a  migration  and  separation  of  changes  in  the  crystal, 
thereby  creating  an  internal  space  charge  electric  field.  This 
field  produces  a  change  in  the  refractive  index  by  the  electro- 
optic  effect  and  creates  a  volume  phase  hologram.  In  photo- 
refractive  materials  the  refractive  index  grating  and  the  two  beam 
interference  pattern  are  shifted  in  phase  spatically.  The  phase 
shift  for  BaTi03  is  typically  90“  out  of  phase;  this  causes  energy 
to  be  transferred  betwen  the  pump  and  probe  beams.  (See  Figure  2  in 
Appendix  C). 

An  experiment  to  investigate  beam  coupling  in  BaTi03  was 
conducted,  a  schematic  of  which  is  shown  in  Figure  4-1.  In  the 
probe  path  were  two  acousto-optic  (A.O. )  modulators.  By  modulating 
the  carrier  frequency  of  the  first  A.O.  modulator  it  is  possible  to 
generate  a  temporal  variation  in  intensity  of  the  probe  beam.  The 
diffracted  order  is  a  replica  of  the  waveform  of  the  envelope  of  the 
acoustic  carrier  frequency.  The  diffracted  order,  however,  is 
shifted  from  the  original  input  laser  frequency  by  the  acoustic 
carrier  frequency.  Because  useful  interference  in  the  BaTi03 
requires  that  the  probe  and  pump  must  be  at  the  same  frequency,  a 
second  unmodulated  A.O.  cell  is  used  to  shift  the  probe  frequency 
back  to  the  original  laser  frequency. 

After  diffracting  from  the  A.O.  modulators,  the  probe  beam 
propagates  to  the  BaTi03  crystal  where  it  interacts  with  the  pump 
beam.  In  these  experiments  the  pump  is  cw  and  is  unmodulated.  The 
pump  flux  is  typically  ten  to  one  hundred  times  the  probe  flux. 

*This  work  was  supported  by  IR&D  funding  and  is  included  here  due  to 
its  relevence  to  this  program. 
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figure  4-1. 


Two-beam  coupling  experiment  schematic. 


After  experiencing  amplification  in  the  crystal,  the  transmitted 
probe  beam  is  detected  by  a  photomultiplier  tube.  The  signal  is 
then  displayed  on  either  an  oscilloscope  or  a  chart  recorder. 

Three  types  of  modulation  on  the  probe  beam  were  investigated. 
The  first  is  a  single  pulse  event,  the  second  is  a  continuous 
unmodulated  beam,  and  the  third  is  a  repetitive  waveform.  An 
example  of  the  single  pulse  event  is  shown  in  Figure  4-2a.  The 
lower  trace  shows  the  transmitted  probe  in  the  absence  of  the  pum£) 
beam.  It  is  seen  to  be  a  square  wave  of  20  mV  amplitude  and 
1300  msec  duration.  The  upper  trace  is  the  transmitted  probe  when 
the  pump  beam  is  present.  Because  of  two-beam  coupling,  some  of 
the  energy  of  the  pump  beam  is  diffracted  in  the  direction  of  the 
transmitted  probe.  This  results  in  amplification  of  the  transmitted 
probe.  In  Figure  4-2a  the  amplification  is  a  factor  of  twenty-three 
at  the  highest  value.  This  result  was  obtained  for  an  internal 
pump-probe  angle  of  3.5“  and  an  angle  between  the  c-axis  and  grating 
normal  of  16.8“.  The  power  in  the  pump  was  33  mW ,  while  the  probe 
power  was  83  wW. 

The  rise  time  of  the  transmitted  probe  pulse  is  related  to  the 
time  constant  associated  with  the  grating  formation.  The  flat 
portion  is  equal  to  the  steady-state  value  when  the  probe  is 
unmodulated*  The  decay  of  the  transmitted  probe  is  characterized  by 
a  very  fast  initial  change,  followed  by  a  slower  decay  associated 
with  grating  erasure.  It  is  believed  that  the  fast  decay  is  a 
coherent  effect  due  to  removal  of  the  input  probe  beam.  If  we 
consider  the  grating  to  have  an  amplitude  reflectivity  n,  and  we 
assume  that  the  transmitted  probe  results  from  the  coherent  addition 
of  the  diffracted  pump  and  the  undiffracted  probe,  then  the 
transmitted  probe  intensity  can  be  approximated  as 


^probe  "  ’’'^ump 

+  n^( 1-n)^  (l  I  h  J  • 

pump  probe 


When  the  probe  is  turned  off,  the  first  and  third  terms  go 
to  zero  in  the  probe  switching  time  (200  nsec).  The  difference 
in  magnitude  between  the  probe  signal  when  switched  off  (lower 
trace)  and  the  amplified  signal  when  switched  off  (upper  trace) 
is  the  contribution  of  the  third  term.  Affer  the  probe  is  turned 
off  the  grating  which  was  formed  continues  to  be  read  by  the  pump 
beam  and  energy  is  diffracted  in  the  transmitted  probe  direction. 
As  the  pump  beam  is  reading  the  grating  it  is  simultaneously 
erasing  the  grating <  Thus,  the  contribution  from  the  second  term 
decays  at  a  rate  characteristic  of  the  refractive  index  grating 
erasure  time  in  photoref ractive  materials.  This  is  demonstrated 
by  the  exponentially  decaying  tail  (x  =  0.35  sec). 

If  the  pulse  length  of  the  probe  beam  is  shorter  than  the 
write  time  (typically  on  the  order  of  the  decay  time,  350  msec), 
the  transmitted  probe  will  not  reach  the  maximum  value  obtained 
in  Figure  4“2a*  The  transmitted  probe  will  build  up  with  the 
same  rise  time  but  will  not  reach  the  steady-state  value,  as 
illustrated  in  Figure  4-2c.  (Note  that  the  vertical  scale  in 
Figure  4“2a  and  4-2b  are  100  mv/division,  and  for  50-2c  the 
vertical  scale  is  50  mV/division.)  The  magnitude  of  the 
discontinuity  is  independent  of  the  pulse  length.  The  decay  time 
of  the  trailing  edge  is  also  seen  to  be  independent  of  pulse 
length . 

The  theory  presented  in  Section  1  and  Appendix  D  describes 
the  charge  dertsities  generated  for  a  pulse  probe  beam  and  cw  pump 
beami  The  amplified  probe  beam  will  qualitatively  follow  the 
charge  density*  The  computer  generated  graphs  in  Section  1 , 
Figures  4,  5  and  6  are  for  probe  pulse  greater  than,  equal  to, 
and  less  than  the  dielectric  relaxation  time.  The  theory  is 
seen  to  agree  well  with  the  data  of  Figures  4-2(a),  4-2(b)  and 


The  second  modulation  regime  occurs  when  the  incident  probe 
is  unmodulated.  When  the  probe  is  first  tuned  on,  the  two  wave 
coupling  signal  rises  to  a  constant  "saturated"  value.  This 
value  is  represented  by  curve  A  in  Figure  4-3.  This  signal  level 
will  be  refered  to  as  the  cw  level.  This  is  the  same  level 
obtained  by  a  single  pulse  if  the  pulse  duration  is  greater  than 
the  write  and  erase  times. 

The  third  regime  occurs  when  the  modulation  of  the  probe  is 
a  repetitive  waveform.  In  this  case  the  transmitted  probe  signal 
depends  upon  whether  the  period  of  the  waveform  is  greater  than 
or  less  than  the  dielectric  relaxation  time  { t ) .  if  the  period 
is  greater  than  x  the  depth  of  modulation  is  one  hundred  percent. 
This  case  is  shown  in  the  data  of  Figure  4-3,  curve  B.  The  signal 
rises  from  zero  with  the  characteristic  material  response  rise 
time  to  the  cw  steady-state  value.  When  the  probe  pulse  is 
switched  off,  the  diffracted  pump  energy  decay  back  down  to  the 
zero  point.  This  waveform  is  repeated  and  is  similar  to  a  con¬ 
tinuous  train  of  individual  pulse  response.  However,  when  the 
waveform  period  is  less  than  t  several  interesting  things  occur. 
These  are  illustrated  in  the  data  shown  in  Fiure  4-4a.  In  the 
first  cycle  the  transmitted  probe  (upper  trace)  is  superimposed 
over  that  of  a  single  pulse  (lower  trace).  However,  instead  of 
decaying  back  down  to  the  zero  level  it  has  a  finite  value  when 
the  next  increasing  portion  of  the  probe  waveform  occurs.  The 
signal  then  rises  until  the  probe  waveform  begins  to  decrease. 
This  increase  and  decrease  in  the  signal  finally  evolves  in  to  a 
stea' /-state  waveform  which  is  comprised  of  a  dc  level  added  to  a 
periodic  component  proportional  to  the  incident  probe  waveform. 

Further  examination  of  Figure  4-3  reveals  other  information 
about  two-wave  mixing.  It  is  evident  that  the  depth  of 
modulation  is  a  ^unction  of  the  period,  and  is  demonstrated  by 
comparin)  traces  iHT=10  sec),  C(T=1  sec),  and  D(T=0.1  sec).  The 
dc  cuinpanent  in  thf?  signal  tends  toward  the  cw  level  (curve  A) 

L  imM:',  tin,'  dut  y  cycle  (50'4).  This  was  verified  in  other 
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Figure  4-3. 


Repetitive  waveform  data  for  two-beam  couplin 
experiment.  A.  Probe  cw.  B.  Probe  0.1  Hz  sg 
wave.  C.  Probe  1.0  Hz  square  wave.  D.  Probe 


experiments,  the  results  of  which  are  not  shown.  These  three 
cases  are  clearly  illustrated  in  Figure  4-4b.  The  upper  trace 
(iii)  is  obtained  when  the  probe  is  cw.  The  middle  trace  (ii) 
results  from  a  continuous  waveform,  while  the  lower  trace  (i)  is 
the  result  of  a  single  pulse. 
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OPTICAL  CONVOLUTION  AND  CORRELATION 

This  portion  of  the  program  was  directed  toward  the 
demonstration  of  optical  convolution  and  correlation.  The  first 
section  addresses  static  two-dimensional  optical  processing. 

The  second  section  describes  efforts  toward  demonstrating  real- 
time  optical  correlation. 

A.  STATIC  TWO  DIMENSIONAL  OPTICAL  CONVOLUTION  AND  CORRELA  I,’ I ( )N 

The  major  thrust  of  these  experiments  was  the  generation  of 
convolutions  and  correlations  of  two  spatially  modulated  input 
waves*  The  geometry  of  the  experiment  is  shown  in  Figure  5-1. 

The  geometry  is  basically  the  standard  degenerate  four-wave 
mixing  (DEV^M)  scheme,  consisting  of  two  counterpropagatinq  pumps 
and  a  probe.  However,  amplitude  and/or  phase  information  are 
now  encoded  on  the  beams.  The  phase  conjugator  is  placed  at 
the  common  back  focal  plane  of  lenses  LI,  L2  ami  L3.  TIi  .> 
at  this  position  are  proportional  to  the  Fourier  transfor;i\s  tjf 
the  inputs  at  input  planes  1,  2  and  3.  The  phase  conjugate 
signal  monitored  at  the  viewing  plane  is  proportional  to  tiie 
convolution  or  correlation  of  the  inputs  under  the  following 
conditions.  The  correlation  operation  is  performed  by  placing 
information  at  input  planes  1  and  2  and  a  point  source  at  plane  3. 
The  convolution  of  input  2  aaC  3  is  generated  if  a  delta 
function  (spatial)  exists  at  input  plane  1. 

Preliminary  correlation  experiments  have  used  naTi03  as 
the  phase  conjugate  medium.  These  experiments  are  similar  to 
those  reported  by  White  and  Yariv.^^  However,  in  our 
experiments  we  are  using  a  specially  generated  character  mask. 

The  characters  are  shown  in  Figure  5-2.  These  ctiaracfers 
chosen  for  ease  evaluation  of  the  two-d i mens iona 1  cross¬ 
correlation  function  defined  as. 


'faB(v)  =  /a(  11  )  (  M+V  )  dM, 
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Figure  5-2.  Character  mask  used  in  2-D  correlation  experiment, 
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0  F  □  D  D  D  D 

F  F  F  F 

□  □  □  □ 


5-2.  Character  mask  used  in  2-D  correlation  experiment 


where  u  =  (M1/U2)  are  the  two-dimensional  vector  coordinates 
and  V  =  ( 1  f  ''2)  two-dimensional  displacement 

coordinates.  The  four  characters  chosen  represent  values  of 
high  correlation,  as  between  the  "C"  and  the  "O" ,  and  relatively 
low  values,  as  between  the  "F"  and  the  "D".  The  normalized 
correlation  values  defined  below  are  given  in  Figure  5-3: 


ap 


a3 


aa  ^33 


A  schematic  of  the  experiment  is  shown  in  Figure  5-4.  A 
Krypton  ion  laser  (6471  A)  has  an  oven-controlled  etalon  for 
single  longitudinal  mode  operation.  The  laser  output  is 
controlled  and  monitored  by  a  scanning  Fabry-Perot  etalon. 
Because  of  orientation,  p-polarization  of  the  input  beams  is 
required  in  our  experiments.  Therefore,  a  A/2  plate  and  Faraday 
rotator  were  used  to  rotate  the  plane  of  polarization  of  the 
Krypton  laser.  The  Faraday  rotator  was  also  used  in  conjunction 
with  a  polarizer  to  prevent  any  return  beams  from  entering  the 
laser  and  destroying  single-mode  operation.  A  beam-expanding 
telescope  and  pinhole  were  used  to  increase  the  beam  diameter  to 
fully  illuiiiinate  the  character  matrix  and  to  spacially  filter 
the  beams.  Beam  splitter  1  reflects  a  portion  of  the  beam  which 
is  used  for  the  forward  pump.  Beam  splitter  2  reflects  the 
probe  beam,  while  the  transmitted  beam  is  used  as  the  backward 
pump.  Lenses  1,  2  and  3  have  a  f=50  cm  focal  length,  with  the 
crystal  placed  at  the  Fourier  transform  plane  of  these  lenses. 
Lach  mask  is  placed  a  distance  f  away  from  the  appropriate  lens. 
The  mask  was  ion-f^eam  etched  in  gold  film  on  glass,  and  the 
characters  are  one  millimeter  in  height.  By  illuminating  one 
oi  the  characters  of  mask  2  and  the  entire  matrix  of  mask  1,  the 
phase  conjugate  signal  will  be  a  matrix.  The  intensity  of  each 
element  corresponds  to  the  correlation  of  the  characters 
illuminated.  The  phase  conjugate  signal  is  directed  by  beam 
splitter  3  tc)  a  convenient  viewing  plane. 
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The  correlation  of  a  single  character  and  the  matrix  of 
characters  was  successfully  demonstrated.  However,  due  to  a 
redirection  of  the  objectives  of  this  contract  toward 
investigating  real-time  optical  correlation  further  work  was  not 
pursued . 

B.  REAL  TIME  CORRELATION 

This  portion  of  the  program  was  directed  toward 
demonstrating  the  temporal  correlation  of  signals  using  four- 
wave  mixing  configurations.  This  signal  processing  concept  was 
proposed  by  William  Micelli  of  RADC. 

The  scheme  uses  acousto-optic  modulators  to  impose  a 
traveling  replica,  of  a  finite  time  segment,  from  an  input 
electrical  signal  onto  the  transmitted  wavefront.  By  using 
barium  titanate  (a  photoref ractive  crystal)  as  the  phase 
conjugator,  the  degenerate  four-wave  mixing  (DFWM)  signal  should 
be  proportional  to  the  time  integrated  product  of  the  two 
electrical  signals. 

The  experiments  to  be  described  use  two  acousto-optic 
modulators  to  generate  the  function.  Pi ( t ) • F2 { t-t ) .  This  is 
done  by  using  an  acousto-optic  delay  line.  The  output 
intensities  and  frequencies  are  shown  in  Figure  5-5  where  the 
diffraction  efficiency  of  the  A.O.  modulator  is  a  and  the 
incident  probe  amplitude  is  A.  The  resulting  outputs  are  also 
shown  where  the  upshifted,  downshifted  or  undiffracted  beams  are 
taken  into  a  second  A.O.  modulator,  with  diffraction 
efficiency,  P.  It  is  seen  that  there  are  several  terms 
proportional  to  Fi(t)*F2(t).  However,  it  is  important  to 
note  that  in  order  for  interference  to  occur  in  the  phase 
conjugator  (which  is  BaTi03)  the  beams  with  the  correct 
temporal  information  must  be  at  the  same  frequency.  This  is 
because  of  the  slow  response  time  of  the  BaTi03.  If  two  beams 
at  different  frequencies  interfere,  a  moving  grating  results. 

If  the  material  cannot  respond  to  this  motion,  then  the  grating 
is  washed  out.  This  is  the  case  in  BaTi03 ,  since  its  response 


Figure  5-5.  Relative  indensity  and  frequency  of  diffracted 
orders  produced  by  acousto-optic  devices 


V 


time  is  a  fraction  of  a  second  (10  Hz),  whereas  the  frequency 
difference  is  40  MHz.  Therefore,  in  order  for  useful 
interference  to  occur  between  the  pumps  and  the  probe  in  the 
phase  conjugator,  they  must  be  at  the  same  frequency.  The 
diffracted  orders  are  shifted  from  the  zero  order  by  so 

care  must  be  taken  to  utilize  the  correct  orders.  (See  Figure 
5-5).  If  the  -1  order  from  the  acousto-optic  modulator  is  used 
then  the  +1  order  from  the  A.O.  delay  line  must  be  used. 
Similarly,  if  the  +1  order  from  the  A.O.  modulator  is  used,  thei 
the  -1  order  from  the  A.O.  delay  line  is  used.  Therefore,  only 
the  bold  face  terms  in  Figure  5-5  are  useful  in  these 
experiments . 

The  first  experiment  performed  using  this  A.O.  cell 
configuration  was  optical  multiplication  of  two  signals.  A 
schematic  of  the  experiment  is  given  in  Figure  5-6.  The  laser 
source  is  a  Krypton  ion  laser  (6471  A),  a  beam  compacting 
telescope  and  pinhole  were  used  to  reduce  the  beam  size  inside 
the  acousto-optic  modulator  and  to  spatially  filter  the  beam. 
Beam  splitter  1  transmits  the  backward  pump  beam,  while  the 
reflected  beam  enters  the  A.O.  modulator.  Both  the  A.O. 
modulator  and  A.O.  delay  line  are  driven  in  phase  at  40  MHz. 
I'he  schematic  of  the  electronics  is  shown  in  Figure  5-7 

A  master  oscillator  generates  the  40  MHz  carrier  frequency 
for  both  A.O.  devices.  The  output  is  split  into  two  branches 
by  a  power  divider.  The  signals  Fi(t)  and  F2(t)  are 
impressed  on  the  carrier  frequency  by  an  electronic  attenuator. 
The  modulated  signal  is  tlien  amplified  by  power  amplifiers  to  a 
level  sufficient  to  drive  the  A.O.  devices,  typically  4W . 

For  phase  matching,  tlie  functions  Fi(t)  and  2^^)  are 
constants.  When  Fi(t)  is  non-zero,  diffracted  orders  exist 
and  the  intensity  is  [jrupor  t  iona  1  to  the  magnitude  of  F-)  . 
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Electronics  used  to  drive  acousto-optic  devices 
bl(3ck  diagram. 


Returning  to  Figure  5-6,  the  diffracted  and  undi f f racte<5  beams 
are  recombined  (for  the  purpose  of  heterodyne  detection  to 
increase  dynamic  range)  using  beam  splitter  2.  The  function 
F2(t)  is  impressed  on  the  carrier  frequency  of  the  delay  line. 
For  these  experiments  the  beam  was  not  expanded  to  fill  the 
delay  line.  The  diffracted  and  undiffracted  orders  from  ttie 
A.O.  delay  line  are  interfered  at  the  phase  conjugator.  Tiie 
resulting  fringe  patterns  (shown  in  Figure  5-8)  were  made  by 
replacing  the  phase  conjugator  with  a  microscope  objective  and 
focusing  the  magnified  image  on  a  film  plane  with  F-;  and  F2 
constant.  To  obtain  the  magnification  of  the  fringe  pattern, 

F-]  was  set  equal  to  zero,  with  F2  equal  to  a  constant;  a 
resolution  chart  was  placed  in  front  of  the  microscope 
objective.  Keeping  the  objective  in  place,  the  chart  was  then 
moved  until  its  image  was  in  focus  at  the  film  plane.  By 
measuring  the  spacing  on  film  between  known  bar  patterns  of  the 
resolution  chart,  the  exact  magnification  of  the  system  could  be 
determined.  The  measured  values  of  fringe  spacing,  0,  agreed 
well  with  the  theoretical  values  predicted  by  D  =  A/(2  siri'j/2), 
where  9  is  the  angle  between  the  pump  and  probe  beams. 

In  Figure  5-6,  the  phase  conjugate  signal  retraces  the 
probe  path  and  is  reflected  by  beam  splitter  3  into  a 
photomultiplier.  The  resultant  signal  is  recorded  by  a  chart 
recorder.  Typical  data  is  shown  in  Figures  5-9A  and  5-9B.  In 
Figure  5-9A  the  function  Fi(t)  is  a  sawtooth  function  of 
.01  Hz.  ^2^^)  3  square  wave  of  0.9  Hz.  Notice  tnat  the 

resultant  signal  waveform  is  qualitatively  similar  to  tiie 
product  of  Fi  and  F2,  modified  by  the  time  response  of  tee 
phase  conjugator,  BaTi03.  In  Figure  5-9B,  F i ( t )  a  0.05  Mz 

square  wave  and  F2(t)  is  a  0.5  Hz  s^juare  wave.  ’I'iiii.i  iata 
clearly  shows  the  multiplication  of  f -|  (  t )  and  fjit). 


Observation  of  this  data  indicates  a  slowly  varying  peak 
amplitude  modulation  due  to  instabilities  in  the  optical 
mounts.  Because  BaTi03  has  a  long  time  response  (~400  msec), 
the  degenerate  four-wave  mixing  process  is  exactly  analogous  to 
making  a  conventional  hologram.  The  write  beam  and  reference 
beam  interfere  in  the  phase  conjugate  material,  forming  a 
grating  which  diffracts  a  portion  of  the  read  beam.  This 
diffracted  energy  is  the  four-wave  mixing  signal.  In  holography 
it  is  necessary  for  the  relative  path  length  differences  between 
the  write  and  reference  beams  to  be  constant  during  the 
exposure.  If  there  are  variations  on  the  order  of  a  half 
wavelength  the  hologram  loses  contrast  and  is  washed  out. 
Therefore,  it  is  important  in  holography  or  DFVJM  using 
photoref ract ive  crystals  for  the  optical  system  to  be 
mechanically  and  thermally  stable  for  extended  periods  of  time 
(minutes).  We  have  undergone  major  modifications  to  our  optics 
mounting  system  to  meet  these  requirements.  A  Michelson 
interferometer  was  used  to  quantify  the  degree  of  stability  for 
various  mounting  schemes.  A  detector  at  the  output  of  the 
interferometer  monitors  an  area  much  less  rhan  one  fringe 
spacing.  The  output  is  then  a  direct  measure  of  fringe  motion 
which  is  proportional  to  the  stability  of  the  optical  system. 

Data  is  shown  in  Figure  5-10  which  charts  the  detector  output 
versus  time  for  three  different  mounting  techniques.  There  are 
typically  ten  or  more  mounts  in  an  experimental  layout. 

Therefore,  any  noise  contributed  by  a  single  mount  can  be 
amplified  because  of  long  optical  paths  and  the  contribution 
from  additional  mounts.  The  stability  problem  was  initially  so 
severe  that  the  phase  conjugate  signal  generated  by  cw  pumps  and 
probe  was  modulated  by  one  hundred  per  cent.  This  modulation 
was  because  the  inotion  caused  the  gratin'"  to  be  continually 
erased  ani]  rewritten.  The  new  mounting  technique  of  mechanically 
rigid,  ma;:s,  ive  post  holders  on  a  magnetic  base  has  significantly 
reduced  the  problem.  The  table  is  presently  suspended  b^ 


A  schematic  o£  aa  experiment  £or  performing  the  correlation 
of  the  functions  Fi(t)  and  using  acousto-optic  devices 

and  phase  conjugation  is  shown  in  Figure  5-11.  The  function, 

F^  (t),  is  impressed  on  the  carrier  frequency  of  the  acousto¬ 
optic  modulator.  When  Fi(t)  is  null,  only  the  undiffracted 
order  is  present,  as  shown  in  dashed  lines.  When  Fi(t)  is 
non-zero,  a  diffracted  order  exists  and  is  proportional  to  the 
magnitude  of  F-)  .  A  beam  expanding  telescope  using  cylindrical 
optics  is  used  to  expand  the  beam  to  "fill"  the  acousto-optic 
delay  line.  The  function,  F2(t)  is  impressed  onto  the  carrier 
frequency  of  the  delay  line.  The  output  from  the  delay  line  is 
given  by  the  product  of  the  input,  Fi(t),  times  F2(t-x/v), 
where  x  is  the  position  measured  perpendicular  to  the  beam 
propagation  direction  and  v  is  the  acoustic  velocity. 

A  lens  is  used  to  demagnify  and  image  the  end  of  the  A.O. 
delay  line  into  the  phase  conjugator.  Although  the  operation  of 
this  system  is  not  completely  understood  at  this  time,  it 
appears  that  the  conjugator  (BaTi03)  performs  as  a  large  array 
of  coherent  detectors  integrating  the  complete  spatial  pattern 
imposed  on  the  A.O.  delay  line.  Furthermore,  because  of  the 
Long  time  response  of  QaTi03  it  also  integrates  the  incoming 
probe  signal  in  time.  The  phase  conjugate  signal  which  is 
generated  retraces  the  probe  path  to  beam  splitter  3,  where  a 
portion  is  deflected  into  a  photomultiplier,  giving  a  signal 
proportional  to  tLie  correlation  of  Fi(t)  and  F2(t). 

The  main  difference  between  this  scheme  and  the 
multiplication  experiment  is  that  the  acousto  optic  modulators 
are  only  in  the  probe  path,  whereas  before,  they  were  in  both 
the  forward  pump  and  probe  paths.  This  new  configuration  has 
the  advantage  that  the  prol^e  beam  is  only  modulated  proportional 
F ^ ( t ) * ( t-^ )  in  intensity,  eliminating  any  cross  terms  (See 
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Figure  5-n.  Real  time  correlation  experiment  schematic, 
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For  these  experiments,  a  single  longitudinal  mode  argon  ion 
laser  (5145^)  was  used.  A  beam  compacting  telescope  and  pinhole 
were  used  to  reduce  the  beam  size  inside  the  acousto-optic 
modulator  and  to  spatially  filter  the  beam.  The  function, 

Fi(t),  is  impressed  on  the  carrier  frequency  (40  MHz)  of  the 
acousto-optic  modulator.  A  beam  expanding  telescope  using 
cylindrical  optics  is  used  to  expand  the  beam  to  "fill"  the 
acousto-optic  delay  line.  The  function,  F2(t),  is  impressed 
onto  the  carrier  frequency,  F(.(40  MHz),  of  the  delay  line. 

Both  the  A.O.  modulator  and  A.O.  delay  line  carrier  frequency 
are  adjusted  to  be  in  phase.  The  output  from  the  A.O.  delay 
line  is  the  product  of  Fi(t)  and  F2(t-x/v).  In  this 
experiment  a  cylindrical  beam  compaction  telescope  was  used  to 
reduce  the  probe  size  smaller  than  the  pump  beams  in  the 
interaction  region  of  the  crystal.  Typical  data  from  these 
experiments  is  shown  for  auto-correlation  in  Figure  5-12.  The 
top  trace  in  Figure  5-1 2A  is  the  waveform  applied  to  the  A.O. 
modulator,  while  the  lower  trace  in  5-12A  is  the  waveform 
applied  to  the  optical  delay  line.  Only  the  envelope  (f=  70  hz) 
1)1  the  40  Mhz  carrier  frequency  is  visible.  Figure  5-1 2B  is  the 
corresponding  degenerate  four-wave  mixing  signal.  This  data 
lo  not  appear  to  be  consistent  with  a  correlation,  although 
(i  le  to  time  limitations  a  complete  interpretation  has  not  been 
riad'.' . 

[t  is  important  to  recognize  that  these  are  the  first 
exper imentdi  results  of  this  type  and  an  additional  amount  of 
analysis  and  experimentation  are  required  before  a  complete 
understanding  can  be  acnieved. 
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Figure  5-12. 


Real  time  auto  correlation  (lata. 

(a)  Fi(t)  waveform  applied  to  acousto¬ 
optic  modulator  lower  trace,  F  ( t ' 
waveform  applied  to  acousto-optic  delay 
time  upper  trace. 

(b)  Resultant  obasa} 


con  )  H'  la  t  r  )r  s  inna  1  . 


SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

We  have  demonstrated  that  phase  conjugate  optics  is  a 
viable  approach  for  information  processing  and  image 
transmission  applications.  The  particular  phase  conjugation 
interaction  which  we  have  found  useful  is  degenerate  Eour-wave 
mixing.  Among  the  materials  and  mechanisms  available  for  OFWM, 
photoref ractive  materials  appear  most  promising.  These 
materials  have  high  efficiency  and  spatial  resolution  at  cw 
power  levels.  The  response  time  in  some  cases  is  slow,  but  our 
analysis  indicates  that  improved  control  of  dopant  valence  state 
and  concentration  should  lead  to  a  significant  increase  in 
speed. 

We  have  demonstrated  that  four-wave  mixing  can  be  used  to 
transmit  images  in  multimode  fibers.  In  our  experiments,  the 
phase-conjugate  image  was  re-transmitted  through  the  same  fiber, 
yielding  a  spatial  resolution  of  15  lines/mm  over  a  1.75  M 
length.  The  goals  of  any  future  effort  in  this  area  would  be  to 
use  improved  collection  optics  to  increase  the  spatial 
resolution,  and  to  use  a  separate  identical  fiber  allowing  for 
forward  transmission  of  the  image. 

We  have  utilized  the  ability  of  phase  conjugation  to 
perform  the  multiplication  of  two  or  three  input  signals  and 
have  demonstrated  optical  correlation  between  various  two- 
dimensional  characters.  In  addition,  we  have  investigated  the 
real-time  correlation  of  two  electrical  signals  by  first 
converting  the  signals  to  spatially  modulated  optical  waves 
which  interact  via  DFWM  in  BaTiOg.  This  unique  experimental 
technique  allows  parallel  processing  using  two  input  signals  at 
rates  which  are  not  limited  by  the  r^'Sponse  time  oE  the 
photoref ractive  crystal.  It  is  recommended  that  further 
experiments  be  conducted  in  this  area  to  increase  the 
understanding  of  this  promising  new  technology. 
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APPENDIX  A 


SECTION  1 
INTRODUCTION 

The  purpose  of  this  study  of  photoref ractive  materials  Is  to  develop 
a  theory  of  photoref ractivity,  with  emphasis  on  simple  physical 
explanations  of  the  results,  in  order  to  Interpret  experimental  results, 
suggest  experiments  to  be  performed,  identify  materials  parameters  that 
should  be  measured,  select  and  rate  types  of  materials  and  specific 
material,  and  seek  Improved  performance.  Photoref ractivity ,  which  is  the 
change  in  the  refractive  index  of  a  material  with  irradlance,  is  analogous 
to  photoconductivity,  which  is  the  change  in  the  conductivity  of  a 
material  with  irradlance.  A  related  effect,  the  photovoltaic  effect,  is 
the  generation  of  a  voltage  by  Irradlance. 

This  interim  report  was  prepared  in  order  to  make  the  mathematical 
details  of  the  treatment  of  photoref ractivity  available  at  the  earliest 
possible  time  and  to  solicit  comments  on  the  analysis.  The  details  of  the 
mathematical  treatment  of  a  simple  photorefractivlty  model  that  we  believe 
contains  the  essential  fractures  required  to  explain  photorefractivlty  in 
LlNb03,  KNb03,  BaTi03,  and  BSO  (B1]^2S1^®20)  presented. 

The  results  of  the  present  report  will  be  applied  to  these  materials,  and 
others  of  current  Interest  in  a  subsequent  report.  The  emphasis  to  date 
has  been  on  lithium  nlobate  (LlNb03)  because  there  is  more  Information 
available  on  this  material  than  on  other  materials,  some  of  which  may  be 
better  candidates  for  the  intended  applications.  This  information  Is 
useful  in  establishing  the  validity  of  the  theory,  which  will  be  done  to 
the  extent  possible  before  applying  the  theory  as  discussed  above. 

In  a  subsequent  report,  the  speed  and  sensitivity  of  the  various 
materials  will  be  estimated,  tables  of  values  of  material  parameters  will 
be  compiled,  materials  will  be  compared,  and  methods  of  improving  the 
performance  will  be  discussed  along  with  the  effects  of  changing  various 
system  parameters.  An  important  feature  of  the  subsequent  rct  ort  will  be 
a  table  establishing  the  operating  regime  for  various  materials  under  diverse 
operating  conditions.  Table  1  lists  the  various  limits  of  laser-pulse 
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duration,  electron-transport  distance,  and  electron-transport  mechanism. 
Finally,  the  present  results  will  be  related  to  previous  results  In  the 
literature. 


As  an  example  of  Improving  performance,  the  most  Important  method  for 
Increasing  the  speed  of  lithium  niobate  is  to  reduce  (the  opposite  of  oxidizing) 
the  sample,  which  Increases  the  ratio  of  the  densities  of  "filled"  (N^q,  for 
Fe^"*")  to  "empty"  (N+q,  for  Fe^)  levels.  This  ratio  Ngo/N+Q,  not 
the  Individual  values  of  NgQ,  determines  the  writing-  and  erasing-speed  (for 
fairly  wide  ranges  of  values  of  Ngg  as  seen  in  the  expression 

for  the  dielectric-relaxation  frequency  in  Equation  3.7. 

The  dielectric-relaxation  frequency  determines  the  write-  and  erase-times. 

(Rewriting  a  new  grating  during  erasure  is  known  to  affect  the  erase  time.) 
There  are  several  characteristic  times  in  the  photoref ractive  effect  including 
the  following:  the  recombination  time  of  (electrons  and  enpty  optical-donor 
levels,  the  chartacteristic  time  for  the  irradiance  to  generate  free  carriers, 
the  dielectric  relaxation  frequency,  the  trapping  time,  and  the  thermal  freeing 
time.  The  time  constant  that  controls  the  writing  and  erasing  times  is  the 
dielectric  relaxation  frequency  for  all  cases,  namely,  diffusion-,  applied- 
field-,  and  photo-voltaic  effect-controlled.  It  is  important  to  use  the  correct 
carrier  density,  which  Is  the  carrier  density  generated  by  the  write-or  erase- 
beams,  in  the  dielectric  relaxation  frequency.  It  is  assumed  for  now  that  the 
explanation  of  the  difference  between  write-  and  erase-times  and  between  erase 
times  for  different  c-axls  directions  is  that  the  writing  of  a  "secondary 
grating"  by  the  erase-  and  dif fracted-beams  is  correct.  This  effect  will  be 
Investigated  in  a  latter  report. 


Table  1.  Operating  Conditions  and  Material  Parameters 


Laser- Pulse 

Short-Time 

Steady-State 

Duration,  t 

Electron-Transport 

Short-Transport 

Long-Transport 

Distance,  A 

Ak  «1 

Ak  »1 

g 

g 

Electron-Transport 

Diffusion 

Applied  Field 

Photovoltaic 

Mechanism 

Effect 
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The  dielectric  frequency  given  by  Equation  3.7,  the  result  for  the  space 
charge  electric  field  E  that  gives  rise  to  the  grating,  and  the  relation  for  the 
grating  efficiency  for  a  given  E  determine  the  performance  of  the  photoref ractl ve 
material.  This  grating  electric  field  is  not  to  be  confused  with  the  electric 
field  of  the  optical  beams*  The  response  time,  which  is  the  read-or  wrlte-tlme 
in  the  absence  of  writing  a  secondary  grating  during  erasure  is  Independent  of 
the  grating  wave  vector  kg  and  diffusion  distance  A.  By  contrast,  the  electric 
field  E  depends  on  A  and  kg  in  general. 

The  photorefractive  effect  is  as  follows:  Illuminating  a  material  with 
light  generates  free  carriers,  which  are  assumed  to  be  conduction  electrons  for 
simplicity  in  the  present  treatement.  Holes  can  contribute  to  the  photo- 
refractive  effect  in  general,  but  are  often  negligible.  The  conduction 
electrons  are  assumed  to  come  from  such  impurity  levels  in  the  electronic  band 
gap  as  Fe^"*"  levels  in  LlNb03,  but  coul<^,  in  general,  come  from  direct 
exitation  of  electrons  across  the  gap,  by  one-  or  two-photon  absorption,  for 
example.  These  Impurity  levels  will  be  called  optical-donor  levels.  The 
electrons  drift,  if  they  experience  an  electric  field  or  diffuse  to  a  new 
location  in  the  crystal  where  they  recombine  with  empty  optical-donor  levels. 

The  space-charge  electric  field  generated  between  the  electrons  in  their  new 
positions  and  '  e  empty,  positively  charged  optical-donor  levels  which  remain  in 
their  orlglnci.  positions  causes  the  refractive  index  to  change  by  the  electro¬ 
optic  effect. 

In  addition  to  drift  in  an  externally  applied  dc  field  and  diffusion,  the 
photovoltaic  effect  can  contribute  to  the  photolnduced  electric  field  E.  The 
components  of  E  resulting  from  diffusion,  E^j,  drift  in  the  net  electric  field 
(E  +  ^app’  ^app  externally  applied  field  assumed  to  be 

spatially  constant],  Ej^.,  and  the  photovoltaic  field,  Ep, 

E  -  Ed  Ejr  +  Ep 

are  discussed  below.  A  typical  value  of  E  is  10^  V/cm.  The  resulting  change  in 
the  refractive  index,  which  is  called  the  grating  index,  is 


"g  ■  2  ”r^ 


(1.1) 


where  n^.  Is  the  refractive  Index  and  r  is  the  appropriate  linear  electro-optic 
coefficient,  having  a  typical  value  of  10“®  cm/V,  which  gives 
Qg  3  0.5  (2.3)3  10-8  10**  »  6  X  lO"** 

In  such  optical  processing  as  four-wave  mixing  used  for  phase  conjugation, 
the  Irradlance  changes  from  one  position  to  another  in  the  crystal.  For  example, 
the  beat  patterns  of  two  overlapping  beams  causes  a  periodic  irradlance.  To  be 
concrete,  we  assume  a  sinusoidal  irradlance  distribution  for  writing  a  grating 

I  »  <I>  (1  +  M(coskgz)  (1.2) 


where  <  >  denotes  the  spatial  average 


(1.3) 


and 

kg  •  4ti  n^.  sin  9/X,  (1*4) 

with  20  the  full  angle  in  the  material  between  the  direction  of  propagation  of 
the  two  beams,  X  is  the  wavelength  of  the  light  in  vacuum,  and  X/n^.  is  the 
wavelength  of  the  light  in  the  material. 

The  electric  field  and  the  resulting  index  of  refraction  ng  have  a  spatial 
frequency  component  at  Che  grating  wavelength  kg.  This  index  grating  can 
diffract  an  incoming  light  beam  with  efficiency^  (for  a  thick  grating) 
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j  /  o  -  irn  n  d 
-ad/cos9  .2  r  a 
■  e  sin  ■  *  v  , 

2^  COSO 


where  a  is  the  optical  absorption  coefficient,  d  is  the  thickness  of  the  sample, 
29  Is  the  full  angle  between  the  two  writing  beams,  n^  Is  the  Index  of 
refraction  of  the  material,  ng  is  the  change  In  the  index  resulting  from  the 
photoref ractlve  effect,  and  X  la  the  vacuum  wavelength  of  the  writing  beams. 

There  Is  general  agreement  In  the  literature  that  this  description  of  the 
photorefractive  effect  la  correct.  However,  there  Is  disagreement  as  to  such 
details  of  the  model  as  which  electron  transport  mechanisms  —  drift,  diffusion, 
or  asymretrlc  photogenetation  of  electrons  (photovoltaic  effect)  —  Is  operative 
in  a  ven  material  and  why  the  grating) efficiency  increases  or  decreases  with 
increasing  pulse  duration  tp  for  a  given  energy  density  Itp,  rather  than 
remaining  constant  as  predicted  by  early  theories. 

The  model  is  discussed  in  Section  2.  In  order  to  make  the  first  analysis  as 
simple  as  possible,  the  first  case  considered  (see  Section  3)  is  that  of  an 
applied  field  only.  Diffusion,  the  photovoltaic  effect,  depletion,  and 
saturation  are  neglected.  The  case  of  diffusion  only  and  of  photovoltaic  only 
are  consider  in  Sections  4  and  5  respectively  in  order  to  determine  the  defects 
of  these  two  transport  processes  acting  Individually.  Depletion  and  saturation 
are  considered  in  Section  6.  The  steady-state  grating  electric  field  with  an 
applied  electric  field,  diffusion,  and  the  photoelectric  field  Included  is 
derived  in  Section  7.  The  results  of  Kukhtarev  and  coworkers ^  and  of  Amodel^ 
are  related  to  the  present  work  and  discrepancies  between  the  three  sets  of 
results  are  resolved  in  Section  8.  Simple  pedantic  models  of  the  photovoltaic 
effect  and  of  dielectric  relaxation  are  given  in  Sections  9  and  10.  A  summary  of 
the  mathematical  results  are  given  in  Section  11.  The  Appendix  contains  a  list 
of  symbols  and  typical  values  of  parameters  for  lithium  nlobate. 
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Electroscaclc  units  (esu)  are  used  in  the  equations.  Numerical  results  are 
given  in  both  electrostatic  units  and  rationalized  aks  units.  Conversion  factors 
are  listed  in  the  Appendix  for  the  convenience  of  the  reader. 

Important  results  are  denoted  by  underscored  equation  numbers  and  used  in 
the  summary  (Section  11). 


SECTION  2 


MODEL  OF  PHOTOREFRACTIVITY 

The  analysis  of  even  the  simplest  model  of  photoref ractlvlty  Is  complicated 

0 

when  carried  out  without  approximations  because  the  coupled  material's  equations 
and  Maxwell's  equations  are  nonlinear,  coupled,  partial  differential  equations. 

In  order  to  develop  the  physical  understanding  required  for  materlals-liiq>rovement 
programs  and  operation-improvement  programs,  a  simple  model  Is  discussed  first, 
then  the  effects  of  the  assumptions  made  In  the  simple  model  are  discussed. 

The  simplest  model  of  photoref ractlvlty  Is  as  follows:  The  material  Is  an 
Insulator  with  Njj  optical  levels  per  cubic  centimeter  In  the  electronic 
bandgap.  See  Figure  1.  The  levels  are  sufficiently  deep  that  they  are  not 
thermally  excited  Into  the  conduction  band.  The  levels  are  called  optical-donor 
levels  because  electrons  from  these  levels  are  photoexclted  Into  the  conduction 
band.  The  density  of  Che  optical-donor  levels  Chat  are  filled  with  electrons  Is 
Ng  and  Che  density  that  are  unfilled  Is  N^.  For  example,  could  be  the 
Fe^"^  density  In  llchlm  nlobate,  Rf.  and  Fe^  density,  and  Nj)  ■  Ng  +  R|.  the 
total  Iron  density  Including  both  Fe^"^  and  Fe^'*".  In  general,  there  are  many 
filled  levels  Ng  and  many  empty  levels  Nl|..  The  ratio  Ng/N+  can  be  varied 
by  oxidizing  or  reducing  the  sample  or  by  doping  with  compensating  impurities. 

In  order  to  make  the  crystal  have  charge  neutrality.  There  nust  be  negative 
impurities  or  Imperfections  present  to  conq>ensate  for  the  positive  charges  of  the 
Iron.  The  density  of  these  levels  that  compensate  for  the  extra  positive  charge 
of  the  Fe^"*"  Is  equal  to  the  value  NL^)  of  NL|.  before  the  Irradlance  Is  turned 
on.  Subscript  zeros  denote  Che  values  of  parameters  at  t  ~  0,  l.e.,  the  time  at 
which  Che  Irrandlance  Is  turned  on. 

Holes  are  neglected  for  simplicity.  Electrons  are  generated  from  Che  filled 
optical-donor  levels  to  the  conduction  band  by  the  Irradlence  and  recombine  with 
empty  optical-donor  levels.  Trapping  Is  neglected  In  this  first  model.  We  used 
the  common  photoconductivity  terminology  that  electrons  In  traps  can  be  thermally 
re-exclted  Into  the  conduction  band.  Thus,  the  capture  of  electrons  by  the  empty 
optical-donor  levels  Is  called  recombination,  rather  than  trapping.  This 
recombination  Into  Che  optical-donor  levels  Is  not  to  be  confused  with 
recombination  of  holes  and  electrons  across  Che  gap. 


Saturation  and  depletion  of  the  optical-donor  levels  are  neglected  when 
possible.  However,  they  must  be  included  in  the  analysis  of  the  steady  state 
because  the  correct  steddy  state  is  not  attained  in  general  when  saturation  and 
depletion  are  neglected. 

The  simple  model  is  believed  to  be  sufficient  to  explain  the  performance  of 
Iron-doped  lithium  nlobate,  and  It  may  be  sufficient  for  the  other  materials  of 
current  Interest.  Extensions  of  the  model  and  the  applicability  to  other 
materials  will  be  considered  In  subsequent  studies  and  reports. 

The  number  of  electrons  per  cubic  centimeter  generated  by  the  light  Is  equal 
to  the  number  of  photons  per  cubic  centimeter  absorbed,  which  Is  equal  to  the 
absorbed  energy  per  cubic  centimeter,  al,  (where  a  Is  the  optical  absorption 
coefficient)  divided  by  the  photon  energy  hm: 

(3n/3t)j  -  (3N|/3E)j  -  al/ha» 

It  was  assumed  that  one  free  electron  is  generated  for  every  absorbed  photon. 

This  generation  rate  can  be  written  In  several  convenient  forms  as  follows: 

al/hw  -  o  N  /hu  -  y,N  -  <Y_>  N  (1  +  Mcos  k  z)  ,  (2.1) 

De  le  le  g 

where  Equation  1.2  was  used  for  I  and 

a  -  ,  (2.2) 

Yi  -  OpI/hm  ,  <Yj.>  -  Op  <I>/hu)  (2.3) 

where  Op  Is  the  optical  absorption  cross  section  of  the  optical  donor  levels. 


The  change  In  the  electron  concentration  resulting  from  movement  of  the 
electrons  away  from  the  generation  region  (and  from  the  motion  of  the  electrons 
In  general)  is  described  by  the  charge  continuity  equation: 

3p/3t  -  7  •  J  -  -dJ/dz  (2.4) 

where  the  charge  density  Is 

p  -  e(N^  -  -  n)  -  e(N^^  -  n)  (2.5) 

where  e  >  |e|,  the  absolute  magnitude  of  the  electron  charge;  - 

and  J  Is  the  electrical  current  (^e  times  the  net  rate  at  which 
electrons  cross  a  unit  area  of  surface).  Substituting  Equation  2.5  Into 
Equation  2.4  gives  the  first  central  equation  ~  Equation  2.8  below.  The 
recombination  rate  Is 

c.co.b  ” 

where  the  recombination  frequency  Is 

■  ’th“"+0  • 

uhere  s  Is  the  cross  section  for  the  capture  of  an  electron  by  an  empty  optical 
donor  level  (having  density  Hf),  and  v^)j  Is  the  average  thermal  velocity  of 
the  electrons,  which  has  the  value: 


at  room  temperature.  Combining  these  contributions  gives  the  basic  equations: 


3n/3t  -  3N^/3t  +  e"^dJ/dz 

3N^/3t  >  -3N^/3t  -  -  Y^nN^/N^ 

J  »  -eunE  +  e03n/3z  +  pN  Ic 

e 

I 

3E/3z  -  4  a“^P 

1TC 

P  -  e(N^  -  -  n)  -  e(N^^  -  n) 

E  »  E  +  E 

app  sc 


(2.8) 

(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 


where  D  is  the  diffusion  coefficient,  p  is  the  photovoltaic  coefficient,  c  is  the 
unit  vector  along  the  photovoltaic  direction,  and  e  -  |e|.  The  contributions  to 
the  current  in  Equation  2.10  are  as  follows:  -eynE  =  is  the 


usual  electrical-conductivity  current  J  *  oE;  e03n/dz  is  the  diffusion  current; 
and  pNglc  is  the  photovoltaic  current,  as  discussed  in  Section  9. 

These  system  Equations  2.8  throu^  2.13  are  nonlinear  because  they  contain 
products  Mfn  and  nE  (where  E  is  a  function  of  n  and  N^.).  Thus,  their 
solution  is  nontrivial,  and  the  complexity  of  the  general  results  would  pbscure 
simple  interpretation.  The  following  limiting  cases  will  therefore  be 
considered: 


(a)  kgA  «  1  and  kgA  »  1 

(b)  drift,  diffusion,  and  photovoltalc-effect-domlnated 

(c)  short-time  limit  and  steady  state. 


where  A  is  the  electron-transport  distance,  as  discussed  below. 

For  all  cases,  a  first  approximation  which  will  be  relaxed  below  for  special 
cases,  the  equations  are  linearized  by  neglecting  depletion  and  saturation. 
Depletion  is  the  reduction  in  the  absorption  coefficient  a  resulting  from 
depletion  of  the  electrons  in  the  optical-donor  levels,  and  saturation  is  the 
reduction  in  the  recombination  rate  Yp  resulting  from  filling  of  the 
recombination  levels  (which  are  the  empty  optical-donor  levels  N+  In  the 
present  model)  with  electrons.  Depletion  is  negligible  when  Ng  s  N^q,  and 
saturation  is  negligible  when  N+  2  N+q*  To  be  more  specific,  define  the 
changes  and  Ng^  as 


and 


N+A  ^  -  «40 


(2.14) 


eA 


2  N  - 
e 


N 


eO 


(2.15) 
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Since  3Ng/3t  *  -3N^/3t  according  to  Equation  2.9 


N 


+A 


-N 


eA 


( 


Thus,  the  approximations  of  negligible  depletion  and  saturation  are, 
respectively, 


«  \o 


( 


For  negligible  saturation,  the  optical  absorption  coefficient  o  *  op! 
is  constant.  Thus, 


a  *  ®D^eO  '  negligible  depletion.  ( 

For  negligible  saturation,  the  recombination  frequency  Vj-j^sN.^  is  constant 
thus 


for  negligible  saturation. 


(? 


SECTION  3 


APPLIED-FIELD  DOMINATED 


In  this  section,  diffusion,  the  photovoltaic  effect,  depletion,  and 
saturation  are  neglected.  The  electron-drift  distance  Is  ' 

Ag  -  WT^E  , 

where  la  the  recombination  time  and  uE  Is  the  drift  velocity, 

two  limits  of  short  drift  distance,  «  1  and  long  drift  distance 

kgAg  »  1,  will  be  considered. 

The  steady-state  solution  of  Equation  2.9  with  Ng  5  N^q  and  N+  = 
N+o  is  ' 


n  -  (1  +  M  cos  k^z) 

where 


<n>  -  «Yi>/Y^)N^q  . 


typical  value 

T  -  1/v^.  sNj_  -  0.2  ns/c,(ppm)  ■  20  psec  (3.3) 

r  th  tu  + 


for  v^h  “  cm/sec  (at  room  teiq>erature) ,  electron-ion  cross  section  for  a 
changed  recombination  center  s  ■  cm^,  and  recombination  center  density 

N+o  “  5  X  10^^  cm  ^  for  a  lightly  doped  sample  with  concentration  c+  ■ 
N+o/^atomic  sites  "  PP™*  heavily  doped  samples,  Tj.  Is  even 

shorter  than  the  value  of  20  psec  In  Equation  3.3.  Thus,  for  pulse  durations  and 
doping  concentrations  of  current  Interest,  «  t  Is  satisfied  throu^out  most 
of  the  pulse  duration.  Furthermore,  n  Is  given  by  the  steady-state  value  In 
Equation  3.1,  except  at  very  short  times  which  will  be  neglected  here. 

Physically,  for  negligible  electron  drift  distance,  the  electrons  are  trapped 
where  they  are  generated.  Thus,  setting  J  •  0  In  Equatl'^n  2.10  as  a  zeroth-order 
approximation  and  solving  for  n  gives  an  exponential  r  ^  lo  the  value  In 
Equation  3.1  with  rise  time  t^. 

Having  calculated  n  under  the  zeroth-order  appro:  itlon,  J  ■  0,  this  value 
of  n  and  the  nonzero  value  of  J  ■  -eunE  In  the  next-or  e  approximation  cause 
to  continue  to  change  after  n  has  reached  Its  steady  state.  For  short 
tlmes(l.e.,  before  the  space  charge  field  Eg^  builds  up  to  substantially 
change  the  electric  field  In  Equation  2.9  from  Its  original  value  Eq  * 

Eapp)>  Equation  2.10  gives: 


J  2  Jq- 


2  „ 
ltm  E 


app 


Substituting  this  J  Into  Equation  2.4  and  Integrating  over  time,  with  p  ■  pq  • 
0  at  t  ■  0,  gives 

p  -  -eu<n>E  k  M  t  sin  k  z  .  (3.4) 

app  g  g 

The  space-charge  electric  field  resulting  from  the  charge  distribution  In 
Equation  3.4,  obtained  by  Integrating  Poissons  equation 


is 


E  »  Yj.  t  E  M  cos  k  z 
sc  dl  app  g 


where  the  dielectric  relaxation  frequency  is 


(3.6) 


Y,.  »  4  TTo/e  “  4  ireM<n>/e  -  4  iiepa_<I>N  -/ehuv  .  sN,.  .  (3.7) 

dl  D  eU  th  +0  - 

One  of  the  important  results  of  the  present  study  is  that  the  response  time 
(the  read-  and  wrlte-tlme  in  the  absende  of  writing  of  a  secondary  grating)  is 
the  dielectric  relaxation  time  =  1/Y<jj[.  Consideration  of  the  expression 
for  Y^i  in  Equation  3.7  can  lead  to  estimates  of  the  response  time  and  to 
suggestions  on  how  to  Increase  the  speed,  as  discussed  briefly  in  Section  1. 

This  electric  field  is  in  spatial  phase  (E  ~  cos  kgZ)  with  the  irradiance,  in 
agreement  with  the  results  obtained  below  by  physical  arguments. 

It  is  instructive  to  rederive  this  electric-field  distribution  in 
Equation  3.5  by  physical  argument  in  order  to  gain  intuition.  Once  the  steady 
state  value  of  n  is  attained,  electrons  that  are  generated  at  z  by  the  irradiance 
are  trapped  at  z  +  on  the  average.  The  trapped-change  density 
(see  Equation  2.12)  can  therefore  be  written  as  the  sum  of  two  terms,  a  positive 
terra  from  the  at  z  where  an  electron  is  generated  and  a  negative  term  of 
equal  absolute  magnitude  but  displaced  by  distance  Aj?,  corresponding  to  the 
recombination  of  the  electron  at  z  +  A^.  In  other  words,  the 


current  causes  a  buildup  of  Che  displaced  crapped  electrons  and  a  corresponding 
build  up  of  the  positive  charges  left  behind  (see  Figure  2).  Integrating  the 
generatioij  equation  ON^/9t)gen  *  C9N+^/3t)gen  -  YiNeO 
for  short  times  gives 

The  corresponding  negative,  displaced,  recombination  term  is 

<"“«>rc.b  ■  “g'"  -  ''dt’i  • 

Adding  these  two  terms  and  using  Che  Taylor  expansion  cos  hg(z  -  Ag)  = 
cos  kgZ  +  ^gkg  sin  kgZ  gives  * 

eN^^  -  -e<Yj>  ’^g*  •  (3*10) 

With  p  =  and  <n>  =  (<Yi>/Yr)NgO,  Equation  3.10  gives 

Equation  3. A.  The  rest  of  the  derivation  of  the  space-charge  field  resulting  in 
Equation  3.6  is  Che  same  as  above.  The  charge  resulting  from  n  and  the 
corresponding  is  easily  shown  to  be  given  by  the  p  in  Equation  3.4  times 
Tj,/t,  which  is  negligible  for  t  »  t^,. 

As  Che  time  increases,  the  space  charge  field  Eg^  continues  to  increase 
according  to  Equation  3.6  until  it  is  no  longer  negligible.  Then  the 
approximation  E  2  E^pp  is  no  longer  valid,  and  Equation  3.5  is  no  longer 
valid.  The  electric  field  for  intermediate  times  is  more  difficult  to  calculate, 
but  the  steady-state  value  is  again  easy  to  calculate.  In  the  steady-state, 

3n/3t  *  0  and  3N+/3t  *  0;  thus.  Equation  2.8  gives  3J/3z  -  0,  or 
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Figure  2.  Schematic  illustration  of  the  positive  and 

negative  charge  distributions  for  the  case  of 
short  diffusion  distance  (A  k  »1)  showing 
the  total  charge  distribution^p  that  is  spatially 
phase  shifted  with  respect  to  the  irradlance. 


J  »  constant  =  Jq 


For  the  present  case  of  drift  only,  setting  D  ■  p  ■  0  in  Equation  2.10  gives 


J  ■  -ewnE 


(3.11) 


From  Equations  3.10,  3.2,  and  J  «  Jq, 


E  -  - 


J^/ep  <n> 

1  -i-  M  cos  Ic  z 

g 


(3.12) 


The  divergence  in  E  for  M  *  1  at  cos  kgZ  “  -1  is  unphysical.  It  arises  from 
the  assumption  that  the  transport  (drift)  length  is  zero.  For  nonzero  transport 
length,  M  is  replaced  by  a  number  less  than  unity,  as  suggested  by  Amodei^  and 
shown  by  Rukhtarev  and  coworkers. ^  The  field  E  is  not  sinusoidal  except  in  the 
limit  of  M  «  1  (where  the  expansion  (1  +  M  coskgS)"^  2  1  +  M  cos  kgZ  is 
valid). 

An  Important  feature  of  the  field  E  (which  gives  rise  to  the  grating)  is 
obvious  when  E  is  expressed  in  terms  of  Jq,  as  in  Equation  3.12,  althou^  it 
has  not  been  discussed  previously  to  knowledge.  For  open  circuit  conditions, 
the  steady-state  field  E  vanishes  because  Jq  ■  0  (no  current  for  an  open 
circuit).  In  order  to  write  a  steady-state  grating  (l.e.,  have  E  *  0)  in  the 
steady-state  with  drift  alone,  there  must  be  an  external  circuit  through  which  a 
current  can  flow.  Otherwise,  the  faces  of  the  crystal  are  charged  by  the  nonzero 
current  before  the  steady-state  is  reached,  and  the  current  approaches  zero  as 
the  time  approaches  infinity  and  the  charges  on  the  faces  approach  the  value 
required  to  make  Eg  =■  -  Egpp  (l.e.,  E  -  0).  At  intermediate  times,  the 
grating  is  nonzero. 


A  rough  approximation  to  the  full  time  dependence  for  this  case  (of 
kgAdr  <<  1.  ^app  don^naCed,  negligible  depletion  and  saturation,  and 
J  “  0)  is  obtained  as  follows:  The  charging  of  the  faces  of  the  crystal  gives, 
from  Appendix  A; 


E  -  E  e 
app 


-^dl*^ 


(3.13) 


where  the  dielectric  relaxation  frequency  Is  given  by  Equation  3.7.  Formally 
replacing  E  =  ^app  Equation  3.6  by  E  in  Equation  3.13  and  Including  the 
spatially  constant  space-charge  field  gives 


-Y..t 

E  =  e  (E  +  E  M  cos  k  z) 

app  dl  app  g 


(3.14) 


The  maximum  of  the  cos  k_2  term  occurs  at 


^  ’  "'dl  '  ’■'^'^dl 


(3.15) 


and  has  the  value  (writing  exp  (-1)  to  avoid  confusion  with  the  electron  charge): 


E  =»  exp  (-1)  E  M  cos  k  z 
mx  r  V  /  gpp  g 


(3.16) 


Since  Che  dielectric  relaxation  time  is  Che  time  at  which  the  steady-state 
is  approached,  as  an  approximation  Co  Che  full  time  dependence,  we  multiply  the 
short-time  electric  field  by  exp  (“Ydit)  and  add  the 

result  to  f^  =  [l  -  exp(-Yjj^t)  cimes  Che  steady-state  value  Egg 


B  2  E 


-"di*^ 


short-t 


+  E  (1  -  e 
88 


■^di^ 


(3.17) 


The  square  of  f  was  used  rather  chan  f  Co  insure  Chat  E  = 

at  short  times.  This  approximation  in  Equation  3.17  does  not  correspond  to  an 
exact  solution  of  Che  systems  equation.  It  Is  simply  an  Intuitive  approximation 
to  the  full  time  dependence.  It  does  give  the  correct  short-time  field  and 
steady-state  field,  and  it  does  develop  with  time  constant  ~  I/Ydi. 

The  value  of  E  at  Intermediate  times  may  not  be  highly  accurate. 

With  Eg^Qp^.j.  given  by  Equation  3.6  and  Egg  given  by  Equation 
3.12,  substituting  In  Equation  3.17  gives 


(3.18) 


Next  consider  the  case  of  »  1,  where  » 

PTrEgpp.  Since  the  electrons  then  drift  a  distance  long  with  respect  to  the 

grating  distance  kg“^,  the  electron  density  is  approximately  constant  spatially 

for  t  >  T  .  Since  the  f ree-electron  generation  rate  (3n/3t)  ••  YtN  and 

recombination  rate  (3n/3t)  v  ”  Y  n  are  the  same  as  for  the  case  of  It  A, 

rcmb  r  g  dr 

<<  1,  the  value  of  this  constant  n  is  just  the  average  value  of  n  from 

Equation  3.1 


I 


n  2  <n>  -  «Yi>/Yj.)N^O  .  (3.19) 

The  recombination  then  occurs  evenly  over  all  z  (rather  than  at  a  short  distance 
from  the  generation  region  as  was  the  case  for  kgA^j.  <<  1).  Figure  2  (for 
kgAdr  <<  1)  is  then  replaced  by  Figure  3  for  the  case  of 

kgAdrO  '>'>  ^•)»  *”<1  the  second  term  in  Equation  3.9  Is  replaced  by  zero,  cor¬ 
responding  to  spreading  the  charge  cos  kg(z  -  Aj^)  evenly 

throughout  the  crystal.  The  resulting  value  of  eN^./^  *  (®%A)gen 
(eN^.^)^gg^^  from  Equation  3.8  and  the  modified  Equation  3.9  Is 

eNw.  -  e<T_>  N  .Mt  cos  k  z  .  (3.20) 

'TV  L  eu  g 

With  p  =  eN+Q  and  Equation  3.5,  this  gives 

E  “  4ite<YT>  N  -(k  e)tM  sin  k  z  (3.21) 

sc  I  eO'  g  '  g 
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Figure  3.  Schematic  illustration  of  the  positive  and  negative 
charge  distributions  for  the  case  of  long  diffusion 
distance  (A  kg  >>1)  showing  the  total  charge 
distribution  p  that  is  spatially  in  phase  with  the 
irradiance . 


which  Is  written  as 


E  Y.-tE.Msln  k  z  v 
sc  dl  A  g 


(3.22) 


where 


E*  =  T_/Mk  (3.23) 

A  r  g  - 

Is  the  field  at  which  the  drift  velocity  (lE/^  Is  equal  to  kg“VTp. 

The  steady-state  value  of  E  for  arbitrary  Jq  Is,  from  Equation  3.11  and  J 

=■  Jq. 

I 

E  -  -J^/eua  .  (3.25) 

The  corresponding  steady-state  value  of  the  drift  distance  is 

“  V*V  • 

For  kgAjpgg  It  n  Is  approximately  constant  as  discussed  above.  Then 
Equation  3.25  gives  E  >  constant;  thus,  there  Is  no  grating  (which  requires  a  z 
dependence  of  E)  In  this  order  of  approximation.  Keeping  higher  order  terms  In 
gives  a  nonzero  grating. 

For  kgA^pgg  «  1,  n  Is  given  by  Equation  3.1  and  there  Is  a 
nonzero  grating,  which  Is  given  by  Equation  3.26.  Thus,  Equation  3.25  gives  the 
steady-state  field 


E  » 


jQ/ew<[i> 


(3.27) 


1  +  a  M  cos  k  z  ’ 
g 

where 

a  -  1  for  JQkg/eYj.  «  1 

a  *  0  for  J-k  /ey  »  1 

0  g  r 

Equations  3.17,  3.22,  and  3.27  give 

"Y  t  * 

E  3  Y..te  E.M  sin  k  z 
dl  A  g 

-  [(jQ/ey<n»/(l  +  a  M  cos  k  z)](l  -  e  ^  )  .  (3.28) 

O 
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SECTION  4 


DIFFUSION  DOMINATED 


As  In  the  previous  case  of  drift  only,  the  charge  density  for  the  case  of 
kgA^  «  1  (where  Is  the  diffusion  distance  discussed  below)  for 
diffusion  only  can  be  obtained  from  physical  arguments  similar  to  those  used 
above  in  the  discussion  of  Figure  3.1.  For  diffusion,  half  of  the  electrons 
diffuse  to  the  left  and  half  diffuse  to  the  right,  in  contrast  to  the  drift  case 
above  in  which  all  electron  drift  to  the  right  in  Figure  2.  Thus,  Equation  3.9 
is  replaced  by 


(4.1) 


(eN  ) 

^  +A''rcmb 


“e<Y  >N  -t[l  +  cos  k  (z  - 
I  eO  '■  2  g 

+  "T  M  cos  k  (z  +  A  )  ] 


Using  the  Taylor  expansion 

cos  k  (z  ±  A.)  3  cos  k  z  ±  k  A,  sin  k  z  -  "T 
g  d  g  g  d  g  2 

and  adding  Equations  4.1  and  3.8  gives 

eN  -  e<Y_>  N  „k  ^A  ^t  cos  k  z 

+A  2  I  eO  g  d  g 


1,  2.  2 

k  A  cos  k  z 
8  d  g 


(4.2) 
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where  D  is  the  diffusion  coefficient.  Using  n  given  by  Equations  3.1  and  3.2, 
Equation  A. 2  and  p  =  eN+^  give 


P 


2 

<n>Mt  coe  k  z 
g  g 


The  corresponding  electric  field  is,  from  Equation  3.5 


E  -  E 

sc 


2  ^di*^ 


I 

E.M  sin  k  z 
d  g 


(4.4) 


(^•5) 


where 


E  5  Dk  /u 
g 


(A. 6) 


is  the  electric  field  for  which  the  drift  velocity  uE^j  is  equal  to  the  grating 
distance  divided  by  the  diffusion  time  (Dkg^)~^. 

The  steady-state  value  of  E  for  the  present  case  of  kgAj  <<  1  and 
diffusion  only  is  also  easy  to  obtain.  For  diffusion  only,  the  current,  which 
must  be  spatially  constant  in  the  steady  state,  is  zero,  J  “  Jq  =  0,  because 
there  is  no  applied  field  and  no  surface  charges.  Setting  J  -  0  and  p  ■  0  in 
Equation  2.10  gives 

•-‘i 

IJ 

•"'v 

*.  i 

’.v’ 

b 

A- 2  8 

•: 


unE  “  D3n/3z 


(^.7) 


Since  J  »  0,  n  Is  given  by  Equations  3.1  and  3.2.  With  Equation  4.7  this  gives 


E  .Main  k  z 
_d _ g _ 

1  +  M  cos  k  z 
g 


(4.8) 


For  M  <<  1,  E  5  Ej  M  sin  kgZ  Is  in  spatial  quadrature  with  the  Irradiance. 

When  M  is  not  small,  E  is  not  sinusoidal,  just  as  it  was  not  for  the  case  of 
drift.  By  symmetry,  the  Fourier  component  in  spatial  phase  is  zero. 

The  electron  transport  stops  when  the  drift  in  the  space-charge  field 
increases  to  the  value  at  which  this  drift  cancels  the  transport  from  diffusion. 
Thus,  the  characteristic  time  for  the  steady  state  to  be  approached  is  just  the 
characteristic  time  for  charge  buildup,  which  is  the  dielectric-relaxation  time 
”  ^/^dl*  Tdi  is  defined  in  Equation  3.7.  This  time 

constant  can  be  verified  as  follows:  For  convenience,  set  sin  kgZ  -  1  and  cos 
kgZ  «  0  in  Equations  4.5  and  4.8.  Equating  the  short-time  solution  in  Equation 
4.5  to  the  steady-state  value  in  Equation  4.8  and  solving  for  t  gives  one 
definition  of  the  characteristic  time  (that  would  equal  the  exponential  time 
constant  if  the  E  field  increased  exponentially).  This  gives  t  =  which 

is  approximately  equal  to  the  dielectric-relaxation  time  as  expected. 
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From  Equations  3.17,  4.5,  and  4.8, 


^  ~=  ^  i  ^d  ^ 

-Y,,t  , 

+  [E  M  sin  k  z/(l  +  M  cos  k.  z)  ]  (l  -  e  ^  ) 

8  8 

Next  consider  the  case  of  kgA^j  »  1.  In  this  long-diffusion-distance 
limit  for  short  times,  the  electron  density  is  constant  and  is  given  by 
Equation  3.17,  just  as  for  the  case  of  kgA^j^.  »  1  considered  in  Section  3. 

The  electric  field  at  short  times  is  therefore  given  by  Equation  3.20. 

sc  di  rg  g  ’  -i - - 

where 


(4.9) 


as  in  Equation  3.21. 


E 

rg 


Yp/uk 


g 
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The  steady-state  value  of  E  from  Equation  4.7  with  n  =  const  Is  E  =  0.  A 
nonzero  value  of  E  in  the  steady-state  results  from  the  next  order  correction  to 
n.  An  example  of  the  calculation  of  the  steady-state  value  of  E  is  given  In 
Section  8  in  the  discussion  of  the  results  of  Kukhtarev  and  coworkers. 

From  Equations  3.17  and  4.10 


E  5 


-^di'^ 


E  M 
rg 


sin  k  z 
g 


+  E 


Dss 


(1  - 


-W.; 

e  J 


(4.11) 


from  Section  8  for  the  case  of  N+q  «  NgO, 


As  an  example  of  Epgg, 
from  Equation  8.10 


-  (4neN^i/ek  )  M  sin  k  z 


(4.12) 


SECTION  5 


PHOTOVOLTAIC-EFFECT  DOMINATED 

The  photovoltaic  effect  Is  discussed  In  Section  9.  It  Is  easy  to  show  that  the 
results  In  Section  3  (for  applied  field  only)  are  valid  for  the  photoelectric  effect 
if  Egpp  Is  replaced  by  the  steady-state  photovoltaic  field: 

E  “Ay/v**pN  <I>/eu<n>  (5.1) 

p  p  r  e 


from  Equation  9.8.  This  prescription  of  replacing  E^pp  by  Ep  can  be 
understood  as  follows:  The  central  feature  of  the  Glass ^  explanation  of  the 
photovoltaic  effect  Is  that  the  electrons  generated  by  the  Irradiance  are  emitted 
preferentially  In  a  direction  given  by  ^  unit  vector  c.  The  electrons  are  displaced 
through  an  effective  distance  Ap  in  the  ^  direction,  as  discussed  in  Section  9. 

Thus,  Equation  3.9  is  valid  If  A^j^  is  replaced  by  Ap.  The  prescription  of 
replacing  E^pp  by  Ep  In  the  expression  A^j.  <■  uE^pp/Yj.  Is 

equivalent  to  replacing  A^p  by  Ap,  as  seen  by  making  the  replacement  (5.1)  in 
the  Ajp  expression.  Replacing  Ejjpp  by  Ep  in  Equations  3.4  and  3.5  and 
Integrating  Equation  3.5  for  short  times  gives,  for  Apkg  «1, 

E  5  Yj.t  E  M  cos  k  z.  (5.2) 

di  p  g  - 


This  result  can  also  be  obtained  from  the  continuity  equation  3p/3t  ■  3J/3z,  with 
I  »  -eynE  +  pNgl  =  pNgl  for  short  times.  Notice  that  kgAp  «  1  is  always 
satisfied  for  cases  of  Interest  (because  Ap  2  0.1  nm  and  kg  "  1/500  nm). 


As  the  time  increases  further,  the  space  charge  builds  up  until  J  «  Jq  - 
constant,  which  is  the  condition  for  the  steady  state  as  discussed  above.  Setting 
J  =»  Jq  in  Equation  2.10  (with  D  «  0),  solving  for  the  steady  state  field,  and 
using 


pN  <I>/ep<n>  “  Ay  /u  ■  E 
eo  p  r  p 


gives 


E  -  [e  (1  +  M  cos  k  z)  -  Jrt/eu<n>]/(l  +  M  cos  k  z).  (5.3) 

P  g  u  g 

I 

The  first  term  Ep  is  independent  of  z.  Thus,  there  is  no  steady-state 
grating  for  open  circuit  conditions  (Jq  =  0),  just  as  for  an  applied  field. 
However,  there  is  a  nonzero  correction  to  this  result  when  diffusion  and/or  an 
applied  field  are  added.  The  <n>  in  the  numerator  in  Equation  5.3,  which  comes 
from  n  s  <n>(l  +  M  cos  kgz)  is  a  zeroth  order  approximation.  Kukhtarev  and 
coworkers  obtained  the  following  first  order  correction  for  the  case  of  kgA  «  1: 

n  5  <n>  [l  +  M  (1  -  5)  cos  k  z]  (5.4) 

g 
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where 


+  E  )/E 
«PP  g 


(5.5) 


with 


E  -  4ireN^/ek^ 


(5.6) 


Expanding  In  powers  of  C  gives 


1  +  M  cos  k  z 

_ g _ 

1  +  M  (1  -  5)  cos  k  a 

g 


C  M  cos  k  z 
s  1  + - E- 


1  +  M  cos  k  z 
g 


Thus,  the  z-dependent  terins  In  Equation  5.3  for  the  case  of  kgA  «  1  are 


E  s  (E  C  M  cos  k  z  -  J„/ev<a>)/(l  +  H  cos  k  z).  (5.7) 

p  g  u  g  - 


[My  preliminary  results  give  a  different  correction  to  n  from  the  value 
in  Equation  5.4  by  Kukhtarev  and  coworkers.  The  discrepancy,  which  is 
unimportant  in  the  present  cases,  in  which  the  photovoltaic  effect  is 
negligible,  will  be  resolved  in  a  later  note. ] 


Combining  Equations  3.17,  5.2,  and  5.7  gives,  for  Apkg  «  1, 


-Y,.t 

E  5  Tj.t  e  EM  cos  k  z 
dl  p  g 

E  CM  cos  k  z  -  J-/eu<n>  “Y,.t  2 

*  "(W  »  CO.  k  -  n  -  .  "S  .  for  »  1 

s 


(5.8) 
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For  Egpp  and  D  included  and  KgA  »  1,  n  is  approximately  constant. 

Then,  replacing  1  +  M  cos  kgZ  by  1  in  the  denominator  in  Equation  5.3  gives, 
for  kgA  »  1 


E  2  E 

P 


(1  +  M  cos  k  z)  -  J/ew<n> 

g 


(5.9) 


Equations  3.17,  5.2,  and  5.9  give  for  this  case  of  kgA  »  1, 


— Y  t  “Y  t  2 

E  2  Y..te  E  M  cos  k  z  +  E  Mfcos  k  z)  fl  -  e  )  .  (5.10) 

dl  p  g  P  g  - - 


As  an  interesting  side  point,  for  an  applied  field  the  voltage  at  small 
times  is  E^pp,  and  the  voltage  decays  to  zero  as  the  space-charge  field 
builds  up  as  time  approaches  infinity.  The  space-charge  field  builds  up  to 
cancel  the  applied-field,  given  a  net  zero  field.  For  the  photovoltaic  case, 
the  electric  field  is  small  at  small  times  and  approaches  Ep  as  the  space- 
charge  field  builds  up  as  time  approaches  Inflnflty.  (The  field  Ep  is  the 
space-charge  field  generated  by  the  photovoltaic  effect).  The  space-charge  field 
builds  up  and  cancels  the  photovoltaic  current,  giving  a  nonzero  field. 
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SECTION  6 


DEPLETION  AND  SATURATION 

The  conditions  for  no  depletion  of  the  electrons  In  the  Ng  levels  (which 
determine  the  optical  absorption  coefficient  S  ■  OpNg)  and  no  saturation  of 
the  available  recombination  levels  N^  (which  determine  the  effective 
recombination  frequency  Vj.|,sN+)  are  given  in  Equations  2.17  and  2.18, 
respectively.  The  first  condition  required  for  no  depletion  and  no  saturation  is 
that  the  steady-state  value  (<Yio^/Tr)NgO  of  the  aiiq>lltude  of  n,  from 
Equations  3.1  and  3.2,  oust  be  sufficiently  small  for  Equations  2.17  and  2.18  to 
be  valid.  By  conservation  of  electrons,  the  number  of  electrons  NgQ 
available  to  be  photoexclted  at  t  *  0  must  always  be  equal  to  the  total  number  of 
electrons 

I 

N  .  -  <N  >  +  <n>  (6.1) 

eu  e 

where  <  >  is  spatial  average,  as  defined  in  Equation  1.3. 

Taking  the  spatial  averages  of  Equations  2.17  and  2.18  gives 

|<N^^>|  «  N^  ,  for  negligible  saturation  (6.2) 

and 

|<N^^>|-|<N^^>|<<  N^q,  for  negligible  depletion.  (6.3) 

Prom  Equation  6.1  and  <Ng>  -  NgO  ■  <Ng  -  NgO>  ■  <Ng^>, 

<n>  -  <N  .>  . 

eA 


(6.4) 


Thus,  Equations  3.2  and  6.4  give 


# 

which  requires 

<Yj>  «  Yj.  ,  for  negligible  depletion  (6.6) 

as  seen  by  cooqiarlng  Equations  6.3  and  6.5. 

The  corresponding  condition  for  negligible  saturation  Is  obtained  by  writing 
Equation  6.5  as  ^ 

Comparison  of  Equations  6.7  and  6.2,  with  1n+^|  -  [Ne^l  from  Equation 
2.16,  gives 

<ri>  «  .  fox  negligible  saturation.  (6.8) 

The  value  of  <I>  corresponding  to  Equation  6.6  Is  <!>  «  Igat*  where 

I  *  Y  hoj/a 
sat  r 

-  2  X  102  8-1  (4.0  X  10-19j)/5  x  lO"!®  cm^  (6.9) 

■  1.6  X  10 1 1  W/cm^  . 
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This  very  large  typical  value  of  Isat  “  10^^  W/cm^,  compared  with  typical 
cw  values  of  10  W/cm^  and  rather  large  pulsed  values  of  10®  W/cm^  (1  J/cm^  and  10 
msec)  shows  that  saturation  and  depletion  are  not  likely  to  be  problems  In  fairly 
heavily  doped  samples  even  under  short-pulse-duratlon  operation.  It  Is  possible 
that  Ngo/N+0  Equation  6.8  could  be  small,  thus  lowering  the  saturation 
threshold.  However,  In  most  current  high-speed,  real-time  applications  it  is 
desirable  to  make  small  and  N^q  large,  so  that  saturation  Is 
negligible  (Equation  6.8),  If  depletion  Is  negligible  (Equation  6.6). 

In  addition  to  these  restrictions  (Equations  6.6  and  6.8)  on  the  maximum 
irradlance  I,  there  are  restrictions  imposed  on  the  values  of  N's  In  order  that 
the  values  of  p  =  at  longer  times  (t  »  t^)  can  be  obtained.  For 

example,  the  maximum  value  of  E  In  Equation  3.16  for  the  case  of  Jq  0  and 
applied  field  only  is,  to  within  a  numerical  factor  of  order  1,  E^pp  M,  and 
the  corresponding  value  of  In^aI  fifo™  Equation  3.10  with  t  “  Ydi~ ^  aud 
^dr  “  =  HTfEappM  is 


N  5  ek  E  M^/Aire 
sat  g  app 


With  Equations  6.2  and  6.3,  this  gives 


eO 


(6.10) 


(6.11) 


where  Is  defined  in  Equation  6.10. 

According  to  Equations  6.11  and  6.10,  both  and  NqQ  must  the 
greater  than 


N 


sat 


„  6(3  X  10^  cm~^’>fl0^V/cm)fl/300)fSV/V) 
4it  4.  X  10  (ergcm) 


-3 

10  cm 


(6.12) 
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for  Egpp  -  10  kV/cm,  e  «  (2.3)^,  sin  0  *  0.5,  X  *  500  nm,  and  M  5  1.  This 
Is  a  rather  small  value  of  Ng^j.  —  corresponding  to  ^  10^®/5  x  10^^  ■  2  x 
10“^  0.2  part  per  million  of  Fe^"*"  and  Fe^.  Large  dopings  are  commonlj  used. 
Thus,  the  usual  doping  are  sufficiently  large  to  avoid  the  long-time  satulration 
and  depletion.  However,  for  small  values  of  M,  the  factor  in  Equation  6.10 
can  make  Ngg^  small. 

In  the  general  case,  Ngg^  must  satisfy  Equation  6.10,  with 

(6.13) 

where  Egn,,  is  the  maximum  value  of  the  z  dependent  electric  field  that 
gives  the  grating  amplitude. 


N  -  ek  E  M/4xe  , 
sat  g  gmx 


t 


SECTION  7 


STEADY-STATE  RESULTS 


The  steady  state  has  been  considered  for  Individual  processes  In  the 
previous  sections.  Since  the  value  of  kgA  for  all  processes  affects  the 
Individual  processes,  It  is  convenient  to  consider  the  steady  state  with  all 
processes  included.  The  key  equations  for  determining  the  steady  state  field  E 
are 


E  -  (eD3n/3t  +  pN^I  -  JQ)/epn 

E  »  (eD3n/3z  +  pN  <I>M  cos  k  z  +  ew<n><E>)/eyn 
e  g 


open  circuit 
short  circuit 


(7.1a) 
(7.1b) 
(7.1c) 
(7. Id) 


E  *  (eD3n/3t  +  pN  I)/eim  , 
e 

I 

E  ■  (eD3n/3z  +  pN  <I>M  cos  k  z)/ejin  , 
e  g 

E  *  (eD3n/3t  +  pN  <I>M  cos  k  z  +  eii<n>E  )/epn,  for  <E>  =■  E  (7.1e) 
e  g  app  app 

3E/3z  •  4iree  ^(N^  -  n) 


(7.2) 

(7.3) 


and 


-  n>  -  0  .  (7.4) 

Equation  7.1  is  obtained  from  the  steady-state  condition  J  *  Jq  >  constant 
*  f  (z,t).  The  average  field  <E>  from  Equation  7.1a  is 

<E>  -  (pN  <I>  -  J_)/eu<n>  .  (7.5) 

e  u 

Solving  Equation  7.5  for  Jg  and  substituting  the  result  into  Equation  7.1a 
gives  Equation  7.1b.  For  an  open  circuit,  setting  Jg  “  0  in  Equation  7.1a 
gives  Equation  7.1c.  For  a  short  circuit  <E>  ••  0,  Equation  7.1b  reduces  to 
Equation  7.  Id.  For  a  voltage  V^pp  applied  with  ohmic 

contacts,  <E>  ■  EJ^ppk,  and  Equation  7.1b  reduces  to  Equation  7.1e. 


Equation  7.2  is  Equation  2.11  with  a  p  given  by  Equation  2.12.  Equation  7.3  is 
obtained  from  Equation  2.9  with  3N^/3t  *  0. 

Equation  7.4  is  the  quasi  charge-neutrality  condition,  which  arises  as 
follows:  In  general,  charge  neutrality  means  that  the  charge  density  p  is  zero. 
Then,  integrating  Poisson's  equation  7  •  eE  ■  Airp  gives  E  <■  constant  for  ' 
the  current  case  In  which  e=  constant.  (It  is  easy  to  show  that  the  spatial 
dependence  of  e  -  ng  corresponding  to  the  spatial  dependence  of  grating  Index 
ng  is  negligible  In  the  cases  considered  here).  Quasi-charge  neutrality  in  the 
context  of  index  gratings  means  that  the  constant  term  <p>  in  p  *  <p>  +  pA  must 
equal  zero  in  order  to  avoid  the  unphyslcal  term  E  ~  <p>z  in  the  electric  field. 
Thus,  Equation  2.12  and  <p>  •«  0  give  Equation  7.4. 

The  case  of  short  steady-state  electron  transport  distance  is  simple.  [The 
condition  for  this  short  transport  distance  will  be  considered  later. 

Kukhtarev's  condition  5  «  1  and  Sg  «  1  are  sufficient  if  <n>  <  NL^q,  but 
are  not  necessary  in  general.  Physically  I  would  expect  kgA^g  <<  1  with 
Agg  =  I'^sa'^r  The  electron  concentration  from 

Equation  3.1  and  3.2  la  n  ••  <n>(l  -  M  cos  kgz).  Substitution  Into 
Equation  7.1a  gives 

E  *  (E  CM  cos  k  z  +  E  .  M  sin  k  t  -  J„/ey<n>)/(l  +  M  cos  k  z)  (7.6) 

n  p  g  d  g  0  g 


where 


E  hAy/m“pN  <I>/eu<n>  ,  E.-Dk/u  .  (7.7) 

p  p  r  e  ’  d  g 

+  M  cos  kgZ)  arises  from  the  deviation  of  n  from 
<n>  (1  +  M  cos  kgZ),  as  discussed  In  the  derivation  of  Equation  5.7.  [The 
consistency  of  E^  In  Equation  7.6  with  Equations  7.2  to  7.4  will  be 
demonstrated  below.  This  consistency  establishes  the  conditions  required  for 
n  =  <n>  (1  +  M  cos  kgZ)  to  be  valid]. 

The  case  of  long  steady-state  electron-transport  distance  requires  some 
care.  For  example,  formally  setting  n  *  <n>  In  Equation  7.1  would  give 


The  term  EpCM  cos  kgz/(l 


E  ■  E  (1  +  M  cos  k  z) 


for  n  »  <n> 


(7.8) 


r 


This  result  Is  Incorrect  because  the  deviation  In  n  from  <n>  gives  terms  in 
Equation  7.1  that  are  much  greater  than  Ep  in  general,  as  will  now  be  shown. 
Setting 


n  ■  <n>  +  n.  ,  N»<N>  +  N, 

0  e  e  eo 


in  Equation  7.3  and  using  gives 


+A 


-  n 


\i  +  "s/<'»A  V">/ 


M  cos  y)(l  +  »^5/<N^» 


-  N 


+0 


-  <n>  -  n, 


with  n^  «  <n>,  which  will  be  verified  ^below,  and 

N  .  -  r<N  >  M  cos  k  z  .  (7.10) 

e<S  e  g 


and 


(1  +  M  cos  k  z)(l  +  r  M  cos  k  z)  ■  1  +  •"  r  +  (1  +  r)  M  cos  k  z 

8  g  2  g 


1  2 

+  •r  r  M  cos  2k  z 
2  g 


(7.11) 


Equation  7.9  becomes 


+A 


-  n 


<N 


+A 


-  n>  + 


(7.12) 


where 


<N  -  n>  =■  N  (g/n  -  1  -  n) 


(7.13a) 


<n>  =  (<Yj>/Y^)<N^>(1  -  -  r  M^)  <<  N_^,  for  <Yj>  <<  N_^Yj./<Ng>  (7.17b) 

We  consider  only  the  low-lrradlance  limit  g  «  1.  Equation  7.14  and  g  «  1 

give 

<Yj>  «  N^V^e^  (1  +  ^  r  ,  (7.18) 

In  which 

<n>  -  (<Yj>/Yj.)<N^>  (l  +  ,  (7.19) 

t 

which  agrees  with  Equation  3.2  In  the  limit  of  negligible  depletion,  <Ng>  = 

NgO,  and  r  <<  1.  Equations  7.14  and  7.19  give 


(N^.  -  n)  5 
+A  z  +0 


(1  +  r)  M 


,  ~  cos  k  z  + 

1  +ir  8 


i  r  m2 


cos  k  z 


Substituting  Equation  7.20  into  Equation  7.2  and  integrating  gives 


E  -  E 


(1  +  rl 


1  + 


rM 


^  sin  k  z 
2  g 


rM^/4 


sin 


2k  z 
g 


M 


(7.21) 


where 


i 


E  “  4iTeN  /ek 


(7.22) 


For  negligible  depletion  the  Inequality  r  <<  1  is  valid,  l.e., 

<<  NgO  from  Equation  7.15.  With  Ng^  »  from 

Equation  2.16  and  M  cos  kgZ  from  Equation  7.20  with 

n^  «  the  term  Ng5/<Eg>  in  Equation  7.9  is  negligible  for 

either  «  NgO  or  small  modulation.  In  this  case  of  r  «  1,  ' 

Equation  7.21  gives 

E  -  E  M  sin  k  z,  for  g,  r  «  1  .  -  (7.23) 

q  g  .  - - 

and  Equation  7 . 14  and  N^.  -  N+q  +  give 

N+  -  (1  +  M  cos  kgZ),  ^  for  g,  r  «  1  ,  (7.24) 

where  g  <<  1  is  the  small-lrradlance  limit  of  Equation  7.17b  and  r  <<  I  Is  the 
negligible-depletion  limit.  Thus,  for  M  «  1  not  satisfied,  depletion  Is  not 
negligible  (N+^/N+q  =  M  cos  kgZ  Is  not  much  less  than  unity). 

The  result  In  Equation  7 .8  obtained  formally  (and  Incorrectly)  by  using 
n  =  <n>  In  Equation  7.1  Is  negligible  with  respect  to  the  correct  result  In 
Equation  7.23,  l.e., 

E  /E  -  k  A  GY  MirepN^  -  3.7  x  lO”^  «  1  .  (7.25) 

p  q  g  p  r  +0 

From  Equation  5.5,  we  have  Cj  ■  E^j/E^.  Thus,  Eq  «  Ej  for 
»  1.  With  Equation  7.25  this  gives 

Ep  <<  E^  <<  E^  for  long  transport  distance  .  (7.26) 

The  following  check  on  the  consistency  of  Equations  7.23  and  7.1  Is  not 
essential  to  the  understanding  of  the  control  results.  The  small  deviations  In  n 
from  the  large  term  <n>  are  retained  and  shown  to  be  small  and  to  account  for  the 
result  In  Equation  7.23.  The  field  E  in  Equation  7.23  comes  from  the  first  term 


in  Equation  7.1.  The  small  deviation  in  n  from  <n>  must  be  retained  in  Equation 
7.1  in  order  to  obtain  Equation  7.23.  First  consider  the  short-circuit  case  of 
Equation  7.1c.  The  second  term  in  Equation  7.1c  la  negligible  as  shown  in  Equation 
7.26.  Setting  n  =»  <n>  (1  +  n^  M  cos  kgz)  in  Equation  7.1c  and  equating  the 
resulting  E  to  the  E  in  Equation  7.23  gives,  with 
Dkg/w  -  Ej, 


E  ,n  M  sin  k  z  »  E  M  sin  k  z 
d  c  g  q  g 


Solving  for  n^.  gives 


n  -  E^/E.  =  AireN  '  u/eDk  ^  H  1/E  «  1 

c  q  a  +u  g  T 


(7.27) 


where  is  Debye-length  parameter  of  Kukhtarev  and  coworkers. 

For  the  open-circuit  case  of  Equation  7.1c  defining  n  =•  <n>(l  +  nj,g), 
where  n^,g  ■  n^  M  cos  kgZ  +  ng  M  sin  kgZ,  and  using  (1  +  n^g)"^  2 
1  -  n^g  for  ngj,  <<  1  and  setting  the  resulting  E  equal  to  the  E  in 
Equation  7.23  gives 


E  ,n  M  sin  k  z  -  E  .n  M  cos  k  z  +  E  M  cos  k  z  -  E  n  M  cos  k  z 
dc  gds  gP  gpc  g 

-  E  n  M  sin  k  z  ■  E  M  sin  k  z 

P  3  g  q  g 


Setting  the  sum  of  the  coefficients  of  the  cos  k^  terms  equal  to  zero  gives 


n 

s 


"p^d 


A  k  /A  A  ^ 
P  g  d  g 


«  1 


(7.28) 
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where  Apkg  «  1  Is  always  satisfied  and  »  1  Is  satisfied  in  the  long 

transport  limit.  Setting  the  sum  of  the  coefficients  of  the  sin  kgZ  terms  equal 
to  zero  gives 


n 

c 


E  +  E  n 
q  P  g 


«  1 


where,  Eq/Ejj  «  1  from  Equation  7.27,  ng  «  1  from  Equation  7.28,  and 
Ep/Ed  «  1  from  Equation  7.28. 


SECTION  8 


RESULTS  OF  KUKHTAREV  AND  COWORKERS  AND  OF  AMODEI 

Several  results  of  Kukhtarev  and  coworkers ^  and  of  Amodel^  are  discussed 
in  order  to  show  the  relation  between  the  present  results  and  the  resdlts  of 
these  previous  Important  papers  and  moreover  to  resolve  inconsistencies  among  the 
three  sets  of  results. 

Kukhtarev  and  coworkers  consider  only  the  steady  state.  We  set  their 
thermal  generation  factor  8  equal  to  zero  because  the  levels  of  interest  here  are 
too  deep  for  thermal  generation  of  free  carriers  to  be  Important.  We  also 
neglect  the  small  correction  to  M  for  the  case  of  kgA  «  1  even  though  It 
should  be  simple  to  obtain  this  correction  from  the  present  model.  The  relations 
between  Kukhtarev* s  notation  and  ours  are  as  follows: 


t  -Yi 

*  V'"40 
T  >  kgT 


Cos  (2k  X  +  <fr)  cos  k  z 

■K  g 

X  u 


uT/|e|  ■••D 


C-  (E  +  E  )/E 
E  app  p  q 


Their  condition 


Yj  «  Y,n/N^ 


(8.1) 


is  valid  only  in  samples  having  many  more  filled  levels  than  eiqpty  levels,  that 
Is, 


(8.2) 


A-A9 


This  can  be  seen  as  follows.  Setting  SN+Zst  ■  0  and  -  Np  +  N4^  in 
Equation  2.9  and  solving  for  gives 


Y  <Yj/Y,)N„H^/« 


(8.3a) 

# 

(8.3b) 


where  the  approximate  equality  in  Equation  8.2b  is  valid  when  the  Kukhtarev-  and 
coworkers  unequallty  in  equation  8.1  is  valid.  Again  setting  BN^/St  -  0  in 
Equation  2.9,  but  retaining  N^,  and  using  Equations  2.14  to  2.16  gives 


n 


/  \ 

V  ""“7  V'  “+/"«) 


(8.4) 


where 


Equation  8.3b  can  be  rewritten  as 


n 


”unsat^^D^^eO^  1  +  N. ./N. 


+A'  +0 


(8.6) 


The  exact  result  in  Equation  8.4  reduces  to  the  approximate  result  in  Equation 
8.6  if  Np  =  Ngo*  Then  N+q  «  N^q,  the  factor 

N+A/Ngo  is  negligible  in  the  negligible-saturation  limit,  and  the 
factor  Np/NgQ  s  1  in  Equation  8.5. 

Equation  9  of  Kukhtarev  and  coworkers  with  m  =»  M  is  in  agreement  with 
Equations  3.12,  4.8,  and  5.3.  Thus,  our  results  and  theirs  agree  for  the  short- 
transport  distance  limit  of  kgA  «  1.  Their  Equation  9  is  valid  even  when 
their  approximation  «  Nq  is  not  satisfied. 
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Next  consider  the  limit  k^A  »  1.  From  Equation  8.6  with  the  Kukhtarev 
and  coworkers  approximation  N^/N^q  s  1,  in  order  to  have  n  =  <n>  we  must 


N.  *  ■  M  cos  k  z 
+A  +0  g 


(8.7) 


so  that  the  factor  (1  +  N4.^/N^)~^  in  Equation  8.6  cancels  the  factor 
(1  +  M  cos  kgZ)  in  (see  Equation  8.5).  Notice  that  in  this  case 

of  N^^  «  Nj),  the  steady  state  is  not  reached  until  N^^  is  "fully 
saturated"  for  M  5  1,  that  is, 


N^  -  (1  +  M  cos  kgZ)  . 


(8.8) 


r  •**  ♦••  •  •  a  ' 


The  charge  density  from  Equation  2.5  is 


■■  ^'V-*  V 


P  -  fi(\^  “  '!)• 


(8.9) 


With  n  «  N+a,  substituting  Equations  8.9  and  8.8  into  Equation  2.11  and 
integrating  gives 


E  -  <E>  -  E  M  sin  k  z 

q  g 


(8.10) 


r. V  ^  ' 


where  <E>  is  the  obvious  constant  of  integration  and 


%  ■  - 


(8.11) 


.'-•'•.‘■-'a .L-' 


is  In  agreement  with  Equation  7  of  Kukhtarev  and  coworkers. 
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Finally,  notice  that  Equation  8.10,  which  Is  Equation  7  of  Kukhtarev  and 
coworkers.  Is  Incorrect  for  the  case  of  an  applied  field  only  and  an  open  circuit 
because  Equation  2.10  with  J  ••  D  ■  p  =■  0  gives  E  *  0  for  n  ^  0.  The  diffusion 
term  In  Equation  2. 10  and  the  small  deviation  oq  in  n  from  <n>  are  needed  to 
establish  the  field  In  Equation  8.10  even  throu^  D  does  not  appear  In  the  result 
(because  the  1/D  In  n^  cancels  the  explicit  factor  D  in  Equation  8.13). 

Amodel^  obtained  the  diffusion-dominated  result 


Ej  M  sin  k  z 

2  _ 8 _ 

“  1  +  M  coa  k  z. 


(8.12) 


even  In  the  limit  of  A^^fkg  »  1.  This  result  disagrees  with  the  present 
result  in  Equations  7.21  or  7.23  and  with  the  result  in  Equation  8.10  from 
Kukhtarev  and  coworkers.  The  reason  for  Che  disagreement  Is  Chat  Amodei  assumes 
that  Ngo  (his  N)  and  N+o  (in  his  t  ■  l/7j-)  are  constant.  Then 


n  *  <n>  (1  +  M  cos  k  z) 

g 


(8.13) 


rather  than  n  *  <n>,  as  obtained  here  by  Kukhtarev  and  coworkers 

This  result  in  Equation  8.12  is  not  correct  in  general  because  Che  changes 
in  Che  values  of  and  Ng  cannot  be  neglected  in  Che  steady  state  in 
general.  For  example,  was  fully  saturated  in  Che  case  above,  according 
to  Equation  8.8.  Equation  8.13  is  physically  lo^lausable  for  this  case  of  kgA 
»  1  because  the  transport  distance  Is  greater  than  the  distance  kg  over  which 
CI>  M  cos  kgZ  changes  substantially  and  the  electrons  cannot  remain  in  the 
region  In  which  they  are  generated  if  they  diffuse  or  drift  such  a  long  distance 
before  recombining. 
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The  result  in  Equation  8.12  is  also  inconsistent  with  the  value  of  E 
obtained  from  Poissons  equations.  For  N+q  =  constant,  =  0,  and 
Equation  2.5  gives 


P  *  -en  “  -e<n>(l  +  M  cos  k  z)  . 

S 


(8.14) 


Substituting  the  z-dependence  term  into  Equation  2.11  and  integrating  gives 


-(4ire<n>M/ek  )9ln  k  z 


which  does  not  agree  with  Equation  8.12.,  Finally,  the  constant  term  in  p  in 
Equation  8.14  does  not  vanish,  as  required  by  quasi-damage  neutrality. 


SECTION  9 


PHOTOVOLTAIC  EFFECT 

In  photovoltaic  materials,  a  space-charge  field  is  believed  to  be  generated 

# 

by  asymmetric  photoemlsslon  of  conduction  electrons.  ^  In  such  ferroelectric 
materials  as  LiNb03  and  BaTi03,  ions  can  have  potentials  that  are  asymmetric 
along  some  direction.  For  example,  the  titanium  ion  in  BaTi03  is  displaced 
from  the  center  of  the  lattice  cell  in  the  ferro-electric  state.  Thus,  the 
photogenerated  electrons  are  emitted  preferentially  in  one  direction,  which  gives 
rise  to  the  photovoltaic  field  Ep  as  follows:  By  definition  of  the  momentum 
transport  time  Tj^,  the  initial  momentum  is  lost  in  time  on  the  average. 

The  displacement  of  the  photogenerated  electrons  during  the  time  has  a 
typical  value  ip  =  vxy^  =  (10^  cm/s)(10"^‘*  s)  =  1  nm  (=•  10  A),  All  of  the 
electrons  are  not  photoemltted  along  the  same  direction,  but  there  is  a  greater 
probability  of  emission  along  the  photovoltaic  direction.  Thus,  an  effective 
distance  Ap  is  defined  as 


P+V  ■  "-‘b-  '  *D  -  Ph 


(9.1) 


where  p+  and  p_  are  the  probabilities  of  emission  in  the  +  and  -  directions, 
respectively. 

A  typical  value  of 

A  =-0.08  nm  (-  0.8  A)  (9.2) 

P 


is  consistent  with  typical  observed  photovoltages,  as  shown  below.  Fouatlons  3.1 
and  9.2  give  »  Ap/ip  =  0.08  nm/1  nm  =■  0.08  as  the  relative 
probability  for  emission  in  the  positive  photovoltaic  direction.  The  short 
distance  of  0.08  nm,  which  is  less  than  an  atomic  spacing,  Is  not  the  actual 
transport  distance,  which  would  be  unrealistically  short.  Rather,  It  is  the 
weighted  distance,  as  in  Equation  9.1. 
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In  the  literature, ^  it  is  common  practice  to  model  the  photovoltaic  effect 
by  including  the  term 


J  =■  p.  le 
~p  k 


(9.3) 


in  the  current,  where  pf^  is  a  constant  and  c  is  a  unit  vector  in  the  opposite 
direction  from  the  preferred  electron  emission  direction.  The  following 
derivation  shows  that  Equation  9.3  is  consistent  with  Glass’  preferential 
emission  model  If  the  constant  Pi^  is  replaced  by  p  Ndg,  where  p  Is  a  constant 
and  Is  the  trapped  electron  density: 


J  =  p  N  le 
D  er 


(9.4) 


Equation  9.4,  but  not  Equation  9.3,  is  equivalent  to  Glass'  expression  J  =  ical 
because  the  absorption  coefficient  a  is  proportion  to  Ng.  Equations  9.3  and 
9.4  are  equivalent  In  the  Limit  of  no  depletion  (Ng  ■  NgO  *  constant). 

Equation  9.4  is  obtained  by  considering  the  generation  rate  of  electron  by 
the  field 


(3n/3t)_  -  (o^/h'*))I(z  -  A  )N  (z  -  A  )  .  (9.5) 

ID  pep 


The  electrons  generated  from  optical-donor  levels  at  z  -  Ap  appear  at  the 
displaced  position  z.  To  Taylor  expand  in  Ng  gives 

(3n/3t)T  ■  a_I(z)N  (z)hu)  -  3(A  o  I(z)N  (z)/hu))/9z  , 

t  U  e  p  D  e 

which  is  of  the  usual  form 


3n/3t  -  a„IN  /Tio)  -  e  3J/3z 
u  e 
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T 


If  We  set 


J  -  J  -  pN  Ic 
~  ~p  e 


(9.6) 


where  6  ■  -“z  and 


p  -  eA^Ojj/ha) 


(9.7 


We  shall  use  Jp  in  Equation  9.6,  rather  than  the  previous  expression  in 
Equation  9.3. 

The  photovelcale  field  that  develops  In  the  steady  state  under  open-circuit 
conditions  (that  Is,  J  <n>euE  +  pNg<l'>  ■  0)  is: 

E  -  pN  <I>/<n>eu  -  A  y  /p  •  (9.8) 

P  e  P  r 


The  time  constant  for  Ep  to  approach  the  steady  state  value  in  Equation  9.8  is 

the  dielectric  relaxation  frequency  defined  in  Equation  3.7,  as 

usual. 

A  small  photovoltage 


E  -  4nE<n>A  /e  -  (t^,/y  )E  (9.9) 

psm  p  dl  r  p 

is  developed  in  the  short-time  r^,  as  seen  In  the  following  model.  Figure  4 
schematically  illustrates  the  charge  distributions  resulting  from  the 
displacement  of  the  electrons,  denoted  by  -,  in  the  positive  z  direction  by  a 
distance  Ap.  The  resulting  net  charge  density  In  Figure  4  Is  approximated  by 
5-functlon3. 


p  ■  e<n>A  f6(z)  -  6(z  -  1) ] 
P 


1  1368-4 


0  JL 

DISTANCE  (Z) 


(a)  IRRADIANCE  DISTRIBUTION 

O  O  O  Q 
O  0  0  © 

©  ©  ©  ©  _ 

©  ©  0  © 

DISTANCE  (Z)  Apv 

(b)  RESULTING  CHARGE  ARRANGEMENT 


(C)  RESULTING  CHARGE  DENSITY.^ 


Figure  4.  Model  for  Che  photoelectric  effort  showing  the 
idealized  charge  distribution  resulting  from 
net  electron  emission  in  the  positive  z 
direction. 
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where  Che  electron  concentration  <n>  develops  in  a  characteristic  time 
Integrating  Poisson’s  equation,  edE/dz  =•  Airp  gives, 

E-E  forO<z<£ 

psro 

••  0  for  all  other  values  of  z 

where  E^g^  Is  given  by  Equation  9.9.  Since  Ep  developes  in  time 
Tjjj  and  Epgm  develops  in  time  Tj.,  It  is  reasonable  that  Epg^j 
“  ('^r/'^dl^^p  *  (Ydi/Yr)E,  as  in  Equation  9.9. 
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SECTION  10 


MODEL  OF  DIELECTRIC  RELAXATION 

The  following  model  illustrates  the  significance  of  the  dlelectrl,c 
relaxation  frequency. 

=  1/t^  ■  4iro/e  -  4TTe<n>/e  ,  (10.1) 


which  is  the  time  required  to  change  or  discharge  a  sample.  Consider  a 
dielectric  sample  in  an  externally  applied  electric  field  E^pp,  with  open 
clrdult  conditions,  as  Illustrated  in  Figure  5.  At  time  t  =  0,  the  sample  is 
undharged.  The  electric  field  E^pp  causes  a  current 


J 


0 


eu<n>E 

app 


at  time  t  ■  0.  The  divergence  of  the  current  at  the  end-faces  of  the  sample 
causes  the  end-faces  to  charge  up,  according  to  the  continuity  equation 

3p/3t  -  -3J/3dz  (10.2) 

where  p  is  the  charge  density.  For 

J(2,t)  -  -eu<n>E(t)r(z)  (10.3) 


where 


r(z)  ■  1  for  0  <  z  <  i 

=  0  for  all  other  values  of  z 


(10.4) 


Equation  10«2  gives 


3p/3t  =  -eu<n>E(t)A(z) 


(10.5) 
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Crystal  with  an  applied  field  with 
open-circuit  conditions  at  the  faces 
of  the  crystal  showing  the  current 
density  J,  electric  field  E,  and 
charged  faces  of  the  crystal. 


whefe 


A(z)  »  5(z)  -  6(z  -  1)  . 

Hatfe  <5{*)  Is  Che  Dlrsc  6-functlon,  and  dr(z)/dz  »  A(z). 

From  Poisson’s  equation, 

p  =•  e(ATr)~^3E/at  -  e(d4it)"^3(E(t)r(z) ] /3z 
=■  ed(4ii)"^E(t)&(z)  . 


(10.6) 


(10.7) 


Taking  the  time  derivative  of  p  with  respect  to  time  in  Equation  10.7  and 
substituting  the  result  Into  Equation  10.5  gives 

I 

3E(t)/at  -  -(4ir  ep<n>/e)E(t)  .  (10.8) 


For  E(t)  =*  K-oo  at  t  «  0,  the  solution  to  Equation  10.8  is 


E(t)  -  E  e 
app 


(10.9) 


tfhePe  is  the  dielectric  relaxation  frequency  defined  In  Equation  10.1. 

The  model  Illustrates  that  the  time  constant  for  the  development  of  the 
final  charge  density  (p  ■  Po^(z))»  which  corresponds  to  the  charged  end  faces 
of  the  sample)  is  the  dielectric-relaxation  time,  the  electric  field  and  current 
decay  to  zero  as  the  end  faces  of  the  sample  charge  up. 

U^lng  the  same  model,  but  with  Egpp  =•  0  and  p  “  po4(z)  at  time  t  *  0, 
it  is  easy  to  show  that  an  initial  charge  decays  with  time-constant  l/Ydi^ 

E  -  Eq  e  .  (10.10) 


The  current  and  electric  field  decay  to  zero  as  the  end-face  charge  decays  to 
zero* 
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SECTION  11 


SUMMARY  OF  RESULTS 

A  list  of  symbols  and  typical  values  of  parameters  Is  given  In  the  Appendix 
Key  results  for  the  grating  electric  field  Eg,  which  Is  the  z-dependent  part  of 
E  =•  Bgpp  +  ffgc  Cwhlch  excludes  the  electric  field  of  the  laser 
lrradlance)i  are  listed  below. 

Tile  corresponding  value  of  the  grating  Index  Is  given  by  Equation  1.1  and 
repeated  here 


n 


g 


1 

2 


n  ^rE 
r 


(1.1) 


and  the  corresponding  grating  efflclenc^y  Is  given  by  an  expression  such  as 
Equation  1<5  and  repeated  here,  l.e.. 


n 


-ad/co8 

e 


e 


sin 


iTn  n  d 

_ E_g _ 

2X  cos  6 


(1.5) 


The  notation  for  the  various  cases  of  E  below  Is  as  follows: 
E^ppi  applied  field  dominated 
Dt  diffusion  dominated 

p:  photovoltaic-affect  dominated 
tt  short-time  solution  (Tj.  t  Tjl) 
sst  steady-state 


^Ing- 

^sht* 

ngbl: 


kgA  »  1 
kgA  «  1 

negligible  depletion  and  saturation. 


Combining  Che  results  of  Sections  3  through  7  gives  the  following  key 
results: 

From  Equations  3.18,  4.9,  5.7  and  7.7,  for  kgA  «  1 


E  2 


"dl^ 


[(e 


app 


E  }  M  cos  k  z  +  ■"  E,  M  sin  k  z 
P  g  2  d  g 


C  M  cos  k  z  -  J-/ep<n>  +  E . 

_ g _ S _ a 

1  +  M  cos  k  z 

g 


for  A 


short 


> 


M  sin  k  z 
_ g_ 


negligible 


) 


(11.1) 


where 


jQ/ep<n>  -  0 

-  E 

P 

-  E  -  E 

P  app 


I 

for  open  circuit 

for  short  circuit 

for  E  with  ohmic  contents 
app 


(11.2a) 

(11.2b) 

(11.2c) 


From  Equations  3.28,  4.11,  5.9,  and  7.20,  for  kgA  »  1, 


E  «  Y,-te  ^  (E.  +  E  )  M  sin  k  z 
di  '  A  p'  g 

.  _  I  (1  +  r)  M  .  .  .  r  M^/4  ^  I 

+  E  - ;  — r  sin  k  z  + - - r  sin  2k  z 


(11.3) 


The  results  for  the  phase  shifts  between  the  grating  electric  field  and  the 
gracing  Irradance  <I>  M  cos  kgZ  are  summarized  for  the  various  combination  of 
long-  and  short-  electron-diffusion  distances  and  drift  (In  applied  or 
phocovoltald  field)  and  diffusion-limited  transport  in  Table  2. 


T^ble  2.  Spatial  Phase  Shift  Between  Che  Gracing 
(l.e.,  Photofractive  Electrtic  Field  E)  and  Irradlance  I. 


kgA  «  1 

kgA 

»  1 

Drift 

Diffusion 

Drift 

Diffusion 

Phase 

TT 

-ir/2 

it/2 

n/2 
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I  Introduction 

Some  of  the  most  attractive  materials  for  practical  applications  of 

degenerate  four  wave  mixing  (DFWM)  are  photorefractive  crystals  such  as 

0 

lithium  niobate  (LiNbO^),  potassium  niobate  (KNbO^),  barium  titanate 
(BaTiOj),  strontium  barium  niobate  (SBN),  bismuth  silicon  oxide  (BSO)  and 
many  others.  These  crystals  have  large  nonlinearities  so  that  low  power 
cw  operation  with  small  interaction  volumes  is  possible.  They  are  solids 
and  operate  at  room  temperature.  Large  diffraction  efficiencies  and  phase 
conjugate  reflectances  have  been  achieved  across  the  visible  spectrum. 
Finally  under  some  conditions,  the  nefractive  index  interference  pattern 
in  these  crystals  is  shifted  90®  from  the  intensity  pattern.  This  shifted 
grating,  which  at  present  is  unique  to  photorefractive  materials,  allows 
construction  of  novel  devices  using  four-wave  mixing  [1-3]. 

In  order  to  analyze  the  performance  of  photorefractive  materials  for 
DFWM  applications,  it  is  first  necessary  to  understand  the  formation  of  a 
space  charge  grating  due  to  the  Interference  of  two  waves  within  the  medium 
In  the  analysis  presented  here,  we  assume  that  the  amplitudes  of  the  two 
waves  are  known  in  a  given  volume  element.  The  irradiance  due  to  this 
interference  of  those  waves  is  given  by 

I  =  I^  (1  +  m  cos  ) ,  (1) 

where  m  is  the  modulation  depth,  kg  is  the  grating  wave  vector,  and  z  is 
the  direction  normal  to  the  grating  planes.  If  the  two  input  waves  are 
co-polarized  and  have  intensities  L  and  I  ,  then 


Grating  formation  in  photorefractive  materials  can  be  characterized 
in  terms  of  three  parameters  -  the  steady -state diffraction  efficiency, 
the  write  or  erase  time,  and  the  write  or  erase  sensitivity.  In  the 
following  section  we  review  results  for  these  performance  parameters 
derived  from  the  theory  of  photorefractive  materials  developed  by 
Kukhtarev  et  al .  [6-11].  In  section  3  these  results  are  applied  to 
several  photorefractive  materials  in  order  to  predict  differences  in  four 
wave  mixing  behaviour.  In  Section  4  we  discuss  materials  improvement 
for  optimal  performance.  ' 

II  Basic  Equations 

We  start  with  the  basic  equations  used  by  Kukhtarev,  Vinetskii  and 
coworkers  [6-11]: 

K 

in  _ =  .  i  ii  (continuity)  (3) 

at  at  e  az 

aN"*" 

=  sKN^-nJ)  -  equation)  (4) 

j  =  eunE  -  kgTu|j  (current  equation)  (5) 

II  =  (n+N^-Np)  (Poisson's  law)  (6) 

where  n  is  the  electron  number  density, 

is  the  total  number  density  of  dopants  (Fe  for  example,  ionized 
and  neutral;  Nq  is  a  constant), 

Np  is  the  number  density  of  ionized  donors  (Fe^^  for  example)  , 
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N 

A  is  the  number  density  of  negative  charges  that  compensate 
for  the  charge  of  Np  under  dark  conditions  (N^  is  a  constant 
of  the  crystal), 

Np-Np  is  the  number  of  neutral  (filled)  donors, 
j  is  the  current  density, 

s  is  the  photoionization  cross  section, 

I  is  the  irradiance  of  the  optical  wave, 

is  the  two-body  recombination  rate  coefficient, 

M  is  the  mobility,  ' 

E  is  the  total  electric  field, 
kg  is  Boltzmann's  constant, 

T  is  the  temperature, 

e  is  the  charge  on  the  electron, 

e  is  the  static  dielectric  constant. 


One  important  parameter  which  has  a  strong  influence  on  the 
performance  of  photorefracti ve  materials  is  the  reduction  ratio,  defined 


as 


R  _  ^D'^A  _ 

^  ■ 


This  parameter  is  just  the  ratio  of  filled  empty  donor  states. 
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Equations  (3)  -  (6)  may  be  easily  modified  to  include  dark  conducti¬ 
vity  and  the  photovoltaic  effect  [8].  The  equations  do  not  consider 
secondary  trapping  levels  or  photorefracti ve  centers  [12,  13]  or  holes 
[14],  all  of  which  may  be  Important  under  some  conditions.  These  equations 
appear  to  give  similar  results  to  the  recently  developed  hopping 
model  [4]. 

The  set  of  nonlinear  equations  can  be  solved  by  using  two  approxima¬ 
tions:  linearization  in  the  grating  modulation  depth  and  the  quasi -steady 
approximation  in  the  zeroth-order  ele|Ctron  number  density  [6].  With  these 
approximations  one  can  solve  for  the  grating  space  charge  electric  field 
in  the  write  and  erase  cases. 

The  grating  space  charge  field  during  writing  is  given  by 


where 


6E  =  mE^^  [-cos(kgZ  +  ij;)  +  e 


-t/T. 


'cos  (kgZ+UJgt+4;)] 


(7) 

(8) 


(9) 


2..  2, 


2.1/2 


^  =  arctan{(Eo/E^)[l  +  Ej^/E^  +  ^//(EjjEq)]} 


E  -  S'bT 

“  e  *  (di'f'fusion-linTited  space  charge  field) 


Bq  ~  applied  field  (normal  to  grating  planes) 


Bq  "  .  (space  charge  field  for  complete 

^  9 

separation  of  available  charge) 


and  the  time  scales  are  defined  by 


^di  ^  '^a*  (dielectric  relaxation  time) 


Te  =  ttJe"  ’ 


ykeTk^ 


,  (diffusion  time) 


■^4.  =  ’  ("’edified  recombination  time)  (17) 
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"1  =  <='o  *  Vo>  ■  "8) 

where  a  is  the  conductivity,  and  is  given  by 

a  =  n^eu 
0 


In  eqs.  (17)  -  (19)  n^  is  the  spatially-averaged  electron  density 
Under  normal  cw  conditions  these  approximations  hold: 


^+"<^di 


(20) 


(21) 


sI(Nn-N.) 


Under  uniform  illumination,  the  space  charge  field  during  erasure  is 
given  by 

-t/T 

6E  =  Eg  sin  (k^z  +  oi^t)  e  ®  ,  (23) 

where  Eg  is  the  initial  grating  amplitude  and  and  are  defined  above. 

There  are  several  interesting  features  of  the  expression  for  the 
response  time  x^  given  by  eq.  (8).  First  it  is  obvious  that  the  erase 
time  is  a  complicated  function  of  gra,ting  period,  applied  field  strength, 
temperature,  etc,,  in  addition  to  being  proportional  to  the  dielectric 
relaxation  time  x^^. .  In  the  classical  case  of  a  charged  dielectric  sample, 
(e.g.  a  capacitor)  the  response  time  is  just  x^^ ,  which  can  be  considered 
as  the  intrinsic  RC  time  constant  of  the  material.  The  added  terms  in 
eq. (8)  are  due  to  the  periodicity  of  the  charge  pattern  in  the  crystal; 
note  that  when  kg  -K),  eq  (8)  reduces  to  x^  =  x^j. .  Even  when  a  grating  is 
present,  there  are  important  limits  for  which  the  correction  terms  are 
small,  and  x^  =  Xj.j .  These  cases  will  be  discussed  later. 

Ill  PERFORMANCE  PARAMETERS 
A.  Response  Time 

The  complete  expression  for  the  response  time  x^  is  given  in  eq.  (8). 
As  indicated  earlier,  there  are  important  limits  when  x^  =  ^di ’ 
first  consider  the  dependence  of  x^.  on  important  material  parameters.  By 
combining  eqs.  (14),  (19)  and  (22),  we  obtain 


‘^di  "  ATreysIj, 

-/  )  1_  i 

Mfreys  '  1^  R 


(24) 


The  quantity  in  parentheses  in  eq.  (24)  can  be  considered  as  a  material 

figure  of  merit,  having  units  of  energy/ unit  area.  Once  a  material  is 

chosen,  the  only  experimental  control  over  this  quantity  is  through  its 

( 

variation  with  temperature  (through  y  and  Yq)  with  wavelength 
(through  s).  However,  the  quantity  R  can  be  varied  over  many  orders  of 
magnitude  by  reduction  or  oxidation  of  a  given  sample.  It  is  clear  from 
eq.  (24)  that  increasing  R  through  chemical  reduction  can  significantly 
reduce  ,  Values  of  R  as  high  as  10  are  not  hard  to  achieve  in  many 

3 

cases,  and  in  some  materials  (e.g.  BSO)  values  of  R  on  the  order  of  10 
are  typical  [13].  Note  also  the  dependence  in  eq.  (24),  indicating 
the  direct  relation  between  speed  and  irradiance. 


Returning  to  eq.  (8),  we  now  wish  to  consider  how  the  simple 
relationship  x^  “  modified  by  the  presence  of  a  volume  space  charge 

grating  in  the  sample.  By  examination  of  the  time  constants  entering  into 
eq.  (8),  we  find  that  x^  =  x^^  is  still  a  valid  approximation  when 
x_^«XQ,Xg.  This  condition  is  approximated  when  either  of  the  following 
conditions  hold:  (1)  the  crystal  is  oxidized  (x^  is  minimized),  or  (2) 
the  applied  electric  field  is  small  (large  x^),  and  diffusion  is  slow 
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(large  due  to  low  temperature  or  large  grating  period.  This  is  shown 
more  clearly  in  figures  1-3,  which  are  log-log  plots  of  t  (normalized  to 

C 

Tdi  for  R  =  1)  as  a  function  of  the  reduction  ratio  R,  for  several  values 
of  applied  electric  field.  The  relationship  corresponds  to  a 

straight  line  with  slope  -1.  We  see  first  that  this  relationship  is  most 
closely  approached  for  E  0,  especially  when  R  <  100.  When  E  is 
zero,  there  is  an  increase  in  especially  in  highly  reduced  crystals, 
and  for  small  valves  of  the  grating  period  Ag.  As  pointed  out  by  Kukhtarev 
[  6],  there  are  also  certain  condition's  (low  mobility,  high  temperature) 
where  an  increase  in  E  can  lead  to  a  decrease  in  x^.  However,  these 
conditions  are  not  generally  attained  in  typical  experiments. 

B.  Sensitivity 

Several  definitions  of  sensitivity  have  appeared  in  the  literature. 

The  two  most  common  are:  (1)  the  energy  density  required  to  write  a 
grating  with  1  percent  diffraction  efficiency,  and  (2)  the  incremental 
change  in  refractive  index  per  unit  of  absorbed  energy.  In  both  cases, 
the  magnitude  and  direction  of  the  space  charge  field  must  be  known  in 
order  to  calculate  the  induced  refractive  index  change.  However,  there 
are  certain  experimental  situations,  (e.g.  domination  of  diffusion)  where 
the  space  change  field  is  nearly  invariant  from  material  to  material; 
thus,  no  new  information  is  provided  by  including  E^^  in  the  definition. 

In  other  situations,  the  space  charge  field  may  be  different  from  material 
to  material,  but  accurate  knowledge  of  its  value  may  not  be  important. 
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Figure  1.  Normalized  erase  time  versus  reduction  ratio  for  Indicated 
recombination  rate  and  grating  spacing. 
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In  this  work  we  would  like  to  propose  a  different  definition  for  the 
write/erase  sensitivity.  Since  we  know  Is  proportional  to  ,  and 
^d1  define  the  write  sensitivity  to  be  the  energy 

density  required  to  write  a  grating  to  (1-e*^)  of  Its  final  value  of 
space  charge  field.  Similarly,  the  erase  sensitivity  Se  Is  the  energy 
density  required  to  erase  a  grating  to  e“^  of  the  Initial  space  charge 
field.  Since  in  this  mode  the  write  time  Is  equal  to  the  erase  time,  we 


can  write: 


S  =»  S 

w 


I 


I  T 


0  e 


(25) 


2 

Note  that  the  sensitivity  (in  J/cm  )  Is  numerically  equal  to  the  response 

2 

time  when  I  =  1  w/cm  .  In  the  limit  where  t„  =  x..,  the  sensitivity  S  is 
0  e  (ji 

given  by 


S 


GYr  1 

^  4Trevis  ^  R 


(26) 


Our  definition  results  In  a  very  simple  expression  for  S,  which  includes 
the  same  material  figure  of  merit  (given  In  parentheses)  which  appears  In  the 
expression  for  x^^ .  Of  course  since  no  information  on  the  Induced  space 
charge  field  or  refractive  index  change  Is  contained  in  eq.  (26),  some  care 
has  to  be  taken  when  optimizing  the  sensitivity.  For  example  a  low  value 
of  G  Is  favored  for  improved  sensitivity,  but  materials  with  low  values  of 
e  have  small  values  of  electro-optic  coefficient,  and  may  thus  be  Inefficient. 
Of  course  this  is  accounted  for  in  the  other  definitions  of  sensitivity,  but 
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other  information  is  lost. 


C.  Efficiency 

The  diffraction  efficiency  for  Bragg  scattering  from  a  volume  phase 
grating  is  given  by 


n  =  e 


al 


COS0  . 
sin 


2 


/TtAnL-  \ 
'Acos0^  ’ 


(27) 


where  a  is  the  absorption  coefficient,  L  is  the  interaction  length,  0  is  the 
bragg  angle,  and  An  is  the  amplitude  'of  the  spatially  varying  refractive 
index.  In  materials  with  aL«l  and  with  An  «X/L,  we  obtain 


In  a  ferroelectric  with  the  grating  wave  vector  along  the  c-axis,  we  can 
write 


An  =  1/2  n^r  (29) 

where  the  electro-optic  coefficient  r  is  equal  to  r^^  or  r^^  depending  on 
the  polarization  direction  of  the  readout  beam.  In  a  cubic  crystal  or  in 
a  ferroelectric  with  a  different  grating  orientation,  the  component  of  r 
may  be  different  and  other  terms  may  be  required  [4,  13].  For  the  simple 
case  considered  here,  we  have 


(30) 


^zxcose' 


■)' 


E 


2 

sc  ' 


In  experiments  in  which  no  applied  field  is  used  and  the  bulk  photovoltaic 
effect  is  absent  (e.g.  in  a  cubic  crystal)  or  is  suppressed  (no  external 
current  path  supplied),  then  only  diffusion  can  contribute  to  the  space 
charge  field.  If  we  further  assume  that  as  is  typical  in  all  but 

highly  reduced  materials,  then 


E 


sc 


1 


(31) 


i  .e .  is  invariant  from  material  to  material.  Similarly,  if 
sc 

E^»E^»Ep,  then  again  the  space  charge  field  is  invariant. 

In  these  limits,  the  efficiency  for  a  given  material  is  determined  entirely 
by  the  electro-optic  coefficient,  except  for  the  smaller  variation  due  to 
then^  term.  Clearly,  the  most  favored  materials  for  achieving  high 
efficiency  are  those  with  large  values  of  electro-optic  coefficient. 

IV  APPLICATION  TO  DFWM  IN  SPECIFIC  MATERIALS 

The  concepts  of  response  time,  sensitivity  and  efficiency  are  all 
defined  with  respect  to  the  writing  or  erasure  of  a  grating  in  a  small 
volume  element  where  the  intensity  is  known.  The  modification  of  those 
parameters  for  DFWM  applications  is  not  straightforward,  and  requires  the 
solution  of  coupled  wave  equations  which  can  not  be  considered  in  the 
steady  state  if  the  time  dependence  is  desired.  In  the  typical  four- 
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wave  mixing  geometry,  the  solutions  show  that  the  nonlinear  reflectivity 
and  response  times  are  strongly  influenced  by  beam  coupling.  An  e>^mple 
of  such  an  analysis  in  the  steady  state  is  the  work  of  Fischer,  et  al  [15]. 

In  this  study,  we  will  continue  to  use  concepts  derived  for  grating 
formation  and  decay,  because  the  simple  expressions  we  have  obtained  are 
expected  to  be  qualitatively  accurate,  and  they  provide  useful  insight  into 
the  relative  performance  of  different  materials. 

In  Table  I  we  have  listed  measured  values  of  grating  response  time 

and  sensitivity,  as  well  as  DFWM  refVecti vi ty  for  a  number  of  materials. 

2 

Since  the  response  time  is  given  for  I^  =  1  w/cm  ,  it  is  numerically  equal 
to  the  sensitivity.  The  values  in  the  table  are  chosen  to  be  representa¬ 
tive  of  commonly  available  materials;  clearly  a  large  range  in  values  can 
be  obtained,  depending  on  orientation,  reduction,  etc.  The  values  marked 
with  an  asterisk  were  measured  at  Hughes  Research  Laboratories.  We  have 
also  listed  the  dominant  electro-optic  coefficient  for  each  material, 
because  of  its  relevance  to  the  efficiency  or  DFWM  reflectivity. 

The  most  sensitive  (and  fastest)  materials  in  Table  I  are  those 
which  are  most  highly  reduced  (e.g.  BSD,  BGO).  Furthermore  the  largest 
values  of  DFWM  reflectivity  are  observed  in  those  materials  with  the 
largest  electro-optic  coefficients.  This  is  because  in  most  DFWM 
experiments  (with  the  exception  of  BSD,  BGO),  no  applied  field  was  used, 
and  diffusion  was  the  dominant  transport  mechanism. 

The  important  conclusion  which  can  be  drawn  from  the  results 
presented  here  is  that  the  most  favorable  materials  are  those  with  the 


^  ^  .A.  U  . 
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Table  Comparison  of  photorefracti ve  materials 


largest  electro-optic  coefficient,  and  which  can  be  reduced  to  values  of 

R’>100.  For  example,  the  reduction  state  in  typical  coimercially  available 

BaTiO^  is  not  known,  but  values  of  R  much  greater  than  unity  are  not 

likely,  because  this  material  is  grown  in  an  oxygen  atmosphere.  Thus,  if 

2  3 

the  reduction  ratio  could  be  increased  into  the  range  10  -  10  ,  signifi¬ 

cant  improvements  in  the  response  time  and  sensitivity  could  be  obtained, 
with  little  or  no  loss  in  efficiency. 
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ABSTRACT 


The  charge  transport  model  of  photoref ractivity  is  used  to 
evaluate  four  figures  of  merit  that  can  be  used  to  characterize 
the  performance  of  photoref ractive  materials.  The  figures  of 
merit  are  the  steady-state  index  change,  the  response  time,  t!ie 
energy  per  area  to  write  a  one  percent  diffraction  efficiency 
grating,  and  the  photoref ractive  sensitivity,  or  index  change  per 
absorbed  energy  per  unit  volume.  These  indices  are  evaluated  as 
a  function  of  grating  period  and  applied  external  electric  field 
for  Bij2Si02Q,  a  fast  material  with  a  relativelly  small  electro¬ 
optic  coefficient  and  BaTi03,  a  slower  material  with  a  much 
larger  electro-optic  coefficient.  Methods  for  optimizing  tlie 
materials  are  discussed. 


A.  INTRODUCTION 


Photoref ractive  materials  such  as  lithium  niobate  {LiNb03), 
potassium  niobate  (KNb03),  barium  titanate  (BaTiOj),  strotitium 
barium  niolate  (SBN)  and  bismuth  silicon  oxide  (Bij^BiO^ij) 
attractive  new  candidates  for  real-time  optical  data  processing, 
(ODP);  in  the  past  few  years,  many  ODP  operations  have  been 
demonstrated  in  these  materials.  Reversible  holographic  storage 
in  LiNbO^  was  the  first  ODP  operation  demonstrated.  More 

recently,  edge  enhancement  ^  and  correlation-convolution 
operations  have  been  demonstrated.  Photoref ractive 

materials  have  also  been  used  for  phase  conjugation  via 
degenerate  four-wave  mixing. Finally,  under  some  conditions 
photoref  ractive  materials  can  be  used  as  optical  linage 
amplifiers.  ^ 

The  promise  of  photoref ract ive  materials  for  ODP  is  sub¬ 
stantial.  A  large  number  of  parallel  operations  can  be  processed 
in  a  single  crystal.  Photoref ractive  materials  can  store  optical 
holograms  for  time  scales  of  hours  to  years,  depending  on  the 
dark  conductivity.  Gains  up  to  a  factor  of  15  have  been  observed 
in  LiNb03,^^  and  gains  are  likely  to  be  higher  in  BaTiOj  and  SBN. 
Operation  on  nanosecond  time  scales  has  been  demonstrated  ^ 
with  Q-switched  pulses  from  doubled  YAG  lasers.  Finally,  the 
requirements  on  write  and  erase  energy  density  in  some 
photoref  ractive  materials  (e.g.,  BSO'^)  are  comparable  to  the  best 
photographic  plates  (50  to  100  pJ/cm^). 

There  are  two  major  factors  that  limit  widespread  appli¬ 
cation  of  photoref ractive  materials  at  present.  First,  some  of 
the  most  promising  materials  (e.g.,  BaTiOj,  SBN  and  KNbO.)  are 
not  widely  available  in  large  samples  with  high  optical  quality. 
Secon<i,  the  crystals  that  are  available  are  not  optimal  in  all 
respects.  For  instance,  in  order  to  demonstrate  'nigh  speed,  low 
write  energy,  long  memory  or  large  gain,  it  is  necessary  to  use 
several  different  types  of  materials.  The  purpose  of  this  paper 
Is  to  compare  photore f ract ive  materials  for  several  different 
applications  and  to  determine  optimal  properties  for  ODP.  To  do 
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this,  we  first  give  a  brief  description  of  the  charge  transport 
model  of  grating  formation  and  erasure.  Theii  we  define  four 
separate  figures  of  merit  for  OOP  operations  and  evaluate  them 
using  the  charge  transport  model.  Finally,  we  compare  two  photo- 
refractive  materials  (BSO  and  BaTi03)  that  each  have  certain 
optimal  properties  and  suggest  methods  for  improving  the 
materials.  For  a  discussion  of  additional  figures  of  merit, 
other  materials  and  applications,  the  reader  is  referred  to  the 
recent  comprehensive  review  by  Gunter. 


B.  CHARGE  TRANSPORT  MODELS 


Two  separate  models  for  describing  charge  transport  in 
photoref ractive  materials  have  been  formulated.  The  hopping 
model  '  assumes  that  carrier  transport  occurs  via  hopping  from  a 
filled  donor  site  to  a  neighboring  empty  trap.  The  band 
transport  model  ^  ^ ^ assumes  that  electrons  (or  holes)  are 
optically  excited  from  filled  donor  (or  acceptor)  sites  to  the 
conduction  (or  valence)  band  where  they  migrate  to  dark  regions 
in  the  crystal  by  drift  or  diffusion  before  recombining  into  an 
empty  trap  (see  Figure  1).  In  both  models  the  transported 
charges  result  in  an  ionic  space  charge  grating  which  is,  in 
general,  out  of  phase  with  the  incident  irradiance  (see  Figure  2) 
The  space  charge  grating  is  balanced  by  a  periodic  space  charge 
electric  field  in  accordance  with  Poisson's  equation.  This  space 
charge  field  modulates  the  refractive  index  through  the  electro¬ 
optic  effect.  The  band  transport  model  has  been  developed  to  a 
greater  deqree  and  is  in  wider  use  in  the  literature;  this  is  the 
model  which  we  will  use  for  our  analysis.  A  detailed  comparison 
)f  the  band-t ransport  and  hopping  model  is  beyond  the  scope  of 
this  paper. 


A  mathematical  de.'.cr  ipt  ion  of  the  grating  formation 
process  is  given  by  the  ecjuations  developed  most  fully  by 
Kukhtarev  ‘  '  f and  co-workers  (see  also  earlier  work  by 

Amodei ' ‘  and  others  ) : 
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Figure  1.  Schematic  o£  the  charge  transport  model  of 

photorefractivity.  Photo  ionization  excites 
an  electron  from  Fe'^  which  drifts  or  diffuses 
to  a  new  location  and  is  trapped  at  Fe^"*", 

This  leads  to  an  ionic  space  charge  and  an 
internal  electric  field  that  modulates  the 


(2) 
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3e  4ff  e 


<^z 


e  (n+N^-N+)  , 


(rate  equation) 
(current  equation) 
(Poisson's  law) 
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where 
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is  the  electron  number  density, 
is  the  total  number  density  of  dopants  ( Fe  for 
example,  ionized  and  neutral;  Nq  is  a  constant), 
is  the  number  density  of  ionized  donors  (Fe'"^  for 
example ) , 

is  the  number  density  of  negative  charges  that 

compensate  for  the  charge  of  N"*"  under  dark 

conditions  (N^  is  a  constant  oE  the  crystal), 

is  the  number  of  neutral  (filled)  donors, 

is  the  current  density, 

is  the  photoionization  cross-section, 

is  the  irradiance  of  the  optical  wave, 

is  the  two-body  recombination  ■'ate  coefficient, 

is  the  mobility, 

is  the  total  electric  field, 

is  Boltzmann's  constant, 

is  the  temperature, 

is  the  charge  on  the  electron,  and 

is  the  static  dielectric  constant. 


One  important  parameter  which  has  a  strong  influence  on  the 
performance  of  photoref ractive  materials  is  the  reduction  ratio, 
defined  as 
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N  -N  N 
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This  parameter  is  the  ratio  of  filled  to  empty  donor  states  under 
dark  conditions. 


The  optical  irradiance,  I,  in  Equation  (2)  is  assumed  to  be 
of  the  form 

I  =  lo  (1  +  m  coskgZ)  (5) 

where  m  is  the  modulation  index,  kg  is  the  grating  wavenumber 
and  Iq  is  the  average  irradiance.  For  two  copolarized  input 
waves  of  irradiances  and  I2/ 

2  di  12)^^^ 

m  =  c — — c -  .  (6) 

1 1  +  1 2 

For  erasure  with  a  uniform  beam,  m  =  0. 

Equations  (1)  through  (6)  must  be  supplemented  by  Maxwell's 
equations  for  the  optical  field  in  the  crystal  and  an  electro¬ 
optic  relation  between  the  space  charge  field  and  the  change  in 
the  crystal  refractive  index.  An.  Thus, 


An  =  1/2  n’  E 


where  n^  is  the  background  refractive  index  and  rgff  is 
the  effective  electro-optic  coefficient.  Equations  (1)  through 
(7)  may  be  modified  to  include  dark  conductivity  and  the  bulk 
photovoltaic  effect,'*^  but  these  effects  are  neglected  here.  The 
set  of  equations  also  does  not  consider  secondary  trapping  levels 
or  photoref ractive  centers. 

Equations  (1)  through  (4)  can  be  solved  with  two 
approximat ions :  1  inear i zat ion  in  the  grating  modulation  index 

and  the  quasi-steady  approximation  in  the  zeroth-order  electron 
number  density,  n,, .  " Using  these  approximations  and  omitting 
those  terms  that  are  important  only  for  high  irradiance  (MW/cm^) 
nansecond  pulses, ’  we  obtain  the  space  charge  field,  <5e: 


-t/i 


,  -cos  ( k  z  +  f  e 


'COSvK  z  +  'Ut  +  '*') 


The  transient  response  has  an  overdamped  oscillatory  behavior 
with  a  response  time, 

-  <1  *  *  (’r/'e'" 


e  di  [I  +  *  'r/'d>  *  ''r/'e> 


and  a  frequency, 


i  ^  ^di^  I  ^  ^ 

’di  +  ’7v^'’r/’e' 


In  addition,  the  amplitude,  Eg^,,  is  given  by 


sc  q  E^) 


and  the  phase  shift,  V,  between  I  and  6e  is  given  by 


with 


Ep,  E_  " 

tan  .  =  ^  (1  +  ^  .  -— )  , 
o  o  D  q 


=  applied  field  (normal  to  grating  planes). 
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(diffusion  field),  and 


4TTeN^ 

- — —  (limiting  value  of  the  space  charge  field) 

g 


Finally,  the  time  scales  are  defined  by 


‘di  ^  T^'eiTn”  (dielectric  relaxation  time). 
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(drift  time) 
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(ditfusion  time). 


(recombination  time),  and 


- + — T^rT"  (inverse  o£  sum  of  photoproduction  and  ion 

^0  r'^o  recombination  rates). 


and  the  zeroth-order  electron  number  density  is  given  by 

_SI(N^ 

^0  Y  N  '  ( i 3 ) 

R  A 

For  erasure  by  uniform  illumination  the  field  is  given  by 

-t/T 

■Se  =  E  cos  k  z  e  ^  ,  (14) 

where  Eg  is  the  initial  amplitude  and  ig  is  same  response 
time  obtained  for  writing  a  grating. 

FIGURES  OF  MERIT  FOR  OOP 
1 .  Steady-State  Index  Change 

The  steady-state  change  in  the  refractive  index  (also  called 
maximum  refractive  index '^^)  is  defined  as  the  index  change  reached 
after  illumination  for  a  time  long  compared  to  the  response  time, 

"p.  The  index  is  related  directly  to  the  space  charge  field  through 


'''n  =  1/2  n.  ^  r  p,-  E 

ss  '  b  eft  sc 


The  background  index,  nj^,  is  typically  2.5  for  most 
photore f ract i ve  materials,  so  materials  dependence  is  in  the 
product,  rpf(;  There  are  three  limits  for  Egg. 


li 

a 


i 

J 


il 


V. 


■-  '.^ '.' 


a.  Esc  ~  ^D*  This  occurs  when  E^  =  0  (no  applied 

field)  and  Eq  <  Eq.  For  materials  with  moderate  values  o£ 
dielectric  constant  (LiNb03,  LiTaO^,  QSO,  BGO;  ^  ^  50)  and 

typical  trap  densities  Ag  (N;^  ^10*^  cni"^),  is  smaller 
than  Eq  for  grating  periods  greater  than  about  0.5  um.  In 
materials  with  large  dielectric  constants  (BaTiO^,  KNbOj,  SBN) 
and  similar  trap  densities,  Eq  <  Eq  occurs  for  Ag  >1.5  to 

5  urn.  For  the  case  of  Eq  <  Eq,  the  steady  state  index  is 
proportional  to  rgff  since  Eq  does  not  depend  on  material 
parameters . 

b.  Eq  >  Eq  and  Eq  <  Eq.  Eg^,  is  proportional  to 

Eq ,  and  the  only  materials  dependence  of  the  steady  state  index 
is  through  feff/  as  in  case  (a). 

c.  Eq  or  Eq  >  Eq.  In  this  case,  Eg^  =  Eq  and  tlie 

materials  dependence  of  the  steady  state  index  is  in  the  product, 
r^ff  .  For  most  materials,  rgff  ~  ,  and  thus, 

the  steady  state  index,  depends  primarily  on  the  trap  number 
density,  N^.  Note  that  making  larger  in  order  to  increase 
the  steady  state  index  may  not  be  effective  if  it  changes  the 
dependence  from  case  (c)  to  case  (a)  or  (b).  Also,  increasing 
may  not  be  desirable  becruse  this  increases  the  response 
time,  as  discussed  below. 

2 .  Response  Time 

The  response  time  is  a  useful  figure  of  merit  for 
applications  in  which  available  energy  limits  the  illumination 
time,  or  in  which  the  grating  must  be  written  or  erased  in  a  set 
time  scale.  As  shown  in  Equation  (9),  the  response  time  is  given 
by  the  dielectric  relaxation  time  multiplied  by  a  factor  that 
takes  into  account  variations  in  the  conductiv'  /  due  to  the 
spatially  varying  electron  number  density.  In  order  to  under¬ 
stand  the  materials  dependence  of  the  response  time,  consider  the 
dielectric  relaxation  time,  ,  where  >  is  the 

C-  1  1 
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conductivity  given  by  J  =  euno.  Using  Equation  (13)  for  the 
electron  number  density,  we  obtain 
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(16) 


The  first  factor  in  brackets  is  a  material  figure  of  merit.  It 
is  essentially  a  constant  for  a  given  material,  although  there  is 
some  control  over  the  dielectric  constant,  e,  through  the  crystal 
i^rientation  and  control  over  the  photoionization  cross-section, 
s,  through  choice  of  operating  wavelength.  The  ratio  N;^/{Nq- 
N^) ,  the  inverse  of  the  reduction  ratio,  is  not  a  fundamental 
property  of  a  given  photoref ractive  material  and  can  be  varied  in 
non-destructive  processing  of  existing  crystals  (redn'^tion  or 
oxidation )  . 


3 .  Energy  Per  Area  to  Write  a  1%  Grating 

A  third  figure  of  merit,  which  is  widely  used,  is  the  energy 
per  unit  area  required  to  write  a  grating  of  1%  diffraction 
efficiency  in  a  crystal  1  mm  thick.  This  figure  of  merit  is 
useful  in  OOP  for  applications  where  a  large  efficiency  is  not 
required.  It  allows  comparison  of  fast  materials  with  small 
values  of  electro-optic  coefficient  with  slow  materials  having 
large  values  of  electro-optic  coefficient.  The  figure  of  merit 
has  the  disadvantage  that  some  materials  never  reach  a  1% 
diffraction  efficiency  in  1  mm  and  that  the  diffraction 
efficiency  depi?nds  on  the  energy  per  unit  area  only  for  times 
thaf  are  short  compared  to  the  response  time. 

The  energy  per  unit  area  to  write  a  1%  grating  is  derived  as 
follows.  The  diffraction  efficiency  is  given  by'^^ 

'1  =  exp(-'<ii/cos'))sin'(;^^^^~)  , 


(17) 


where  ^  is  the  absorption  coefficient,  is  the  length  of  the 
crystal,  0  is  half  the  Bragg  angle  and  X  is  the  vacuum 
wavelength.  If  we  ignore  the  oscillations  at  frequency  ■ 
(Equation  (10)),  the  index  change  can  be  written  as 

An  =  1/2  [l  -  exp(-t/i^)'  ,  (18) 


For  t  <<  T 


An  =  1/2  n.^r  ^^E  m  t/i 
b  eff  sc  e 


And  since  is  proportional  to  which  is  inversely 

proportional  to  Iq,  one  can  set  n  =  1%,  i  -  0.1  cm  and  solve 
for  lot,  which  is  the  energy  per  area  required  to  write  a  1% 
grating.  For  cases  in  which  Eg^  =  Eq  or  Egj,  =  E^, 

An  is  proportional  to  ^eff/^e'  that  materials  with 
large  electro-optic  coefficients  and  long  response  times  can  have 
similar  1%  energies  per  area  to  materials  with  small  electro-optic 
coefficients  and  short  response  times. 


PhO'  jrefractive  sensitivity 


The  photoref ract ive  sensitivity  is  defined  as  the  index 
change  per  absorbed  energy,  An/(  alQt )  . '"*’  The  photoref  ract  ive 
sensitivity  is  a  useful  figure  of  merit  because  it  tells  liow  well 
a  material  uses  a  given  amount  of  optical  energy.  Alternative!/, 
it  allows  comparison  of  materials  with  different  absorption 
coefficients  on  an  equal  basis.  Since  the  response  time,  i  , 
can  be  written  as  the  dielectric  relaxation  time  times  a  function 
of  the  grating  period  and  applied  field,  f(',EQ),  Equation  (18) 
can  be  used  to  write  the  photoref ractive  sensitivity  as 


An 

“I  ,t 


(■- 1  t 


since  n^  and  are  roughly  the  same  in  all  photoref ract ive 

mater  ials , ^  and  since  is  often  independent  of  material 

properties  (for  Eg  >  Eq  or  Eg  >  Eq) ,  the  photoref ractive 
sensitivity  is  mainly  determined  by  the  mobility,  u,  and  the 
recombination  time,  Tpj.  Large  mobilities  and  long 
recombination  time  are  desirable  because  these  properties  allow 
the  electrons  to  diffuse  or  drift  longer  distances  before  they 
recombine.  The  mobility  and  recombination  coefficient  are 
properties  of  a  given  crystal  that  appear  difficult  to  modify; 
however,  the  recombination  time  can  be  increased  by  decreasing 
the  number  density  of  traps  in  the  material,  which  is  sometimes 
possible  with  existing  crystals.  Note  that  decreasing  also 
lowers  Eg,  which  means  that  this  limiting  space  charge  field. 
Esc  ’  ®q'  will  ne  attained  for  smaller  values  of  Eq  or 
E,-)  . 


0,  Comparison  of  BSO  and  BaTiO ^ 

In  order  to  gain  a  better  understanding  of  the  figures  of 
merit  introduced  above,  we  have  calculated  them  in  detail  for  two 
attractive  photore f ract ive  materials,  BSO  and  BaTi03.  BSO  is 
known  to  be  a  Past  and  sensitive  material,  with  a  relatively 
small  electro-optic  coefficient.  On  the  other  hand,  BaTi03  has 
large  electro-optic  coefficients  (especially  r^j  =  r32)  and  is 
highly  effi'^ient,  but  it  is  known  to  respond  rather  slowly.  The 
materials  parameters  we  have  chosen  for  illustrating  these 
figures  of  merit  are  discussed  in  the  Section  F.  We  use  an 
optical  wavelength  of  A  =  0.5  urn,  extraordinary  polarization  • 
both  write  beams,  and  a  20°  angle  between  c  axis  and  g;  i»  i  ;  : 
wavevector.  The  la^5t  tvo  parameters  allow  CiJ  I  ci  I  a  ‘  i  •  • 

effective  <Jiolectric  constant,  ■  ,  and  th*^  eft  art  i  .■ 
coefficient.  The  value  of  20°  i.;  a  toj-r  -  ixi 


,  in  BaTiO,  an  I  a  common  oper  ii  in  : 
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The  steady  state  index  as  a  function  of  grating  period  for 
BSO  and  BaTi03  is  given  in  Figures  3  and  4  for  an  applied  field 
of  0  and  10  kV/cm.  The  maxima  in  the  Eq  =  0  curves  at  0.5  urn 
for  BSO  and  at  1.5  urn  for  BaTiOj  correspond  to  the  change  from 
the  Egj,  =■  Eq  limit  at  small  Ag  to  the  Eg^.  Eq 
limit  at  large  Ag.  por  Eq  =  10  kV/cm,  BSO  is  approaching  the 

Esc  “  Eq  limit  near  Ag  =  lo  urn;  in  BaTiOj,  with  Eq  = 

10  kV/cm  the  steady  state  index  increases  linearly  in  Ag , 

because  Eq  >  Eq  even  at  Ag  =  lo  iim.  Note  finally  that  the 
steady  state  index  in  BaTi03  is  20  to  50  times  larger  than  that 
in  BSO. 

Corresponding  curves  for  the  response  time  at  an  irradiance 
of  1  W/cm^  are  given  in  Figures  5  and  6.  Note  that  for  BSO  the 
response  time  decreases  with  grating  period  while  for  BaTiO  the 
response  time  increases  with  Ag.  Also  in  BSO  the  Eo  =  0 
response  time  is  shorter  than  at  Eq  =  10  kV/cm,  while  BaTiOj 
shows  the  opposite  behavior.  BSO  is  about  1000  times  faster  than 
BaTi03,  due  to  its  larger  absorption  coefficient,  smaller 
recombination  coefficient,  smaller  dielectric  constant,  and 
despite  its  smaller  mobility. 

Figure  7  shows  the  energy  per  unit  area  to  write  a  1% 
diffraction  efficiency  grating  in  a  1  mm  thick  crystal.  Since 
the  energy  per  area  is  proportional  to  the  response  time  divided 
by  the  steady  state  index,  the  requirement  for  both  BSO  and 
BaTi03  can  be  plotted  on  the  same  scale. 

Figure  8  shows  the  photoref ractive  sensitivity  or  index 
change  per  energy  per  volume  input  to  the  crystal.  Since  the 
absorption  coefficient  ratio  for  our  illustrative  parameters  for 
BaTi03  and  BSO  is  1.6/0. 3  =  5.3,  the  curves  in  Figure  8  are  a 
factor  of  5.3  closer  than  those  in  Figure  7.  For  small  grating 
periods  the  photoref ractive  sensitivity  of  BaTiOj  exceeds  BSO, 
while  above  A^  =  1  pm  BSO  is  more  sensitive. 
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Figure  5.  Response  time  of  BaTiO^  as  a  function 
of  grating  period  for  applied  fields, 
Eq  =  0,  10  kV/cm. 
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Response  time  of  BSO  as  a  function 
of  grating  period  for  applied  fields 
Eq  =  0,  10  kV/cm. 
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CONCLUSIONS  AND  OPTIMIZATION  OF  BSO  AND  BaTiO 


For  BSO  and  BaTi03  with  the  parameters  used  here,  BaTi03  has 
by  far  the  largest  steady  state  index,  while  BSO  has  a  much 
shorter  response  time  and  a  somewhat  smaller  energy  requirement 
for  a  1%  efficiency;  the  photoref ractive  sensitivities  are 
comparable. 

Other  factors  should  also  be  considered  in  comparing  these 
materials.  Certain  phase  conjugation  processes  and  image 
amplifiers  require  that  the  refractive  index  grating  be  shifted 
from  the  optical  interference  pattern  by  90*.  This  occurs  when 
Bo  =  0  or  when  Bq  <<  Bo .  Thus,  BSO  in  the  applied  field 
limit  cannot  be  used  as  an  image  amplifier  for  grating  periods 
larger  than  a  few  microns.  Some  of  the  phase  conjugation 
processes  require  that  the  steady-state  gain,  which  depends 
directly  on  the  steady-state  index,  exceed  a  threshold.^**  For 
these  purposes  BaTi03  may  be  above  threshold,  while  BSO  is  below 
threshold . 

Optimization  of  these  materials  for  better  performance  may 
be  possible  by  changing  the  crystal  parameters,  u,  N^,  N^/  s, 
and  itr.  Increasing  the  mobility  is  desirable  for  cases  in 
which  the  response  time  is  proportional  to  the  dielectric 
relaxation  time.  Increasing  the  number  density  of  filled  donors, 
Nq,  will  improve  both  the  speed  of  response  and  the  energy  per 
area  requirement  by  increasing  the  absorption  coefficient.  Nq 
can  be  increased  by  reducing  the  crystal;  i.e.,  increasing  N^ 
at  the  expense  of  N;^  or  by  increasing  the  number  density  of  the 
dopant.  The  limit  on  useful  increase  in  Np  is  set  by  the  total 
absorption.  For  typical  3  mm  thick  crystals  an  absorption 
coefficient  of  3  cm~ ^  gives  a  transmission  coefficient  of  0.41. 

If  one  increased  Nq  in  BaTi03  so  that  its  absorption  increased 
from  0.3  to  3  cm” ^ ,  the  response  time  and  energy  per  area 
requirement  would  decrease  a  factor  of  10. 

Optimization  of  performance  by  varying  the  number  density  of 
empty  traps  and  recombination  coefficients  is  more  complicated. 

On  the  simplest  level  one  would  want  to  increase  N;^  to  increase 


the  limiting  space  charge  field.  Eg,  and  one  would  want  to 
decrease  and  iCr  in  order  increase  the  recombination  time, 
increase  the  carrier  number  density  per  unit  irradiance,  and 
hence  decrease  the  dielectric  relaxation  time.  For  the  energy 
per  unit  area  for  a  1%  grating  and  the  photoref ractr ive 
sensitivity,  changing  would  in  this  simple  case  give  no 
effect,  and  decreasing  would  improve  the  response.  Closer 
inspection  of  Equation  (9)  for  the  response  time  yields  a 
different  answer.  Consider  four  cases  at  Ag  =  1  lun:  (1)  BSO 
with  Eo  =  0;  (2)  BSO  with  Eo  =  10  kV/cm;  (3)  BaTiOj  with 
Eo  =  0;  and  (4)  BaTi03  with  Eq  =  10  kV/cm.  In  case  (1) 


T  -  t  , .  T  /  T 

e  di  R'  D 


Since  is  proportional  to  changing  Yr  or  has  no 

effect  on  in  case  (2) 


=  T 


because  >>  1.  Since  xj  =  l/Y^no  and  n^  =  '-^I/YrN;^, 

Y|^  has  no  effect  in  this  limit,  and  decreasing  N;;^  decreases 
the  response  time.  In  case  (3) 


X  X  /x 
I  R 


Substitution  for  and  shows  that  is  proportional  to 


Yj^N  .  In  this  regime,  decreasing  will  decrease  the 


response  time  and  improve  the  1%  energy  requirement  and 
photoref ractive  sensitivity.  In  case  (4) 


X  X  2/t  2 

I  E'  R 


Here  x^  is  proportional  to  Y ^ ,  and  decreasing  and  Yr 
is  even  more  effective  in  improving  performance.  It  should  be 
noted  that  all  of  the  conclusions  in  this  paragraph  are  dependent 
on  grating  period  and  the  other  parameters  of  BSO  and  BaTiO,. 


MATERIALS  PARAMETERS  FOR  BSO  AND  BaTiO 


The  materials  parameters  needed  to  apply  the  charge 
transport  model  of  photorefractivity  are  given  in  Table  2.  The 
BSO  parameters  are  obtained  from  the  work  of  Peltier  and 
Micheron^^  and  the  review  by  Gunter, which  should  be  consulted 
for  the  original  references. 

Parameters  for  BaTiO 3  are  more  difficult  to  obtain.  We  use 
the  large  r^ 2  electro-optic  coefficient  for  our  comparisons 
because  many  of  the  most  exciting  new  processes  discovered  in 
BaTiOj  use  orientations  which  exploit  this  coefficient.  To  our 
knowledge  no  separate  measurement  of  photoionization  cross 
section  and  donor  number  density  has  been  made  in  BaTiO 3,  but  for 
the  calculations  performed  here  the  absorption  coefficient,  u  = 
0.3  cm“^  at  \  -  0.5  wm,^®  is  sufficient.  The  trap  number  density 
is  measured  from  grating  decay  times. 

The  recombination  coefficient  and  mobility  in  BaTiO 3  are 
difficult  to  obtain.  Mobilities  of  2.10"'^  to  0.5  cm^/Vs 
have  been  measured  in  BaTi03.  The  dielectric  relaxation  time 
measured  by  Feinberg  et  al.,^^  gives  the  ratio  Yr/u,  if  e  and  a 
are  known.  Also,  the  short  pulse  observations  in  BaTiO 3 
suggest  that  the  recombination  time  is  of  the  order  of  a 
nanosecond.  The  recombination  coefficient  of  5*10“®  is 
consistent  with  a  mobility  of  0.5  cm^/Vs  and  the  other 
observations.  Mobilities  as  low  as  2*10“**  cm^/Vs  seem  completely 
inconsistent  with  the  photoref ractive  observations.^®/®^  There 
is,  however,  still  one  problem  with  this  set  of  values.  The 
measured  value  of  u  =  0.5  cra^/Vs  is  an  electron  mobility,  while 
Feinberg  et  al.,^®  found  that  the  sign  of  the  moving  charges  in 
their  BaTi03  sample  was  positive,  and  in  general  one  expects  the 
hole  mobility  to  be  considerably  smaller  than  the  electron 
mobility. 


Table  2.  Materials  Parameters  for  BSO  and  BaTiO , 


Parameter 

Symbol 

Units 

BSO 

BaTiO , 

Background  index 

f’b 

2.5a 

2.43 

Dielectric  constant 

G 

- 

56^ 

4300^ 

168 

Electro-Optic 

coefficient 

nj 

pm/V 

1640e( r,2 )* 

80® ( r  j  ,  )* 

Donor  number  density 

Nd 

cm"  ^ 

10 

- 

Trap  number  density 

Na 

cm"  ^ 

10  1  uc 

2*10“’® 

Photocross-section 

S 

cm^ 

1.6*10" 

- + 

Recombination  coefficient 

cm  ^/s 

2*10"  ‘ 

5* 10"- 

Mobility 

u 

cin^/Vs 

0.03^ 

0.5, 

a  Gunter  [26! 

b  Peltier  and  Micheron  [23] 
c  Feinberg  et  al.  [15] 
d  Bursian  et  al.  [29] 
e  Landolt-Bornste in  [32] 

•  Undamped  value.  Clamped  values  820,  28 
+  s(Nd-Na)  =  a  =  0.3  [28]  at  X  = 


0.5  urn 
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The  volume  grating  solutions  developed  by  Kukhtarev  provide 
the  steady  state  and  step  function  response  for  the  space  charge 
density  and  space  charge  field,  but  are  not  directly  applicable 
for  pulsed  excitation.  In  this  study  we  present  a  model  for 
grating  formation  with  pulsed  excitation.  This  model  also 
applies  to  two-wave  mixing  and  four-wave  mixing  where 
nondepleted,  strong  pumps  and  a  weak  probe  wave  are  involved. 

Our  model  relies  on  the  £w  wing  assumptions; 

1.  The  pumps  are  cw. 

2.  The  incident  probe  is  a  pulse  of  duration,  fp,  and 
magnitude,  ^“p. 

3.  The  interaction  geometry  is  as  shown  in  Figure  D-1. 

For  DFWM,  all  three  input  waves  are  present.  For  two- 
wave  mixing,  we  shall  set 

4.  The  pump  waves  are  nondepleted  during  the  nonlinear 
interaction. 


The  response  of  the  material  to  the  three  input  waves  can  be 
described  in  terms  of  the  charge  transport  equations  of  Kukhtarev 
(as  reviewed  in  Appendices  B  and  C),  in  the  absence  of  applied  dc 
electric  fields; 
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=  V  KNn-So.)  - 

V  n  N, 

(5) 

cLinental  ■3eoraeti:y.  Ef,  and  Ep 
the  forward  pump,  backward  pump,  and 
e  waves,  respectively. 


with  The  quantities  n  and  are  the  electron 

and  ion  densities,  respectively,  q  is  the  charge  of  the  electron, 
J  is  the  electronic  current  density,  p  and  D  are  the  electron 
drift  and  diffusion  coefficients,  respectively.  Yj  and  Y^  are 
the  photo-ionization  and  recombination  constants,  respectively.  I 
is  the  intensity  of  the  optical  radiation,  and  c  is  the  static 
dielectric  constant. 

The  initial  conditions  for  our  problem  can  be  stated  as 
follows.  In  the  absence  of  the  probe  wave,  the  space  charge 
field  and  the  electronic  current  density  are  null.  These 
conditions  yield  the  steady  state  electron  and  ion  densities  in 
the  presence  of  the  pump  waves  only  and  are  given  by 


n^*^^  +  N^(0)  -  =  0 


(0)  _ 


N^.(0) 


r 

+ 


where  =  Yj  Iq  is  the  photo-ionization  rate  and  Iq  = 

Ipump  1^1^*  presence  of  the  probe  wave,  the 

charge  transport  equations  are  solved  in  the  perturbation  regime, 
consistent  with  the  assumption  of  non-depleted  pump  waves.  The 
first  order  equations  for  the  electron  and  ion  densities  are 


(vVe)"'”  *  =  ‘■p  r  <9) 


^  (VnK'”  *  ■’e 


where  =  qyn'  is  the  dielectric  relaxation  rate, 

r_  =  D  k,r-k^  is  the  electron  diffusion  rate,  r  =  y  n'  'is 
D  '  r  p  '  Nr 

the  ion  recombination  rate  and  r  =  y  is  the  electron 

e  r  + 

recombination  rate.  Equations  (9)  and  (10)  form  a  set  of  linear 
coupled  first  order  differential  equations  subject  to  the  initial 
condition  that  in  the  absence  of  the  probe  wave, 

,  n^^^  =  O  =  N  The  solution  for  Equation  (9)  and  (10) 

r'  ' 

can  be  obtained  by  means  of  Laplace  transform  techniques.  for 
our  case  we  have  chosen  the  time  dependence  of  the  probe  wave  to 
be  given  by 


{H(t)  -  H(t-Tp)}  ,  (11) 

where  Eq  is  a  constant  amplitude  and  H(t-t')  is  the  Heavyside 
function.  The  application  of  the  Laplace  transform  technique 
yields  the  following  solution  for  the  first  order  space  charge 
field,  Eg^,; 


E  =  -i 
sc 
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1,2 
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o  D  e  p  N 


B  =  (r^+r  j(r  +i’  j  +  r  r 
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(14) 

(15) 

(16) 


To  show  that  the  solution  has  the  appropriate  limits,  we  consider 
the  cases  at  t=0  and  at  t-*’".  At  t=0. 


S1S2 


^1-^2 


(i-  .  1_)  = 


(17) 


And  at  t ♦* , 


„(1)  »  (1)  a 

+  S  -  s 


1^2 


s^S2 


=  0 


(18) 


These  results  are  consistent  with  the  fact  that  in  the  absence  of 
Ep  the  first  order  densities  are  small.  Furthermore,  the 
steady  state  result  can  be  obtained  by  letting  and 

t>>Si~^  and  S2”^.  In  the  case. 


_r  r  ^ — - 

o  0  l_  ‘ 
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(19) 
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The  space  charge  field,  !isc»  gives  rise  to  a  spatial 
modulation  in  the  refractive  index  by  means  of  the  electro-optic 
effect.  The  scattering  of  the  forward  pump,  Cf,  in  the  case  of 
two-wave  mixing  (or  backward  pump,  in  the  case  of  four-wave 

mixing)  from  the  spatial  modulation  generates  the  signal  wave, 
fg.  The  signal  wave  will  have  the  same  time  dependence  as  the 
space  charge  field,  provided  that  the  nonlinear  medium  length,  L 
is  thin  enough  to  satisfy  the  condition. 
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where  c  is  the  speed  of  light  and  nQ  is  the  host  index  of 
refraction . 
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We  have  demonstrated  image  transmission  through  a  multimode  fiber  by  using  optical  phase  conjugation.  By 
using  the  wave-front  reversal  properties  of  degenerate  four-wave  mixing,  we  have  compensated  for  the  degradation 
of  the  image  that  is  due  to  modal  dispersion  in  a  fiber. 


We  have  demonstrated  the  transmission  of  images 
through  multimode  fibers  by  using  phase  conjugation.  ‘ 
The  transmission  of  spatial  information  through  fibers 
is  degraded  because  of  modal  dispersion.  To  com¬ 
pensate  for  this  degradation,  Yariv^  proposed  the  use 
of  the  wave-front-reversal  property  of  certain  nonlinear 
optical  techniques.  In  this  Letter  we  describe  the  Hrst 
reported  demonstration  of  image  transmission  and 
restoration  using  the  nonlinear  optical  process  of  de¬ 
generate  four -wave  mixing  (DPWM).  In  the  past,  other 
techniques  were  used  for  the  transmission  of  images 
through  fibers.  These  approaches  include  using  fibers 
of  a  precise  length,®-*  lenslike  index  variations,''-®  and 
holographic  filters.®-*® 

In  order  to  understand  image  degradation  by  trans¬ 
mission  in  an  optical  fiber  it  is  necessary  first  to  consider 
modal  dispersion.  Following  the  analysis  by  Yariv,®  any 
general  input  picture  field  can  be  expressed  as  the  sum 
of  the  bound  modes  and  unbound  modes  of  the  fiber: 

/oU.y.z  =  o,  f)  =  £  An„E„„(,x,y)expUa)t) 

m.n  *0 

radiation  modes.  (1) 

If  we  let  the  field  propagate  through  a  fiber  of  length  L 
and  assume  that  most  of  the  energy  is  carried  in  the 
bound  modes,  the  resultant  output  field  is  given  by 

MAT 

/lfx,_y,2*L,r)  ^  Afnn^mn  (x ,  y ) 

m,n*0 

XexpfiM  -  dmnLll.  (2) 

The  energy  that  is  lost  because  of  radiation  will  cause 
a  slight  reduction  in  resolution.  Also,  in  this  analysis, 
we  have  assumed  no  differential  modal  attenuation  or 
mode  mixing.  In  order  for  the  output  field  to  be  pro¬ 
portional  to  the  input  field  the  following  condition  must 
be  satisfied: 

0mnE  ~  2)r/„,  (3) 

where  dmn  is  the  propagation  constant  associated  with 
mode  m,  n  and  Im,n  is  an  integer.  Equation  (.3)  cannot 
be  satisfied  for  all  modes  simultaneously  for  an  arbitrary- 
value  of  L.  By  using  pha.se  conjugation  one  can 
reconstruct  the  original  input  field,  ll  the  beam 


transmitted  through  the  fiber  is  first  conjugated  and 
then  transmitted  through  an  identical  fiber,  the  original 
input  image  will  be  formed  at  the  output  of  the  second 
fiber.  This  can  be  seen  by  taking  the  conjugate  of  Eq. 
(2)  and  propagating  the  field  through  an  identical  fiber 
of  length  L  (or  back  through  the  same  fiber): 

f2(x,y,z  =  L,t)  =  ^  Amn*E„„(x,y) 

fr.,n»0 

X  explKwt  +  HmnD],  (4) 

The  output  after  traversing  an  identical  fiber  or  retra¬ 
versing  the  same  fiber  is 

f:t{x,y,z  =  o,t)  =  fAx.y.z  =  2L,t) 

“  ^  ^mn  *Emn(x,yhxp{ia)t),  (5) 

m.n  *0 

which  is  seen  to  be  proportional  to  the  input  field  in  Eq. 

(1). 

In  the  present  experiments  a  single  fiber  was  used. 
The  purpose  of  these  experiments  was  to  show  that 
pha.se  conjugation  can  compensate  for  modal  distortions 
introduced  by  a  fiber. 

The  demonstration  was  accomplished  by  utilizing  a 
multimode  (85-Mm  diameter,  1.75-m-long)  step-index 
fiber.  The  output  of  a  krypton  laser  (6471  A)  was  used 
to  backilluminate  a  portion  of  a  U.S.  Air  Force  resolu¬ 
tion  chart  that  was  imaged  at  the  input  end  of  the  fiber. 
The  fiber  output  was  directed  into  a  crystal  of  barium 
titanate  (BaTiOn)  in  which  DFWM  was  used  to  generate 
the  phase-conjugate  signal.  The  phase-conjugate  signal 
then  retraversed  the  same  fiber,  resulting  in  a  recon¬ 
structed  image  of  the  resolution  chart. 

The  geometry  employed  for  the  four-wave  mixing 
experiments  in  this  work  is  similar  to  the  high-reflec¬ 
tivity  configuration  for  BaTiO;,  of  Feinberg  and  Hell- 
warth."  It  consists  of  two  counterpropagating  pumps 
E/  and  Ei,  and  a  probe  wave  E,,  (which  is  the  lieam  with 
the  image  information  impo.sed  upon  it),  incident  upon 
the  nonlinear  medium.  'Phese  beams  interact  to  gen¬ 
erate  the  phase -conjugate  signal  E',.  A  schematic  of  t  he 
experiment  is  shown  as  Fig.  1.  The  krypton-ion  laser 
operates  in  a  single  longitudinal  mode,  giving  coherence 
lengths  of  the  order  of  I(t  in.  Because  of  crvstal  orien- 
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Fig.  1.  Schematic  of  experiment  for  image  transmission 
through  an  optical  fiber. 


tation,  p  polarization  of  the  input  beams  is  required  in 
our  experiments.  This  polarization  is  desired  to  couple 
into  the  large  r 42  electro-optic  coefficient  of  BaTiO.-i  (>-42 
=  820  X  10“''^  mV“').  Therefore  a  X/2  plate  and  a 
Faraday  rotator  were  used  to  rotate  the  plane  of  polar¬ 
ization  of  the  krypton  laser.  The  Faraday  rotator  was 
also  used  in  conjunction  with  a  polarizer  to  prevent  any 
return  beams  from  entering  the  laser  and  spoiling  the 
single-mode  operation.  Beam  splitter  BSl  picks  off  the 
pump  beams.  An/  =  100-cm  lens  (LI)  was  used  to  re¬ 
duce  the  pump  beam  in  the  crystal  to  a  size  that  opti¬ 
mized  the  interaction  length  and  probe  overlap.  Beam 
splitter  BS2  generated  a  secondary  probe  signal  that  was 
used  to  optimize  the  crystal  orientation.  Beam  splitter 
BS3  was  used  to  generate  the  individual  counterprop- 
agating  pumps.  The  pumps  typically  have  a  power 
ratio  of  2:1  (/«  =  7.1  mW,  //  =  3.2  mW).  The 
throughput  of  beam  splitter  BSl  was  the  probe  beam 
propagated  through  the  fiber.  This  probe  beam  was 
expanded  by  a  beam-expanding  telescope  to  illuminate 
the  Air  Force  resolution  chart  Following  the  resolution 
chart  is  a  beam  splitter  (BS5),  which  is  used  to  pick  off 
the  return  signal  that  forms  the  reconstructed  image  on 
a  film  plane  at  image  plane  1.  An  //I,  7.6-cm  focal- 
length  relay  lens  (L2)  was  used  to  capture  the  light 
diffracted  by  the  patterns  of  the  resolution  chart. 
Lenses  L3  and  L4  are  7X  objectives  with  0.20  N.A.  and 
are  used  to  couple  the  light  into  and  out  of  the  Corning 
fiber.  With  an  N.A.  of  0. 184,  theoretically  on  the  order 
of  10^  modes  propagate  at  the  wavelength  being  used. 
The  second  microscope  objective  (L4)  quasi-collimates 
the  output  beam  from  the  fiber.  A  5-cm  focal-length 
lens  (L5)  is  used  to  reduce  the  probe  beam  to  a  diameter 
smaller  than  the  diameter  of  the  pump  beams  inside  the 
crystal.  TVpicady*  angle  between  the  forward  pump 
and  the  probe  inside  the  crystal  is  7°,  and  the  grating 
vector  makes  an  angle  of  18°  with  respect  to  the  c  axis, 
a  condition  that  yields  the  highest  reflectivity.  By  using 
this  geometry  with  the  crystal  in  air,  we  have  obtained 
reflectivities  in  excess  of  150%.  These  large  returns 
have  enabled  us  to  reconstruct  the  image  and  also  to 
show  image  amplification. 

The  mutual-coherence  requirements  between  the 
prolie  (i?p)  and  the  pump  (E/)  are  important  in  the 
present  experiment.  It  has  been  well  established  that 
a  holographic  analogy  exists  in  four-wave  mixing.  In 
particular,  one  can  think  of  the  interference  of  the  probe 


Ep  and  the  pump  Ef  as  generating  a  volume  hologram 
in  the  BaTiOs,  which  is  then  read  out  by  pump  Eb  to 
generate  the  conjugate  wave  This  would  suggest 
that,  for  the  experimental  setup  used,  the  coherence 
length  of  the  laser  must  be  at  least  as  long  as  the  path- 
length  difference  between  the  probe  Ep  and  the  pump 
Ef  (which,  for  our  experiments,  is  ~4  m). 

In  fact,  the  above  mutual-coherence  requirement 
strictly  applies  only  if  the  time  response  of  the  nonlinear 
medium  is  infinite  and  approximately  applies  when  one 
uses  slowly  responding  materials,  sqph  as  BaTiOs  (100 
msec).  However,  if  nonlinear  materials  that  respond 
instantaneously  are  used,  then  this  requirement  is  al¬ 
tered.  To  understand  this,  recall  that  the  product  of 
the  response  time  t  and  the  frequency  bandwidth  Ai> 
(or  the  pump-probe  frequency  difference)  is  unity. 
Thus,  when  BaTiOs  is  used,  Ai>  ~  10  Hz  is  required.  If 
a  material,  such  as  sodium  vapor,  with  nanosecond  re¬ 
sponse  times  is  used,  a  bandwidth  of  a  few  gigahertz  is 
allowed  (more  precisely,  Ai>  is  given  by  the  generalized 
Rabi  frequency’^).  This  means  that  interference  will 
occur  for  any  length  fiber  as  long  as  the  frequency  dif¬ 
ference  between  the  probe  and  the  pump  is  less  than  a 
few  gigahertz,  which  is  easy  to  achieve  by  using  stabi¬ 
lized  dye  lasers. 

The  results  of  this  experiment  are  shown  in  Fig.  2. 
The  first  photograph  shows  the  input  image  of  a  portion 
of  a  U.S.  Air  Force  resolution  chart.  The  second  pho¬ 
tograph  shows  the  output  at  the  end  of  the  fiber.  The 
output  is  typical  of  a  highly  multimode  fiber  in  which 
a  large  number  of  modes  are  excited.  From  the  output, 
it  is  impossible  to  discern  that  any  spatial  information 
has  been  transmitted  through  the  fiber.  The  last 
photograph  shows  the  image  that  has  been  recon¬ 
structed  after  double  passing  the  fiber  (approximately 
25-fitn  spacings  are  shown). 

An  interesting  issue  related  to  the  fidelity  of  the 
image  reconstruction  is  the  polarization.  Typical 
commercially  available  fibers  do  not  preserve  polar¬ 
ization;  therefore,  when  linearly  polarized  light  is  fed 
into  a  fiber,  the  output  is  randomly  polarized.  In  our 
DFWM  experiments,  only  the  p -polarization  compo¬ 
nent  is  conjugated.  This  polarization  is  scrambled  on 
retraversing  the  fiber.  Because  the  input  image  is  re¬ 
constructed.  we  speculate  that  the  information  carried 
by  each  polarization  component  is  similar. 

Another  test  of  the  fidelity  of  the  phase-conjugation 
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Fig.  2,  Den:()nstration  of  image  transmission  through  an 
optical  fiber:  a,  input  image  resolution  pattern;  h,  output 
after  a  single  pass  through  the  fiber;  c,  reconstruction  of  res¬ 
olution  pattern  after  output  has  been  phase  conjugated  and 
retraverses  fiber. 
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Fig.  3.  Photographs  of  the  far-field  intensity  distribution  forta,  input  beam;  b,  single  pass  through  Tiber;  c,  corrected  beam. 
Oscilloscope  traces  d,  e,  and  f  are  the  corresponding  cross  sections. 


process  is  to  measure  its  ability  to  reproduce  a  Gaussian 
input  wave.  In  these  experiments  the  beam  passes 
through  the  fiber,  is  conjugated,  retraverses  the  fiber, 
and  is  examined  in  the  far  Held.  The  beam  to  be  stud¬ 
ied  was  monitored  at  the  back  focal  plane  of  a  1-m 
focal-length  lens.  This  image  was  magnified  and  then 
projected  onto  the  vidicon  element  of  a  television 
camera.  The  video  signal  was  processed  by  an  image 
analyzer  that  produces  two  outputs.  The  first  is  a 
display  of  intensity  versus  x-y  position.  The  video 
signal  is  level  sliced,  and  each  level  is  assigned  a  color. 
This  output  is  displayed  on  a  color  television  screen. 
The  second  format  is  an  intensity  profile.  A  cursor 
permits  one  to  display  a  cross  section  of  the  far-field 
pattern  on  an  oscilloscope.  The  results  of  the  far-field 
experiments  are  shown  in  Fig.  3.  The  input  beam  is 
monitored  at  image  plane  2  and  is  typical  of  an  apodized 
Gaussian  beam  (Fig.  3a).  The  output  from  the  fiber, 
monitored  at  image  plane  3,  is  shown  in  Fig.  3b.  It  is 
seen  to  be  very  irregular  in  intensity,  indicating  the  se¬ 
vere  aberrations  introduced  by  the  fiber.  This  output 
is  then  conjugated  and  retraverses  the  fiber.  The  re¬ 
constructed  bieam  is  monitored  at  image  plane  1  and  is 
seen  very  nearly  to  reconstruct  the  initial  input  beam 
(Fig.  3c).  The  corresponding  intensity  cross  sections 
3d,  3e,  and  3f  were  taken  at  the  position  of  the  horizontal 
line  in  Figs.  3a,  3b,  and  3c. 

In  conclusion,  we  have  demonstrated  the  ability  of 
phase  conjugation  to  compensate  for  modal  dispersion 
in  optical  fibers.  These  experiments  were  done  for  a 
static  two-dimensional  image,  but  the  technique  has  a 
wider  range  of  applicability.  In  particular,  the  real-time 
correction  capability  for  modal  dispersion  is  limited  only 
by  the  time  response  of  the  phase-conjugate  medium. 
Many  materials  have  the  nanosecond  response  time 


exemplified  by  sodium  vapor.  These  fast  response 
times  not  only  allow  for  transmission  at  a  high  data  rate 
but  also  correct  for  dynamic  thermal  and  mechanical 
stresses  in  the  fiber.  In  addition,  since  gain  can  be 
achieved  by  using  DFWM,  the  phase  conjugator  could 
be  used  as  an  optical  amplifier  in  a  transmission  line. 

We  acknowledge  valuable  discussions  with  D.  G. 
Steel,  A.  Yariv,  M.  Rourke,  and  C.  Giuliano.  We  thank 
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ported  in  part  by  Rome  Air  Development  Center  under 
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